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a  acceleration  vector 

A  any  material  property 

ae  earth»s  gravitational  acceleration 

B  any  external  influence 

c  specific  heat 
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D(z)  dimension  of  physical  quantity  z 
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z  )  dimensional  relationship  between  z  ,  z  ,  and  z  . 
3  i  a  s 

density  function  for  the  physical  quantity  z. 

modulus  of  elasticity  in  the  linear  elastic  range, 
base  of  natural  logarithm 

correlation  function  for  the  physical  quantity  z 
applied  force  (Appendix  A) 
specific  mass 
body  force 

coefficient  of  linear  thermal  expansion 
heat 
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length 

load  per  influence  volume 
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mass 

ultimate  moment  of  section  at  complete 
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applied  surface  stress  (section  2.7) 
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subscript  used  to  denote  prototype  properties 

any  property  of  a  structure 

distance  in  system  of  spherical  coordinates 

rate  of  change  of  slope,  or  the  curvature  of  the 
elastic  curve 

standard  deviation  for  physical  quantity  Z 
time 

temperature 
displacement  vector 
velocity 
influence  volume 
shearing  component  of  load 
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specific  gravity- 
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coefficient  of  thermal  conductivity 

gravitational  constant 

vertical  deflection 
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deflection  along  N 
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rupture 
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friction  coefficient 
Poisson  Ratio 
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NOTATION  (Cont'd.) 


*  angle  in  system  of  spherical  coordinates 

( in  Chapter  2) 

♦  half  the  Joint  rotation  in  the  knee 

(in  Chapters  1,  3  and  ij.) 

$  angle  in  system  of  spherical  coordinates 

( in  Chapter  2) 

rate  of  change  of  slope  of  the  elastic  curve 
O'  stress 


O' 

2 

r 

max 


lower  yield  stress  of  a  material 
upper  yield  stress  of  a  material 
dynamic  yield  stress 
major  principal  stress 
minor  principal  stress 
maximum  shear  stress 
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ABSTRACT 


The  objective  of  this  research  work  is  the  inves¬ 
tigation  of  materials  and  techniques  that  can  be  used  in 
model  studies  of  structures  under  blast-type  loads.  The 
materials  are  selected  on  the  basis  of  theoretical  con¬ 
siderations,  and  their  usefulness  is  demonstrated  by  means 
of  actual  model  structure  tests  which  are  compared  to 
prototype  tests. 

The  theory  of  structural  models  is  reviewed  and 
extended  to  a  very  general  form,  with  emphasis  on  the 
basic  differences  in  the  approach  suitable  for  analysis 
by  means  of  the  mathematical  model  and  the  physical  model 
respectively.  This  extension  of  the  theory  proves  to  be 
necessary  when  the  large  number  of  factors  affecting  a 
dynamic  model  study  in  the  plastic  range  are  to  be  taken 
into  account  in  a  systematic  manner.  A  description  of 
the  general  properties  of  a  structure  is  followed  by 
a  dimensional  analysis  leading  to  a  comparison  between 
model  and  prototype.  The  criteria  for  compatibility  between 
model  and  prototype  materials  are  derived  in  detail. 
Modelling  is  discussed  as  a  random  process,  and  the  theory 
is  applied  to  massive  structures  in  particular. 

It  is  shown  that  plastic  deformations  during 
dynamic  response  impose  severe  restrictions  on  the  equation 
of  state  of  the  model  material  which  is  to  be  compatible 
with  steel  as  a  prototype  material.  The  properites  of  phos¬ 
phor  bronze  and  ethyl  cellulose  and  their  usefulness  as 
modelling  materials  for  steel  structures  under  dynamic  loads 
are  discussed.  Manufacturing  techniques  for  phosphor  bronze 
and  the  loading  devices  used  in  the  studies  are  described 
in  detail.  The  results  of  model  studies  of  beams,  joints 
and  portal  frames  under  static  loads  are  compared  to  results 


obtained  from  tests  on  full  scale  structures,  and  a  tech¬ 
nique  of  testing  models  under  dynamic  loads  is  described 
Since  model  and  prototype  results  generally  show 
good  agreement,  it  is  concluded  that  the  metal  phosphor 
bronze  can  be  used  as  a  modelling  material  in  small  scale 
studies  of  welded  steel  structures  in  all  cases  where  the 
model  and  prototype  surface  loads  do  not  have  to  be  of 
equal  intensity.  This  will,  for  instance,  be  the  case 
where  the  blast  loading  is  mainly  of  the  drag  type.  The 
use  of  plastics  in  the  modeling  of  steel  structures 
presents  problems  which  cannot  be  overcome  without  a 
reliable  error  analysis  of  the  model  since  in  this  case 
the  model  will  be  distorted  above  the  proportional  limit. 
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CHAPTER  1 


INTRODUCTION 


The  use  of  physical  models  in  the  study  of  engineer- 
lng  problems  has  become  prevalent  In  the  past  decade  (1 ,2,30  . 
Two  general  classes  of  problems  lend  themselves  particularly 
to  a  model  solution.  First,  In  problems  where  the  number  of 
variables  or  the  complexity  of  the  situation  make  the  analyt¬ 
ical  solution  difficult  or  particularly  cumbersome  and  second, 
in  areas  where  the  general  laws  governing  the  phenomenon 
have  not  yet  been  sufficiently  developed  and  a  general  theory 
is  lacking.  The  study  of  the  dynamic  response  of  structures 
is  one  area  of  structural  mechanics  where  the  modeling 
method  of  analysis  may  be  used  to  advantage.  It  is  the  aim 
of  this  report  to  discuss  the  problems  of  using  structural 
models  to  study  the  dynamic  behavior  of  structures  in 
general  and  to  present  specific  examples  in  the  choice  of 
modeling  materials  and  the  modeling  techniques  developed 
for  their  use. 

1.1  THE  STRUCTURAL  MODEL 

A  structural  model  is  a  device,  which  being  a 
reproduction  of  the  prototype  following  certain  rules,  can 
be  used  to  find  particular  unknown  quantities  about  the 
prototype.  Two  large  classes  of  models*  emerge  from  this 
definition.  One  encompasses  all  those  model  studies  dealing 
with  elastic  problems.  In  this  case  both  model  and  prototype 
are  within  the  elastic  range  and  the  results  can  be  compared 
with  those  obtained  analytically  If  such  methods  exist.  A 
second  group  representing  a  more  recent  approach  to  the  use 
of  models  is  in  closer  agreement  with  present  day  concepts 
of  structural  analysis  and  design^.  Here,  the  model  Is 

Numbers  in  parenthesis  indicate  the  particular  reference 

cited  as  listed  at  the  end  of  the  report. 
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a  faithful  reproduction  of  the  prototype  following  the  laws 
of  similitude  and  boundary  conditions  from  which  the  ultimate 
behavior  of  the  prototype  can  be  predicted  as  well  as  the 
factor  of  safety  which  is  defined  as  the  ratio  of  ultimate 
load  to  the  working  load  acting  on  the  structure.  Although 
the  second  approach  by  the  nature  of  Its  scope  places  more 
restrictions  on  the  model,  it  is  the  only  one  feasible  in 
certain  problems  where  plastic  deformations  are  tolerated. 

Such  problems  arise  in  the  study  of  the  dynamic  response 
of  structures  where  energy  absorption  criteria  allow  the 
use  of  some  plastic  deformations  in  design.  Two  basic 
types  of  models  arise  in  our  study  of  structural  problems. 

One  is  the  stress  model  or  a  model  which  falls  because  of 
rupture  of  the  material  at  local  points,  and  the  buckling 
model  which  becomes  unstable  at  certain  critical  loads  and 
results  in  large  deflection.  Our  main  concern  here  is  the 
stress  type  of  model  and  we  shall  attempt  to  investigate 
Its  behavior. 

1.2  A  HISTORICAL  SURVEY  OP  MODEL  STUDIES 

Physical  models  have  been  used  to  some  extent  In 
structural  analysis  and  design  for  many  years.  In  general 
the  principal  external  forces  involved  in  these  studies 
have  been  static  in  nature.  Further,  only  the  elastic 
behavior  of  the  model,  hence  the  prototype,  was  of  Interest. 
Model  studies  to  verify  the  analytical  results  of  compli¬ 
cated  structural  problems  are  recorded  in  the  literature^0* 

In  most  of  these  cases  the  Importance  of  the  project  made 
It  feasible  to  make  a  model  study  but  in  all  these  cases 
the  model  was  used  as  a  check.  It  was  inconoelvalbe  then 
to  base  the  design  on  the  results  of  a  thorough  model  study. 

As  teaching  aids  in  Indeterminate  structural  analysis,  models 
have  found  their  place  in  the  laboratories  of  many  schools^^  . 
These,  however,  have  been  predominantly  of  the  first  class 
or  elastic  models. 
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More  recent  model  studies  especially  in  the  field 
of  arch  daws  hsve  been  done  in  Europe  J  '  .  Techniques 

have  been  developed  and  refined  to  the  point  where  the  model 
study  plays  a  major  role  in  the  design  of  arch  and  domed 
dams.  Besides  the  extensive  work  on  dams  however,  model 
studies  of  complicated  buildings  and  earthquake  effects  as 
well  as  pioneering  work  in  model  studies  beyond  the  elastic 
limit  have  also  been  done^^'.  All  being  said,  the  European 
model  studies  have  set  a  good  precedent  and  have  prompted 
many  engineers  to  base  their  designs  on  a  model  study. 

In  the  United  States  great  Impetus  has  been  given 
to  model  studies  in  the  postwar  years.  This  revival  stems 
from  two  main  sources.  Architectural  taste  for  more  compli¬ 
cated  but  pleasing  structural  forms  and  the  intensive  efforts 
to  develop  techniques  for  the  design  of  structures  and 
structural  elements  to  resist  dynamic  loads  have  introduced 
some  very  interesting  but  difficult  problems.  The  models  in 
studies  dealing  with  dynamic  response  have  been  loaded 
to  failure  in  shock  tubes,  by  Impulsive  loading  machines 
and  in  some  cases  by  high  explosive  detonations^  8. 18 *26 # 
In  the  majority  of  cases  reinforced  concrete  structures 
were  modeled  by  the  steel  wire-mortar  technique.  Dynamic 
tests  on  model  or  prototype  steel  structures  are  far  less. 
Many  problems  are  encountered  in  such  studies  but  the  results 
up  to  the  present  (1962)  seem  encouraging  and  it  is  antici¬ 
pated  that  with  future  research  enough  data  will  be  collected 
to  make  the  modeling  technique  an  indlspensible  tool  at  the 
disposal  of  the  structural  engineer. 

1.3  LIMITATIONS  OP  PRESENT  DAY  MODEL  STUDIES 

Although  the  general  theory  of  structural  models 
under  the  most  general  conditions  has  been  known  for  a  long 
time  its  application  to  specific  cases  has  been  limited. 

The  reasons  are  many  but  among  them  one  stands  out  as  a 
limitation  which  must  be  overcome  before  the  modeling  tech¬ 
nique  can  be  applied  with  confidence  to  non-linear  problems 
in  structures  and  to  distorted  models.  This  limitation  is 
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the  leek  of  a  comprehensive  method  of  evaluating  errors  In 
models  which  may  be  "distorted".  Thus  far  model  studies 
are  limited  to  elastic  cases  and  to  a  few  known  rheologically 
compatible  model  materials.  In  many  cases  the  same  material 
Is  used  In  the  model  (to  be  on  the  safe  side)  and  external 
means  are  employed  to  satisfy  similitude  thus  introducing 
errors  of  a  different  nature.  An  error  analysis  of  models 
would  be  a  great  contribution  to  model  theory  and  would 
open  new  horizons  in  the  applications  of  the  modeling  tech- 
nlque. 

More  work  is  needed  to  establish  criteria  on  the 
scale  limits  of  model  studies.  It  is  obvious  that  the  smaller 
the  model  the  bigger  the  chance  of  increasing  the  errors  in 
the  results.  Thus  if  there  is  a  correlation  between  percent¬ 
age  error  and  scale  for  a  particular  class  of  problems  the 
experimenter  would  be  able  to  pick  a  scale  which  will  be 
within  the  tolerance  of  his  overall  allowable  error. 

Another  area  which  needs  more  study  is  the  need  of 
a  method  of  making  the  model  and  prototype  materials  com¬ 
patible  while  satisfying  the  other  similitude  requirements. 

It  is  hoped  that  some  light  will  be  thrown  upon  this  prob¬ 
lem  in  the  sections  that  follow. 


CHAPTER  2 


A  GENERAL  THEORY  OF  STRUCTURAL  MODELS 


The  basic  aim  of  the  analysis  of  a  structure  is 
the  prediction  of  its  behavior  at  a  time  different  from 
that  at  which  initial  conditions  are  specified.  For  the 
purpose  of  the  investigation,  the  nature  of  the  structure 
may  be  described  either  abstractly  by  means  of  a  mathematical 
model  or  concretely  by  means  of  a  physical  model,  frequently 
called  "the  structural  model".  Each  of  the  two  methods 
has  its  own  specific  advantages.  The  mathematical  method 
tends  to  use  a  model  in  which  some  of  the  physical  events 
occurring  in  the  actual  structure  are  either  simplified  or 
neglected.  As  a  result,  a  speedy  analysis  is  possible  with* 
out  excessive  loss  in  accuracy.  The  physical  model,  on  ■ 
the  other  hand,  permits  the  investigator  to  account  accurately 
for  a  large  number  of  physical  events  and  variables  even 
if  they  have  complex  space-time  distributions. 

This  important  property  of  the  structural  model 
can  be  used  to  full  advantage  only  if  the  experiments  are 
based  on  a  very  general  theory.  Indeed,  direct  application 
of  formulations  appropriate  for  use  in  current  methods  of 
mathematical  analysis  will  unnecessarily  limit  the  scope  of 
structural  model  investigations.  This  reasoning  has  led 
to  the  unusual  formulations  encountered  in  some  of  the 
later  sections. 

The  theory  of  structural  models  forms  the  basis 
for  the  planning,  execution  and  interpretation  of  laboratory 
studies  on  a  physical  model  which  permit  the  prediction 
of  the  behavior  of  a  prototype  structure  under  static  or 
dynamic  conditions.  It  establishes  the  rules  according  to 
which  the  geometry,  material  properties,  initial  conditions 
and  boundary  conditions  of  the  model  and  of  the  prototype 
have  to  be  related  so  that  the  behavior  of  the  one  can  be 
expressed  as  a  function  of  the  behavior  of  the  other.  The 
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degree  of  perfection  with  which  these  relationships  can  be 
maintained  during  the  actual  experiment,  and  the  effect  of 
any  possible  deviation  on  the  accuracy  of  the  prediction  of 
the  prototype  behavior,  are  discussed  briefly  in  Section  2.$. 

The  physical  events  considered  in  this  theory  are 
described  in  Sections  2.1  and  2.3.  If  an  event  occurs  in 
either  model  or  prototype  which  involves  physical  laws  not 
considered  in  the  theory  present  ed  here,  the  model  studies 
will  lead  to  erroneous  results.  So,  for  instance,  electro¬ 
magnetic  effects  are  not  taken  into  account  since  they  play 
no  role  in  the  vast  majority  of  structural  problems. 

2.1  DESCRIPTION  OF  THE  MODEL  AND  PROTOTYPE  STRUCTURES 

Before  a  one  to  one  correspondence  between  model 
and  prototype  can  be  established,  the  properties  of  and 
external  influences  on  each  of  them  individually  must  be 
specified.  The  following  physical  quantities  in  most  cases 
fully  describe  the  nature  of  a  structural  problem: 

2.1.1  Geometric  Properties.  If  time  is  not  considered 
to  be  a  basic  dimension,  a  point  in  a  structure  is  fully 
described  by  its  three  space  coordinates.  A  general  theory 
must,  however,  take  the  variation  of  geometry  with  time 
into  account,  so  that  time  must  be  regarded  as  a  coordinate 

in  the  same  manner  as  the  space  coordinates. 

* 

If  no  mass  is  added  to  or  removed  from  the  struc¬ 
ture  during  the  time  interval  in  which  its  behavior  is  to 
be  studied,  its  geometry  needs  to  be  specified  at  one  time 
t0  only.  The  geometry  at  any  other  time  t  will  then  be 
uniquely  determined  by  the  displacement  vector  uv defined  later. 
But  if,  for  instance,  the  erection  period  of  the  structure  is 
to  be  Included  in  the  model  study,  a  complete  time  history 
of  the  geometry  during  that  interval  must  be  given. 

The  position  of  any  mass  particle  of  a  struc¬ 
ture  in  space  may  be  described  in  any  one  of  the  many  well- 
known  coordinate  systems.  It  is  thus  necessary  to  prove 
that  the  theory  loses  none  of  its  generality  if  it  is 
restricted  to  a  single  system  of  coordinates.  Assume, 
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for  instance,  that  spherical  coordinates  are  chosen,  so  that 
the  position  of  point  P  at  time  t  is 

P  5  p(r,  ♦,  9,  t)  (2.1.1) 

In  all  other  coordinate  systems,  P  would  also  be  fixed  by 
at  least  one  length  and,  in  some  cases,  dimensionless  ratios 
such  as  angles.  Thus  no  new  physical  variables  will  be  intro¬ 
duced  by  a  coordinate  transformation,  and  the  theory  may  be 
limited  to  the  system  of  equation  (2.1.1). 

The  extent  of  the  structure  is  defined  by  its 
enclosing  surfaces  which  can  be  expressed  by  equations  of 
the  form 

Vt  (r,  ♦,  f )  =  0  (2.1.2) 

The  number  i  of  these  equations  varies  from  one  upward  and 
is  equal  to  the  number  of  surfaces  by  which  the  structure 

is  bounded.  As  time  varies,  their  number  may  increase  or 

decrease  and  the  region  of  validity  of  existing  functions 
may  change.  Thus  the  limits  of  any  two  of  the  three  space 
coordinates  must  be  given  in  terms  of  the  fourth  coordinate, 
time: 

rj(t)  £  rt  *  rk(t> 

♦j(t)  «  ♦1  *  *k(t)  (2.1.3) 

where  i  is  the  same  index  as  in  equation  (2.1.2).  The  limits 
of  the  third  space  coordinate,  f ,  can  be  calculated  from 
equation  2.1.2  and  the  geometry  of  the  entire  structure  is 
now  fully  described. 

2.1.2  Material  Properties,  A  structure  will,  in  the 
most  general  case,  consist  of  several  materials,  each  of  which 
has  its  own  particular  physical  properties  that  may,  or  may 
not,  change  with  time  and  temperature.  For  the  purpose  of 
this  theory  it  is,  however,  more  convenient  though  still 
perfectly  general  to  treat  the  structure  as  if  it  consisted 
of  a  single  material  whose  properties  change  from  point  to 
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point  in  space  as  well  as  with  time  and  temperature.  The 
intrinsic  properties  of  the  material  are  then  given  by  the 
functions: 


0 


Specific  Mass: 

Specific  Heat: 

Coefficient  of  Thermal 
Conductivity: 

Coefficient  of  Linear 
Thermal  Expansion: 


g  =  g(r,  ♦ ,  ?  ,  t,  T) 
c  =  c(r, e ,  f ,  t,  T) 

x  =  x(r,  &  ,  f  ,  t,  T) 

h  =  h(r,  ft  ,  <f  ,  t,  T) 


(2.1. lj.) 

A  material  is  not  fully  specified  by  its  intrinsic  properties. 
Those  of  the  material  properties  which  are  dependent  on 
external  Influences  are  interrelated  by  the  equation  of 
state  of  the  material,  which  will  be  presented  in  Section  2.3. 

2.1.3  Initial  Conditions.  When  the  geometry  of  the 
structure  is  specified  at  time  t0>  the  material  may  already 
be  subject  to  initial  stresses  due  to  prestressing  or  other 
processes  that  have  taken  place  before  time  t0.  There  will 
also  be  an  initial  temperature  distribution  so  that  it  is 
necessary  to  specify 


Initial  Stress:  5b  =  ir(r,«,f  ,  to) 

Initial  Temperature  :  T0  3  T(r, * , f ,  t0) 

(2.1.5) 

If  the  history  of  a  structure  is  known  from  the  beginning 
of  its  erection  period,  end  all  events  (  such  as  shrinkage 
of  concrete)  during  and  immediately  following  this  period 
are  taken  into  account  in  the  analysis,  no  specific  initial 
conditions  need  to  be  specified.  In  some  cases  where  the 
molecular  processes  occurring  in  a  prototype  (e.g.,  during 
curing)  cannot  be  reproduced  correctly  in  the  physical  model, 
the  specification  of  initial  stresses  may  become  mandatory 
even  if  the  entire  history  of  the  structure  is  known. 

2.I.I4.  External  Influences  on  a  Structure.  The  influences 
acting  on  a  structure  from  the  outside  fall  into  two  categories 
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the  forces  and  displacements  imposed  on  its  boundary  surfaces 
and  the  body  forces  acting  as  a  result  of  the  attraction 
between  the  mass  of  the  structure  and  the  mass  of  the  earth. 

At  every  location  on  the  face  of  a  structure, 
the  imposed  surface  force  (i.e.  force  per  unit  area)  is 
defined  by  the  vector 

P  =  p(r,  ♦  *  V  *  t)  (2.1.6) 

In  the  literature,  a  distinction  is  sometimes  made  between 
point,  line  and  surface  forces.  This  distinction  is  arti¬ 
ficial  since  all  forces  necessarily  have  to  be  applied  over 
an  area.  It  is  not  desirable  to  introduce  into  a  struc¬ 
tural  model  study  without  justification  all  the  simplifi¬ 
cations  that  are  convenient  for  mathematical  methods. 

At  some  points  on  the  surface  ( e.g.  at  a  per¬ 
fectly  rigid  support)  displacements  rather  than  forces 
may  be  applied  to  the  structure.  They  are  specified  by 
the  displacement  vector 

u1  =  u(r,  ♦ ,  f  ,  t)  (2.1.7) 

where  superscript  1  distinguishes  the  externally  imposed 
displacements  from  those  which  are  the  result  of  the  response 
of  the  structure  to  other  external  influences  such  as 
forces.  Frequently  structures  are  subjected  to  a  combina¬ 
tion  of  boundary  forces  end  displacements. 

The  body  forces  caused  by  the  earth's  attraction 
on  the  structure  are  described  by  the  vector 

O1  =  6(r,  *,f)  (2.1.6) 

where  time  variations  in  the  earth's  field  of  attraction  are 
neglected.  The  nature  of  these  body  forces  is  quite  different 
from  that  of  the  inertia  body  forces  caused  by  the  response 
of  the  structure  to  dynamic  surface  loads  or  support  accel¬ 
erations,  as  will  be  shown  in  more  detail  in  Section  2.6.1. 
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The  structure  will  also  respond  to  temperature 
changes  at  its  surfaces.  They  are  specified  by  means  of 
the  scalar 

T1  =  T(r,  *,  ?  ,  t)  (2.1.9) 

2.1.5  Description  of  Structural  Behavior.  Once  a 
structure  with  given  geometric  and  material  properties  and 
given  initial  conditions  is  subjected  to  a  set  of  external 
influences,  its  response  has  to  be  registered  by  measuring 
certain  physical  quantities  which  are  indicative  of  the 
nature  of  the  structure.  Three  different  aspects  of  the 
behavior  of  a  structure  are  of  particular  interest  to  the 
engineer: 

a)  Stress.  Strain  or  Displacement  Distribution  -  The 
stress  or  strain  distribution  in  a  concrete  structure, for 
instance,  is  of  importance  in  the  determination  of  the 
positioning  and  dimensioning  of  the  steel  reinforcement. 

The  displacement  distribution  is  of  importance  in  flexible 
structures  where  deformations  are  restricted  by  building 
codes. 

Since  all  the  properties  of  and  influences  on  a 
structure  are  space  and  time  dependent,  the  stress,  strain  . 
or  displacement  distributions  will  be  known  in  terms  of 
the  space  and  time  coordinates  once  the  following  functional 
relationships  have  been  determined: 


Displacements: 

u 

=  5(vt, 

8, 

c. 

X, 

h. 

Oo  , 

To, 

P, 

u1,  G1, 

T1) 

Stresses: 

a 

=  *  (Vj, 

8, 

c. 

X, 

h. 

<To, 

To, 

P, 

u1,  g1; 

T1) 

Strains: 

6 

-*  (Vlf 

8, 

c. 

X, 

h, 

ffo. 

To, 

P, 

(2 

51,  o1, 
.1.10) 

T1) 

This  is  the  most  general  form  of  solution  since  it  is  valid 
for  all  space  time  distributions  of  the  variables  in  the 
brackets.  For  any  specified  space  time  distribution  of  the 
variables,  u,  5 -  and  e  are  found  as  functions  of  r,  ♦  , f  and  t 
by  direct  substitution. 
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Frequently,  the  analysis  is  simplified 
considerably  if  the  distribution  of  5,5,  and  6  has  to 
be  found  for  just  one  specific  space  time  distribution 
of  the  variables.  These  investigations  are  referred  to  as 
case  studies,  and  the  results  are  obtained  directly  as 

u  =  u(r,  *  ,  ,  t) 

*  =  o  (r,  e  ,  <p ,  t) 

t  =  6  (r,  ? ,  t)  (2.1.11) 

u.  o  and  &  are  actually  related  through  the  equation  of 
state  of  the  structural  material  and  the  definition  of 
strain: 


dy  (2.1.12) 

where  y  is  the  distance  measured  in  the  direction  of  u  and 
C.  Since  the  stress-strain  relationship  is  not  always  unique 
and  in  addition  changes  with  the  stress  history  of  the 
material  as  will  be  shown  in  Section  2.3,  and  since  dis¬ 
placements  have  to  be  obtained  from  strains  by  a  frequently 
difficult  integration,  it  is  customary  in  model  studies  to 
observe  them  as  if  they  were  independent.  The  possibility 
of  obtaining  u  and  e  separately  in  a  physical  model  is  one 
of  the  main  advantages  of  that  method  over  current  mathe¬ 
matical  approaches. 

b)  The  Ultimate  Load  -  A  knowledge  of  the  stress,  strain 
or  displacement  distribution  at  a  load  level  below  that 
required  for  failure  in  general  does  not  constitute  sufficient 
information  to  predict  the  failure  load  itself.  The  deter¬ 
mination  of  the  ultimate  load,  which  is  of  paramount  impor¬ 
tance  to  the  safety  of  the  structure,  is  thus  a  problem  of 
its  own.  It  requires  the  specification  of  an  exact  failure 
criterion,  which  in  most  cases  depends  on  the  purpose  of 
the  structure  studied.  It  may  be  decided  to  impose  limits  on 
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case  studies,  and  the  results  are  obtained  directly  as 

u  =  u(r,  £  ,  $  ,  t) 

9  =  a  (r,  e  ,  <p ,  t) 

fe  =  6  (r,  *,  t)  (2.1.11) 

u,  9  end  £  are  actually  related  through  the  equation  of 
state  of  the  structural  material  and  the  definition  of 
strain: 

_  du 

dy  (2.1.12) 

where  y  is  the  distance  measured  in  the  direction  of  u  and 
£.  Since  the  stress-strain  relationship  is  not  always  unique 
and  in  addition  changes  with  the  stress  history  of  the 
material  as  will  be  shown  in  Section  2.3,  and  since  dis¬ 
placements  have  to  be  obtained  from  strains  by  a  frequently 
difficult  integration,  it  is  customary  in  model  studies  to 
observe  them  as  if  they  were  independent.  The  possibility 
of  obtaining  u  and  C  separately  in  a  physical  model  is  one 
of  the  main  advantages  of  that  method  over  current  mathe¬ 
matical  approaches. 

b)  The  Ultimate  Load  -  A  knowledge  of  the  stress,  strain 
or  displacement  distribution  at  a  load  level  below  that 
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1)  the  absolute  displacement  of  points  In  the 
structure  to  prevent  excessive  deflections  associated 
with  excessive  strains. 

ii)  the  relative  displacement  of  points  in  the 
structure  to  prevent  rupture. 

iii)  spacing  and  width  of  cracks  in  reinforced 
concrete  structures. 

It  will  be  shown  in  section  2.5  that  the  choice  of  the 
failure  criterion  has  considerable  influence  on  both  the 
outcome  and  the  reliability  of  a  model  investigation. 

c)  The  Buckling  Load  -  Thin  structures  such  as  shells 
and  thin  structural  elements  such  as  flanges  of  I-beams 
sometimes  undergo  large  deflections  without  excessive 
strains  if  they  are  subjected  to  compressive  loads.  They 
are  then  said  to  have  buckled.  The  determination  of  a 
satisfactory  general  buckling  criterion  for  the  mathematical 
model  is  at  present  still  one  of  the  major  problems  of  struc¬ 
tural  engineering.  The  physical  model  here  shows  a  very 
marked  superiority.  It  is  possible  in  most  cases  to  recog¬ 
nize  a  buckling  failure  by  direct  observation.  Most  shells, 
for  instance,  show  a  sudden  snap-through.  If  the  material 
is  still  in  its  elastic  range  when  buckling  occurs,  the 
original  shape  of  the  structure  is  regained  once  the  load 
is  removed. 

2.2  THE  RELATIONSHIP  BETWEEN  MODEL  AND  PROTOTYPE 

It  is  essential  to  keep  in  mind  at  all  times  that 
the  information  eventually  to  be  obtained  from  a  structural 
model  study  concerns  the  behavior  of  the  prototype,  not  that 
of  the  model.  Since  all  the  measurements  are,  however,  carried 
out  on  the  model,  it  is  necessary  to  develop  laws  according 
to  which  the  model  results  can  be  extrapolated  to  the 
prototype.  For  the  derivation  of  these  laws,  assume  that 
both  model  and  prototype  are  fully  described  by  the  physical 
quantities  defined  in  Section  2.1.  Model  and  prototype  will 
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then  be  related  If  the  relationship  of  each  of  these  quan¬ 
tities  individually  is  described  by  means  of  a  "correlation 
function"  f.  In  the  literature  (3*31)  ^  these  corre¬ 

lation  functions  are  frequently  defined  as  constant  ratios. 

In  a  more  general  theory,  they  may  be  made  functions  of  the 
basic  coordinates  T,  ^  and  t.  It  will  then  be  possible, 
for  instance,  to  use  distorted  models  without  violating  any 
of  the  similitude  requirements  derived  later  in  this  section. 
The  correlation  function  concept  will  also  be  found  useful 
in  Section  2.5  when  structural  model  analysis  i3  considered 
as  a  random  process. 

2.2.1  Geometric  Properties.  The  distance  r  and  time 
t  are  the  only  physical  quantities  required  for  the  descrip¬ 
tion  of  the  geometric  properties.  Thus  model  and  prototype 
are  related  if 

rm  "  Mr-  «•  ?•  rp 

**  =  ft(r-  t-  lp  (2.2.1) 

where  the  following  conventions  are  observed: 

1)  Subscript  m  denotes  model  property 
ii)  Subscript  p  denotes  prototype  property 
Hi)  f_  denotes  correlation  function  for  the 
physical  quantity  z  and  is  evaluated  at  r,«-,  9  and  t  of 
the  prototype. 

Since  angles  are  dimensionless  ratios,  it  follows  that 

♦  =  * 

m  p 


and 


=  9 


p 


2.2.2  Material  Properties.  Model  and  prototype  are 
related  as  follows: 


«m  =  t)gp 

cm  =  fc(r'*'?'  t)cp 


Specific  Mass: 
Specific  Heat: 
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Coefficient  of  Thermal  x  =  f_(r,«>,f,t)x 

Conductivity!  p 

Coefficient  of  Thermal  h^  =  fw(^t*  ♦  »  9*  t)h_ 

Expansion:  m  * 

(2.2.2) 


2.2.3  Initial  Conditions.  Model  and  prototype  are  related 
if 


Initial  Stress 

*°m  = 

^(r,  ♦  »  9  ,  t0)?op 

Initial  Temperature 

T*m  = 

fT(r,  *• ,  9 ,  to )Top 

2.2.4  External  Influences 

on  a  Structure. 

(2.2.3) 

Model  and 

prototype  are  related  if 

Surface  Forces: 

= 

fp(N,  ♦,  9,  t)Pp 

Displacements : 

ii 

If 

fu(r»  e’»  9* 

Body  Forces: 

• 

<f  ,  t)  Gp 

Temperature: 

n 

fT(r,e  ,  f,  t)T* 

(2.2.4) 


It  should  be  noted  that  the  correlation  functions  defined 
in  Sections  2.2.1  to  2.2.4  are  independent  of  the  state  of 
stress  or  strain  of  the  materials  involved. 

2.2.5  Comparison  of  Structural  Behavior,  In  addition 
to  the  quantities  described  in  Section  2.1.5*  the  acceleration 


a  = 


a2u 


The  required  correlation  functions  thus  are 
Displacements: 

Stresses : 


(2.2.5) 

also  be 

considered. 

are 

= 

Vp'f  »  *  » 

t)Up 

O’  = 
m 

*  ?» 

t)rp 
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Coefficient  of  Thermal 
Conductivity: 

*m 

*” 

fx(r'*»  *  »t)jcp 

Coefficient  of  Thermal 
Expans icn: 

S 

V*'*'  t)hp 

(2.2.2) 

2.2,3  Initial  Conditions. 

Model  and  prototype  are  related 

if 

Initial  Stress 

s 

f0-(r»  ♦»  f  »  t0)?o 

Initial  Temperature 

T‘m 

= 

fT(r,  *  ,  ?  ,  to) To 

(2.2.3) 

2.2,4  External  Influences 

on  a  Structure 

.  Model  and 

prototype  are  related  if 

Surface  Forces: 

S 

fp(N,  *,  t)Pp 

Displacements : 

“i 

= 

fu(r»  <t>  t)uj 

Body  Forces: 

s 

fff(  r»  &  >  ,  t)  6p 

Temperature: 

s 

i>Qi(  **  *  9  >  <f  t  t)  Tp 

(2.2.4) 


It  should  be  noted  that  the  correlation  functions  defined 
in  Sections  2.2.1  to  2,2.4  ar®  independent  of  the  state  of 
stress  or  strain  of  the  materials  involved. 

2.2,5  Comparison  of  Structural  Behavior.  In  addition 
to  the  quantities  described  in  Section  2.1.5»  the  acceleration 


a  = 


si3 


Displacements: 


(2.2.5) 

also  be 

considered. 

are 

ii 

t)Up 

O’  a 

m 

Vr»#  *  f  *  t)frp 

Stresses : 


Accelerations: 


fa(r,  «■#  <?,  t)ap 

(2.2.6) 

(2.2.7) 

so  that  no  correlation  function  for  strains  needs  to  be 
specified. 

2.2.6  Dimensional  Analysis.  If  all  the  correlation 
functions  defined  in  the  previous  sections  were  independent, 
the  prototype  behavior  could  not  be  predicted  from  the  model 
behavior.  An  analysis  of  the  dimensions  of  the  physical 
quantities  defining  a  structure  will,  however,  reveal  that 
only  a  few  of  the  correlation  functions  can  be  considered 
as  independent.  Thus  the  correlation  functions  in  equa¬ 
tion  (2.2.6)  comparing  model  and  prototype  behavior  can  be 
expressed  in  terms  of  the  other  correlation  functions,  and  a 
prediction  of  the  prototype  behavior  is  possible. 

Before  the  dimensional  analysis  is  carried  out, 
a  special  type  of  model  study  that  does  not  require  this 
approach  will  be  discussed  for  the  sake  of  completeness 
of  presentation.  Suppose  that  quantity  y  describing  the 
structural  behavior  is  known  as  a  function  of  other  quanti¬ 
ties  x^,  •  ••»  i.e.,  the  physical  event  controlling  y  is 

known.  Say 

X.X-  X. 

y  =  -1  —  +  Ji  -  ...  (2.2.8) 

x3  x5 

Chose  x^,  Xg  ...  in  the  model  and  prototype  in  such  a  manner 
that 

*lm*2m  =  k  *lp*2p 
*3m  X3P 


m 


It  follows  from  equation  (2.1.12)  that 


m 


'm 


f  y_ 
r  “  p 
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k 


3a 

X5m 


x5p 


etc. 


Then  equation  (2.2.8)  can  be  used  to  predict  the  proto¬ 
type  behavior  yp  from  the  model  behavior  yn  since 


^m 


(2.2.9) 


The  usefulness  of  this  method,  which  has  been  very  popular 
In  the  past.  Is  restricted  since  equation  (2,2.8)  must  first 
be  found  by  mathematical  investigation.  The  physical  model 
is  degraded  to  become  an  "analog  computer"  for  the  mathe¬ 
matical  model  and  in  this  manner  is  deprived  of  its  main 
advantage  of  generality. 

Full  generality  is  maintained  if  the  model  to 
prototype  relationship  is  based  on  an  analysis  of  the 
dimensions  of  the  physical  quantities  describing  the  struc¬ 
ture.  The  following  concepts  are  of  importance: 

Dimension:  The  dimension  is  that  part  of  the 
description  of  a  variable  which  specifies  the  type  of  physi¬ 
cal  measurement  that  has  to  be  carried  out  to  determine  its 
magnitude.  The  dimension  is  dictated  by  nature  and  is  there¬ 
fore  independent  of  the  basic  unit  adopted  for  measurement, 
which  may  be  arbitrarily  chosen  by  man.  If  the  respective 
magnitudes  of  two  physical  quantities  with  the  same  dimen¬ 
sion  are  determined  by  comparison  to  two  different  basic 
units,  they  cannot  be  compared  directly.  Yet  the  quantities 
are  not  totally  unrelated  since  their  magnitudes  can  be  com¬ 
pared  once  the  ratio  between  the  basic  units  is  specified. 
This  principle  forms  the  basis  for  dimensional  analysis  of 
structural  models. 

Dimensional  Relationship:  Some  basic  physical 
laws  such  as  Newton's  second  law  relating  mass  M,  length  L 
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•  and  time  t  to  force  P  make  it  possible  to  regard  certain 

quantities  as  physically  completely  equivalent  to  a  pro¬ 
duct  of  other  terms,  e.g., 

P  =  2%  (2.2.10) 

t2 

The  magnitude  of  each  of  the  four  quantities  in  this  equation 
is  determined  by  expressing  it  as  a  multiple  of  a  chosen 
basic  unit,  such  as  a  foot,  a  pound,  a  meter,  etc.  The 
numerical  values  that  are  to  be  substituted  into  the  equation 
thus  change  with  the  basic  units  of  measurement  that  are 
chosen.  The  type  of  measurement  to  be  performed,  e.g., 
determining  a  length,  determining  a  mass,  is,  however, 
independent  of  the  basic  units  that  are  chosen.  The  type 
of  measurement  required  to  determine  the  magnitude  of  a 
given  variable  is  called  its  dimension.  If  we  let  P  represent 
the  magnitude  only  of  the  variable  (e.g.  10  feet,  20  pounds) 
and  D  (P)  the  dimension  only  (e.g.,  length,  mass)  then,  in 
addition  to  equation  (2.2.10),  we  can  write 

D(F)  =  DfiJL'W  ■  -  P-ffi  (2.2.11) 

\t  D(t)2 

The  left  and  right  hand  sides  of  this  equation  are  independent 
of  the  basic  units  chosen  for  the  variables,  so  that  the 
equation  is  more  characteristic  of  the  phenomenon  than 
equation  (2.2.10)  was.  Equation  2.2.10  could  also  be  written 

Ft2  =  ML 

so  that  equation  (2,2.11)  would  become 

D(P)  [D(t)]  2  =  D(M)  D(L) 

Thus  multiplication  and  division  of  dimensions  are  both  defined. 

Fundamental  Physical  Quantities:  There  exist  a 
number  of  quantities  such  as  length,  mass  and  time  whose 
magnitude  can  be  determined  by  one  specific  type  of  physical 
measurement  only.  They  will  be  called  the  fundamental  phys¬ 
ical  quantities.  Their  number  depends  on  the  type  of  exper¬ 
iment  performed  since  heat,  for  Instance,  may  be  considered 
a  fundamental  quantity  only  if  thermodynamic  transformations 
are  not  studied  in  the  experiment. 
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If  a  physical  quantity  is  not  fundamental,  its 
magnitude  can  be  determined  by  several  types  of  measure¬ 
ment.  Force,  for  instance,  can  be  determined  from  measure¬ 
ments  of  mass,  length  and  time  or  mass  and  acceleration 
or  directly  by  comparison  to  a  unit  force. 

Dimensional  Independence:  It  follows  from  the 
definition  of  fundamental  physical  quantities  that  the 
magnitude  of  none  of  them  can  be  determined  by  a  combina¬ 
tion  of  the  types  of  measurement  used  to  determine  the  other 
fundamental  quantities.  Thus  the  fundamental  quantities  are 
considered  dimensionally  independent.  The  set  of  all  the 
fundamental  quantities  for  a  specific  experiment  will  thus 
contain  all  the  types  of  measurement  to  be  carried  out 
during  the  experiment,  and  is  thus  a  basis  for  the  dimen¬ 
sions  of  the  other  quantities  describing  the  phenomenon. 

Consider  the  dimensions  D(q^)  of  all  the  quanti¬ 
ties  that  play  a  role  in  a  specific  experiment.  Let  i 
range  from  1  to  n.  If  Kj  with  j  ■  1,  2,  . . . ,  r  is  the 
complete  set  of  fundamental  physical  quantities  for  the 
experiment  (  ••.  n  jg  r),  each  of  the  dimensions  D(<^)  can  be 
expressed  as  a  function  of  all  the  D(Kj),  in  the  manner 
demonstrated  by  equation  (2.2.11).  Thus 

D^)  =  D1(K1,  K2,  ...,  Kp)  (2.2.12) 

where  i  =  1,  2,  ...,  n 

Since  the  processes  of  multiplication  and  division  of 
dimensions  are  defined,  the  laws  for  the  inversion  of 
mathematical  functions  can  be  applied  to  equations  (2.2.12). 
Thus  if  any  r  of  the  n  equations  (2.2.12)  are  selected  in 
such  a  manner  that  the  Jacobian 

'b  ^  2*  ^  ^  0 

d(kr  k2,  ...»  Kr)  (2.2.13) 
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it  ia  possible  to  write  explicitly 

D( K j )  —  Dj  ^  2’  9  r^  (2.2.14) 

where  J  =  1,  2,  ...»  r 

Since  the  Kj  form  a  complete  set  of  dimensionally  inde¬ 
pendent  quantities,  equation  (2.2.1k)  guarantees  that 

q  ...»  also  form  a  complete  set  of  dimensionally 
independent  quantities  for  the  experiment.  Thus  functional 
Independence  is  not  restricted  to  the  fundamental  quantities. 
Substitution  of  equation  (2.2.14)  into  the  remaining  (n-r) 
equations  of  expression  (2.2.12)  shows  that 

®(  S  |j)  “  D^(q  (2.2.15) 

where  h  =  r,  ( r  +  1 ),...,  n. 


Dimensional  Analysis  -  Equation  (2.2.15)  indicates 
which  quantities  have  to  be  measured  before  q  h  is  known. 

As  discussed  in  the  section  on  dimensional  relationships, 
equation  (2.2.15)  must  be  augmented  by  an  equation  relating 
the  magnitudes  of  the  quantities  concerned: 

^  h  —  Ph^  i»  ^  2'  •  •  *  *  9  (2.2. 16 ) 


D^  in  equation  (2.2.15)  and  in  equation  (2.2.16)  are 
functions  of  identical  form,  i.e.,  if 


Dh^l'  ^  2'  ^ 


Ki>] a  b 

(p<,3>r 


then  it  follows  that 

(qJ8  (o2)b 

ph^l»^2»^3)  =  - — 7c~ 

'13} 

The  number  of  variables  in  the  event  has  in  effect  been 
reduced  from  n  to  r.  While  the  basic  units  of  measure¬ 
ment  for  q  q  2*  •  can  still  be  chosen  arbitrarily, 

the  basic  unit  of  measurement  for  must  now  be  calculated 
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from  equation  (2.2.16).  This  is  the  basic  principle  of 
dimensional  analysis. 

2.2.7  The  Laws  of  Similitude.  Suppose  that  the  same 
physical  events  occur  in  both  the  model  and  the  prototype. 
Then  the  form  of  function  D^,  and  therefore  also  the  form 
of  function  F^,  is  the  same  for  both  structures.  Let  sub¬ 
scripts  m  and  p  denote  model  and  prototype  properties  re¬ 
spectively.  Then 

^  hm  _  Fh(^lm>  ^  2m*  J 

%p  Ph(<llp'  <l2p*  •**'<*rp) 

Since  the  relationship  between  model  property  and  proto¬ 
type  property  qp  has  already  been  defined  as 

-  *1  m  “  p 

the  preceding  equation  can  be  rewritten 

fqh^hp  B  V^l^lp*  f«2^2p»  fqr<^rp) 

^hp  Fh^lp»  ^2p*  *'*'<lrp) 

Since  is  always  a  product  or  quotient  of  the  variables 
in  the  brackets,  it  will  always  be  possible  to  factor 

Fh(  W  fH2*  •••»  fq Ph^lp*l2p* 

Ph'llp'  ^ 2p’  •••»  Srp) 

Fh ( ^*<j i *  fqg*  •••*  r )  (2,2,17) 

where 

h  r,  (r  +  1),  •  •  • ,  n. 

The  original  goal  of  the  dimensional  analysis,  to  establish 
the  relationship  between  the  correlation  functions,  has 
been  achieved  in  equation  (2.2.17).  The  correlation  func¬ 
tions  may  simply  be  regarded  as  the  relationship  between  the 


so  that 


qh 


Lqh 
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basic  units  of  measurement  in  model  and  prototype.  It  is 
thus  possible  to  conclude  thst 

i)  there  are  as  many  independent  correlation  func¬ 
tions  as  there  are  dimensionally  independent  physical 
quantities. 

ii)  the  relationship  between  the  correlation  func¬ 
tions  is  of  a  form  identical  to  that  of  the  relationship 
between  the  dimensions  of  the  corresponding  physical  quanti¬ 
ties. 

The  laws  of  similitude  are  defined  as  the  relations  which 
must  exist  between  corresponding  physical  quantities  if  a 
one  to  one  correspondence  between  model  and  prototype  is 
to  exist.  They  are  easily  derived  from  equation  (2.2.17) : 

^ hm  =  ft|h^hp  =  ^hpPh^f9l'  f92>  ••••  f9r^ 

(2.2.18) 

where  h  =  r,  ( r  +  1)  ,  . . . ,  n. 

If  the  (n-r)  laws  of  similitude  (2.2.18)  are  obeyed ■,  model 
measurements  can  be  used  to  predict  prototype  behavior, 

2.2.8  Laws  of  Similitude  for  Quantities  in  Section  2.1. 
Both  the  number  of  variables  and  the  number  of  fundamental 
physical  quantities  vary  from  experiment  to  experiment  and 
have  to  be  guessed  intuitively  by  the  experimenter.  If 
superfluous  variables  are  considered,  unnecessary  restric¬ 
tions  are  imposed  on  -the  model  since  the  number  of  laws  of 
similitude  (2.2.18)  will  be  larger  than  actually  required. 

If  an  essential  variable  is  neglected,  the  model  will  not 
satisfy  one  of  the  essential  laws  of  similitude  and  the 
prediction  of  the  prototype  behavior  will  be  erroneous. 

The  reasoning  outlined  in  Sections  2.2.7  and 
2.2.8  will  now  be  illustrated  by  the  study  of  a  specific 
case:  a  structure  fully  described  by  the  quantities  defined 
in  Section  2.1.  To  conform  with  practice,  it  is  assumed 
that  the  model  will  not  be  used  to  investigate  thermodynamic 
transformations.  Then  all  the  types  of  physical  measurement 
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employed  In  the  experiment  are  contained  among  the  dimensions 
of  the  following  five  fundamental  physical  quantities: 


Table  2,2.1  -  The  fundamental  physical  quantities  for  the 
properties  described  in  Section  2.1. 


Quantity 

Symbol 

Ft,  Lb,  Sec  Unit 

Mkgs  Unit 

Length 

L 

Ft 

Meter 

Mass 

M 

Lb  Sec2/Ft 

Kg  Sec2/M 

Time 

t 

Sec 

Sec 

Temperature 

T 

°F 

°C 

Heat 

H 

B.T.U. 

Calorie 

Equations  (2.2.12)  are  summarized  in  Table  2.2.2. 

Table  2,2,2  -  Powers  to  which  quantities  in  equation  (2.2.12) 
have  to  be  raised. 


r 

t 

c 

X 

h 

O' 

T 

G 

a 

L 

1 

0 

-3 

0 

-1 

0 

-1 

0 

-2 

1 

M 

0 

0 

1 

-1 

0 

0 

1 

0 

1 

0 

t 

0 

1 

0 

0 

-1 

0 

-2 

0 

-2 

-2 

T 

0 

0 

0 

-1 

-1 

-1 

0 

1 

0 

0 

H 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

It  will  be  noted  that 

1)  The  equation  of  state  does  not  contribute  any 
additional  variables  to  Table  2,2.2,  as  will  be  shown  in 
Section  2.3. 

ii)  P  and  u  are  not  listed  in  Table  2.2.2  since  they 
are  dimensionally  identical  with  O’  and  r  respectively.  This 
guarantees  that  f  =  f.and  f  =  f  .  (where  p  is  the  surface  stress 
The  independent  quantities  must  now  be  chosen 

subject  to  condition  2.2.13.  Suitable  choices  are  discussed 
in  Section  2.3}  in  this  example  r,  c,  h,  0  and  a  are  used: 
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D(r) 

=  D(L) 

D(0) 

=  D( 

D(c) 

=  D(S_) 

DU) 

=  D( 

MT 

D(h) 

=  D(i) 

T 


M  x 


(2.2. 
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Condition  2.2.13  becomes 


d(rf  c,  h,  G,  a) 
fc(L,  M,  t,  T,  H) 


0 

0 


2M 


7 


M  L  t5 


H 

~7~ 
M  T 

0 


~2  2 


t  0 


0 

0 

2M 

T"T 

ITt-5 

2L 

‘  ? 


-  -5p 
mt£ 
i 


Thus  equations  (2.2.19)  can  be  inverted  to  yield 


D(L)  =  D(r) 

D(M)  = 

xad) 

a 

D(T)  =  D(i) 
h 

D(H)  = 

D(Sfitl) 

ah 

D(t)  =  D(-/f  ) 

(2.2.20) 

The  dimensionally  dependent  variables  can  now  be  expressed 
in  terms  of  the  dimensionally  independent  variables  as  in 
equation  (2.2.15) - 

D(t)  =  ) 

D(<r)  = 

D(  Or) 
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D(g) 

=  D(-) 

D(T)  =  D(i) 

a 

h 

D(x) 

=  D(2£j!> 

a* 

(2.2.21) 

The  form  of  functions  is  thus  known.  The  form  of  func¬ 
tions  F,  is  therefore  also  known  so  that  the  laws  of  simili- 
h 

tude  can  be  written  down  directly: 


t 


m 


m 


=  Wp 


«m  = 

'T 

a 

a 

T  = 

f  f  f  a 

c  g  r 

m 

C 

m 


(2.2.22) 


Equations  (2.2.22)  are  the  equivalent  of  equation  (2.2.18) 
in  the  general  theory. 


2.3  THE  EQUATIONS  OF  STATE  OF  THE  MODEL  AND  PROTOTYPE  MATERIALS 


The  equation  of  state  relates  those  of  the  quantities 
describing  the  material  which  are  dependent  on  the  external 
influences  acting  on  the  structure.  This  relationship  may 
be  time  and  temperature  dependent,  it  may,  or  may  not,  vary 
from  point  to  point  in  the  material  and  it  may,  or  may  not, 
be  different  for  different  directions  at  the  same  point.  In 
addition,  the  events  occurring  in  one  direction  msy,  or 
may  not,  affect  the  events  occurring  in  all  other  directions 
(Poisson  effect).  There  is  some  question  as  to  whether  the 
events  occurring  at  one  point  affect  the  events  occurring 
at  neighboring  points,  eg.  whether  the  stress  gradient 
has  an  effect  on  the  crack  formation  in  concrete^^.  The 
general  form  of  the  equation  of  state  will  be  developed  in 
stages. 
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FIGURE  2.3.1 

STRESS  —  STRAIN  DIAGRAM  OP  A  FICTITIOUS  MATERIAL 


CONCRETE  0X466  10  YEARS. 

TIME  AFTER  LOADING 

FIGURE  2.3.3 

VARIATION  OP  CREEP  STRAIN  WITH  TIME  FOR  CONCRETE,  STEEL  AND 
PLEXIGLAS.  NOTE  THAT  TME  CURVES  ARE  AFFECTED  CONSIDERABLY  EY 
MATERIAL  COMPOSITION,  EVlRONMfNT  ANP  OTHER  FACTORS. 
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STRftSt  STMLEE 


COMPATIBLE 


INCOMPATIBLE 


FIGURE  2.3.2. 

EXAMPLES  OF  COMPATIBLE  AND  INCOMPATIBLE  MODEL  AND  PROTOTYPE  MATERIALS 
-  MODEL  MATERIAL  -  PROTOTYPE  MATERIAL 
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2.3.1  The  Stress-strain  Relationship.  Assume  that  a 
material  is  very  gradually  subjected  to  an  external  influence 
for  the  first  time  in  its  history  in  such  a  manner  that  all 
material  properties  except  the  state  of  stress  and  strain 
are  kept  constant.  The  event  may  be  described  by  the  graph 
in  Figure  2.3.1.  This  curve,  which  is  particular  for  the 
material  studied,  can  be  described  mathematically  as 

<r  =  ae  +  be2  +  ...  \2.3.1) 


The  physical  nature  of  the  event  is  not  evident  from  this 
equation,  since  it  is  not  dimensionally  homogeneous  if  a 
and  b  are  constants.  A  better  form  is 


OO  K 

f  J— 

n! 


n*l 


(2.3.2) 


gr-o 


where  n  is  a  constant  and  the  equation  is  dimensionally 
homogeneous.  If,  in  a  more  general  case,  the  derivatives 
are  known  at  point  A  =  (  *.,  ),  this  equation  becomes 


(2.3.3) 


where  (-OmflJt  *  6  *  (+«)maJc. 

Both  equation  (2.3.3)  and  its  limits  of  validity  (i.e., 

-€max^  may  vary  from  point  to  point  in  the  structure. 

The  superiority  of  equation  (2.3.3)  over  equation  (2.3.1) 
will  become  evident  when  the  equations  of  state  of  model  and 
prototype  materials  are  compared  later. 

2.3.2  The  Effect  of  Time,  The  time  variable  affects 
the  equation  of  state  in  two  different  manners: 

a)  Aging  may  change  the  shape  of  the  curve  in  Fig.  2.3.1 


(2.3.10 
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and 


b)  The  rate  at  which  the  external  influences  are  applied 
may  affect  the  shape  of  the  curve  in  Pig.  2.3.1.  The  rate 
of  application  is  measured  by  means  of  the  strain  rate  ^ 


2.3.3  The  Effect  of  Temperature.  The  effect  of  tem¬ 
perature  on  the  equation  of  state  is  twofold  like  that  of 
time: 

a)  It  may  change  the  shape  of  the  curve  in  Pig.  2.3.1: 


*cr 


and 


e-i 


(t-t.)£ 

0  or 

t! 

St"  ST1 

(T-T.)1 

a1*™. 

l\ 

T»T. 


T.T. 


(2.3.6) 


b)  A  change  in  temperature  may  cause  the  material  to 
expand  or  shrink,  setting  up  internal  stresses  and  strains. 

This  fact  has  already  been  taken  into  account  by  the  defini¬ 
tion  of  the  coefficient  of  linear  expansion. 

2.3.14-  Materials  Subjected  to  Several  Load  Cycles. 

If  the  material  in  Pig.  2.3.1  is  stressed  to  point  A  =  (60*^0) 
and  the  load  is  then  removed,  it  may  follow  curve  AB  rather 
than  curve  AO.  Thus  the  stress-strain  curve  at  A  cannot  be 
described  uniquely,  but  is  dependent  on  a  logical  decision: 
Follow  OA  on  first  loading  cycle,  continue  along  AC  if  load 
is  Increased  further,  follow  AB  if  load  is  decreased.  AB 
can  be  described  in  a  manner  similar  to  the  description  of  OA: 
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O’*®; 


where  the  bars  in 
is  described. 


indicate  that  the  "curve  when  unloading" 


The  permanent  deformation  associated  with  this 
phenomenon  may  lead  to  a  considerable  reduction  of  the  initial 
stresses  caused  by  the  manufacture  of  the  structure.  There¬ 
fore,  if  a  prototype  structure  is  put  through  several  load 
cycles  before  it  is  subjected  to  the  maximum  load  of  its 
lifetime,  and  the  model  has  been  manufactured  in  a  manner 
analogous  to  the  prototype  construction,  the  loading  history 
should  be  faithfully  copied  in  the  model. 

2.3.5  The  General  Equation  of  State.  The  results 
obtained  so  far  can  be  combined  if  it  is  assumed  that  the 
effects  of  time  and  temperature  are  independent.  The 
general  equation  of  state  then  is 


0*(r,  *,?,  t,  T,^l)  =  O'  o(r,-e-,t,  tc,  T0,  ^ 

fct  Bt 

oo  ,  n"  f  oo  L  -."•j 

’T  («-«.)  J  Bo-  /,  t  Be  \  VU-t.V,,  ^  <y 


)  + 


f.  fat)1  jcr 

fa  t! 


3T4  +  k!  L  >e"-l 

t*t,  'ot'  ; 


(2.3.7) 

where  the  limits  of  the  maximum  positive  strain,  for  instance, 
are  given  by: 


"’It1  = 

t  ^ 


max' 


)  + 


f  (t-g*  ay—  f  ct-o  bc— 

—  Bti  fc  t!  BT1 


't»t. 


t.t. 


(2.3.8) 
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Similar  equations  can  be  written  for  the  "curve  when  unloading". 
They  are  far  too  complicated  for  use  in  mathematical  investi¬ 
gations,  but  can  be  used  to  advantage  in  the  detennination 
of  the  similitude  requirements  for  the  equations  of  state 
of  the  model  and  prototype  materials.  It  should,  however, 
be  noted  that  ma'terial  properties  auch  as  creep,  relaxation, 
the  Poisson  effect,  etc.  which  are  discussed  in  Sections 
2.3.7  to  2.3.10  are  not  reflected  in  these  equations. 

It  has  been  shown  in  Section  2.2  that 


«»  *  % 


m 


f  o* 

*  p 


t  =  ft 
m  t  p 


m 


Vp 


(2.3.9) 


Term  by  term  comparison  of  equation  (2.3.7),  written  once 
for  the  model  and  once  for  the  prototype,  yields  the  follow¬ 
ing  similitude  requirements  if  equation  (2.3.9)  is  used: 


°ia. 


f 


r  % 


Ap  .  r  A* 


V  (tf&A  V  I  Wf\ 

St l\s&ry[9 

^p 


2 

* 


(2.3.10) 

Similarly,  the  similitude  requirements  for  equation  (2.3.$) 
are 
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^  (»— .1-  .  4-  . 

*2;  if  8t£ 

8*  feiwftAw  .  P  **  ^  C^mw)f 

ac&T  »(*;? 


Once  the  requirements  of  equations  (2.3.10)  and  (2.3.H) 
are  satisfied  for  all  points  (V,  and  all  values  of 

n,  j,  k  and  the  model  and  prototype  materials  are 
compatible. 

2.3.6  Examples  of  Equations  of  State. 
a)  Suppose  equations  (2.3.7)  and  (2.3.8)  for  a  particu¬ 
lar  prototype  material  have  the  form 

*V*[«-*btf  +  cTi»‘  +  cl  Hj]e» +  [4  +%Tf] 

(*««).  *  kit.  +  e.T»  ♦  di  J 

It  is  desired  to  determine  the  equation  of  state  for  the 
model  material  if  the  values  of  ft  and  fT  are  prescribed. 
Comparison  to  equations  (2.3.7)  and  (2.3.8)  shows  that  for 
the  prototype 


te  III)*  * 

1r 

**  -b 

]>*">** ,  w 

*t  ' 

&6^T*-*Zc 

lllMM  .  C. 
VT  C‘ 

HI.) '2* 

^£hmu«  a  A, 

3<St>  *■ 

According  to  equations  (2.3.10) 
lng  values  for  the  model  are 


fee* 


T 

M 


and  (2.3.11)  the  correspond- 


*  *c 

mr'Tr\ 
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*  .il  b 

c 

(If)  *  f»{tA 

Thus  the  equation  of  state  for  the  model  material  is 


b)  Assume  that  due  to  technological  considerations, 
both  model  and  prototype  material  are  specified.  It  is 
desired  to  determine  under  which  conditions  the  two  materials 
will  be  compatible.  Suppose 

%  *  «-«p  +  [  b*  c  C  |f '  ~  <•> 

«•»-  t  A+  *  (T»-T„)]  e„  +  (ii) 

According  to  equation  (2.3.10) 

f .  =  -  f.  =  =  SL. 

9  d  r  cfg.  ac 

Equation  (i)  does  not  contain  a  term  corresponding  to 

e(Tm  -  Tom) €  m  In  equation  (ii).  Thus  the  model  material 

will  satisfy  the  similitude  requirements  only  if  the  model 

is  kept  at  a  constant  temperature  T_  =  T0«  so  that  this 

n  n 

term  vanishes.  Equation  (li)  does  not  contain  a  term 

p 

corresponding  to  b«p  in  equation  (i).  Since  «p  ^  0,  the 
two  materials  are  not  compatible. 


at 

3T 

'd  6  mu 


■t 

'% 
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c)  Fig.  2.3.2  shows  several  examples  of  compatible  and 
incompatible  model  and  prototype  materials.  Proofs  of 
compatibility  and  incompatibility  proceed  as  in  Sections 
(a)  and  (b)  above. 

2.3.7  Creep  and  Relaxation.  Creep  and  relaxation  are 
time  dependent  phenomena  which  have  not  been  taken  into  account 
in  the  general  equation  of  state.  Creep  describes  those 
phenomena  which  occur  in  the  material  if  an  external  force 
of  constant  magnitude  remains  applied  for  some  period  of 
time,  and  to  some  extent  after  the  influence  has  been 
removed  (creep  recovery).  According  to  F.  Leonhara^ 
the  total  creep  strain  in  concrete  is 

6k  =  $  x  ki  x  k2  *  6el  (2.3.12) 

where 

=  elastic  strain 

J  =  variable  dependent  on  relative  humidity. 

k,  =  variable  dependent  on  the  duration  of  load 
application. 


kp  =  variable  depending  on  the  composition  of  the 
concrete  and  the  thickness  of  the  member. 

Whereas  it  appears  possible  to  control  the  variables  on  which 
§  and  k,  are  dependent  in  a  model  study,  research  remains 
to  be  done  on  factor  kg. 

Relaxation  occurs  where  a  constant  displacement 
is  applied  to  a  structure.  In  steel  prestressing  cables, 
for  instance,  it  is  dependent  on  both  the  initial  stress 
and  the  composition  of  the  material.  Thus  if  f v  /  1, 
relaxation  may  be  very  difficult  to  model.  Plexiglas,  a 
frequently  used  model  material  also  displays  creep  prop¬ 
erties.  Fig.  2.3.3.  compares  typical  creep  properties  of 
steel,  concrete  (taken  from(2^.  Fig.  2.54  and  2*2°)  and 
Plexiglas  (taken  from  Fig.  1,  p.  316). 


33 


2.3.8  The  Polaaon  Effect.  The  effect  of  the  events 
occurring  in  one  direction  at  a  point  in  a  material  on  the 
events  occurring  in  another  direction  at  the  same  point  is 
described  by  Poisson’s  Ration).  Since  this  ratio  is  dimen¬ 
sionless,  it  must  have  the  same  value  in  both  the  model  and 
the  prototype  structure.  Different  values  for  ^  in  model 
and  prototype  can  lead  to  serious  errors: 

Consider  a  case  of  plane  strain  where  strains 

«  and  e  have  been  measured  in  orthogonal  directions  at  a 

*  y 

particular  point.  Then  the  stress  in  the  direction  of 
6  x  for  a  linear  elastic  homogeneous  material  is 

- - -  f(l  -d)fi  +^6  1 

1  +*  )(1  -  3d)  L  x  yJ 

0.40  (eg.  a  plastic)  =  0.15  (eg.  concrete) 

1000  yUin/in,  €y  =  500  yULin/in 

=  2860 
=  1150  Ep 

E 

=  — ,  the  model  leads  to  a  prediction  which 

EP 

is  in  error  by  In  cases  where  the  shape  and  boundary 

conditions  of  a  structure  do  not  indicate  if  it  will  be  in 
a  state  of  plane  stress  or  plane  strain,  the  Poisson  Effect 
makes  it  necessary  to  measure  strains  at  a  point  in  three 
mutually  orthogonal  directions  before  the  stress  distri¬ 
bution  can  be  determined. 

2.3.9  Crack  Formation  and  Deflections,  It  has  already 
been  shown  that  similitude  of  deflections  in  model  and  proto¬ 
type  can  be  maintained  only  if  the  strains  in  the  two 
structures  are  the  same.  The  concrete  structures  where  the 


If 


and 


then 


m 


€*  = 


(o-  ) 

1  x'm 


«Vp 


Since  f»  =  — 


deflections  are  affected  considerably  by  the  formation  of 
cracks,  further  requirements  have  to  be  met  before  defor¬ 
mations  will  be  similar: 

a)  Strain  Level  at  Which  Crack  Formation  Starts  - 

A  crack  will  start  at  a  point  in  a  structure  if  6  =  «  max 
where  t  has  been  defined  in  Section  2.3.5.  If  relation- 

UlfiX 

ships  2.3.11  are  satisfied,  similitude  requirements  for 

the  initiation  of  cracks  will  be  satisfied.  This  is  due 

to' 

to  the  fact  that  the  stress  gradient  ,  which  in  the 

most  general  case  will  be  different  in  model  and  prototype, 
has  no  effect  on  the  crack  formation,  as  discussed  in  Ref  .  10 

b)  Distance  Between  Cracks  -  It  is  shown  in  Ref.  10 
that  for  ffi  J,  the  distance  between  cracks  in  the  model 
will  be  f  times  the  distance  in  the  prototype,  as  required 
by  similitude,  if 

i)  the  diameter  of  the  individual  reinforcing 
bars  is  scaled  correctly. 

ii)  the  percentage  reinforcement  over  the  entire 
section  is  the  same  in  model  and  prototype.  Hence,  if 
several  bars  in  the  prototype  are  replaced  by  a  single  bar 
in  the  model,  similitude  requirements  with  respect  to  crack 
formation  may  not  be  satisfied. 

c)  Width  of  Cracks,-  The  width  of  the  cracks,  which  has 
considerable  influence  on  the  deflection,  is  scaled  by  fp 

if  rr  s  i. 

2.3.10  Damping.  If  a  structure  is  subjected  to  dynamic 
loads,  its  behavior  is  considerably  affected  by  the  damping 
properties  of  its  material.  Since  the  ratio  of  the  actual 
damping  in  the  material  to  the  critical  damping  is  dimension¬ 
less,  it  must  have  the  same  value  in  model  and  prototype. 

2.4  THE  DESIGN  AND  USE  OP  STRUCTURAL  MODELS 

The  theory  presented  in  Sections  2.1,  2.2  and  2.3 
can  be  used  to  make-  the  design  of  structural  models  a  systematic 
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process.  This  section  presents  two  different  methods  of 
design,  one  for  cases  where  the  model  material  is  specified, 
the  other  for  cases  where  the  model  material  can  be  chosen 
freely.  In  both  cases  it  will  be  assumed  that  the  influences 
on  and  properties  of  the  prototype  structure  are  known. 

Where  this  is  not  the  case,  as  discussed  in  Section  2.6.2, 
additional  problems  are  encountered. 

Method  I  -  Suppose  that  technological  considerations 
make  the  use  of  a  particular  model  material  mandatory. 

Then  the  value  of  f^.  will  be  prescribed,  as  demonstrated 
in  example  (b)  of  Section  2.3.6.  The  value  of  ffc  and  fT 
may  or  may  not  be  prescribed.  The  design  then  proceeds  as 
follows : 

1)  Decide  which  fundamental  physical  quantities 
are  involved  in  the  experiment.  Let  their  number  be  r. 

ii)  Decide  which  variables  are  involved  in  the 
experiment. 

iii)  Since  there  are  r  fundamental  quantities, 
r  of  the  properties  of  the  model  may  be  chosen  arbitrarily 
as  a  dimensionally  ind  •  pendent  set  as  long  as  they  satisfy 
equation  (2.2.13).  If  fr,  f^  and  fT  are  prescribed  by  the 
compatibility  requirements  of  the  equations  of  state,  as  in 
example  2.3.6  b,  only  (r  -  3)  of  the  model  properties  may 
be  chosen  arbitrarily. 

iv)  Calculate  the  correlation  functions  corre¬ 
sponding  to  the  r  dimensionally  independent  properties  of 
step  (iii)  above. 

v)  Find  the  laws  of  similitude  as  demonstrated 
in  the  example  of  Section  2.2.0. 

vi)  Find  the  value  of  the  remaining  model  prop¬ 
erties  from  the  laws  of  similitude  and  the  known  values  of 
the  correlation  functions  from  step  (iv). 

Method  II  -  Suppose  that  no  particular  material  has 
been  specified  for  the  model  study.  The  design  then  pro¬ 
ceeds  as  follows: 
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i)  Decide  which  fundamental  physical  quantities 
are  involved  in  the  experiment.  Let  their  number  be  r. 

ii)  Decide  which  variables  are  involved  in  the 
experiment. 

iii)  Choose  r  of  the  model  properties  arbitrarily, 
but  so  that  equation  (2,2.13)  is  satisfied.  These  quanti¬ 
ties  then  form  a  complete  dimensionally  independent  set 

for  the  experiment. 

iv)  Calculate  the  correlation  functions  corre¬ 
sponding  to  the  quantities  in  step  (iii). 

v)  Find  the  laws  of  similitude  as  demonstrated 
in  the  example  of  Section  2.2.8. 

vi)  Find  the  value  of  the  remaining  model  proper¬ 
ties  from  the  laws  of  similitude  and  the  known  values  of 
the  correlation  functions  from  step  (iv). 

vii)  Find  the  equation  of  state  for  the  model 
material  from  the  known  values  of  f^.,  ft  and  fT  and  the 
known  equation  of  state  of  the  prototype  material,  as  demon¬ 
strated  by  example  2.3.6  (a).  The  method  usually  breaks 
down  at  this  stage  because  a  model  material  with  the  required 
equation  of  state  cannot  be  found. 


It  should  be  noted  that  dimensional  analysis  can  be  used 
in  case  studies  only,  i.e.,  it  can  be  used  to  obtain  equations 
(2.1.11).  The  example  of  Section  2.2.8  will  be  used  to 
demonstrate  that  dimensional  analysis  cannot  be  used  directly 
to  obtain  equations  (2.1.10),  i.e.,  the  general  solution 
to  the  problem.  Equations  (2.2.22)  show  that 


But  it  will  be  equally  correct  to  use  the  remaining  equations 
of  (2.2.22)  to  write 
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where  a,  to  a4  are  constants.  In 
the  sane  value  as  ^  another  factor 


addition. 

(ir 


since  fp  has 
could  be 


multiplied  into  the  right  hand  side.  The  value  of  the 
constants  a^  to  a^  cannot  be  found  by  dimensional  analysis 
since  the  value  of  each  of  the  brackets  according  to  equa¬ 
tions  (2. 2. 22)cannot  be  different  from  one.  The  model  can, 
however,  be  used  to  find  the  relationship 


O'  =  F  (G,  r,  t,  a,  g,  x,  c,  T,  h,  £)  (2.4.2) 

This  is  achieved  by  varying  the  terms  on  the  right  hand  side 
one  by  one,  keeping  the  others  constant,  and  noting  the 
corresponding  changes  in  O’  .  If  expression  (2.4. 2)  is  in 
the  form  of  a  product,  it  can  be  used  to  determine  a4  to 
ag  in  equation  2.4.1  by  forming  the  quotient 

f  -  ®JL_  .  F(Qm,. . . .  ,Pm)  ^ 

*'  P(Gp . ,Pp) 


2.5  MODEL  ANALYSIS  AS  A  RANDOM  PROCESS 

In  the  derivation  of  the  Laws  of  Similitude  it 
was  assumed  tacitly  that 

i)  physical  phenomena  in  the  prototype  can  be  repre¬ 
sented  by  a  "certain"  scheme  of  knowledge,  I.e.,  each  prop¬ 
erty  can  be  described  by  a  fixed  number  at  a  particular 
point  in  space  at  a  particular  time, 

ii)  physical  phenomena  in  the  model  can  be  repre¬ 
sented  by  a  "certain"  scheme  of  knowledge, 

ill)  the  model  to  prototype  relationships  are  "certain" 
numbers. 

None  of  these  assumptions  hold  in  actual  situations  since 
such  random  processes  as  manufacture  and  loading  of  model 
and  prototype  are  always  involved.  It  therefore  becomes 
necessary  to  use  the  concepts  of  mean,  variance,  standard 
deviation  and  statistical  independence  as  defined  in  the 
Theory  of  Statistics,  Ref.  16.  Instead  of  equations  (.2.1.10) 
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and  (2.1.11)  it  becomes  necessary  to  find 

i)  the  variance  of  the  mean  of  the  prototype  out¬ 
come  as  predicted  from  the  mean  outcome  of  the  model,  its 
variance  and  the  variance  of  the  correlation  functions. 

ii)  the  standard  deviation  of  the  prototype  outcome 
as  predicted  from  the  standard  deviation  of  the  model 
outcome  and  the  correlation  functions. 

While  a  statistical  investigation  is  very  essential  to  the 
usefulness  and  reliability  of  a  model  test,  particularly 
in  cases  such  as  blast  loading  where  errors  may  occur  in 
the  modelling  of  a  large  number  of  properties,  investiga¬ 
tions  so  far  have  been  very  limited  due  to  the  great  number 
of  difficulties  encountered.  This  section  will  be  restricted 
to  a  presentation  of  basic  principles  and  the  outline  of 
a  method  of  solution. 

2.5.1  Determination  of  Stress,  Strain  or  Deflection 
Distribution.  The  deviation  of  the  actual  prototype  behav¬ 
ior  from  that  predicted  by  a  model  study  will  be  due  to 
a  combination  of  the  deviations  of  the  assumed  model  and 
prototype  properties. 

i)  Deviations  in  Geometry  -  The  position  of  appoint 
may  deviate  from  its  mean  position  in  any  direction  in  3- 
dlmenslonal  space.  This  can  be  represented  by  a  3-dimen¬ 
sional  density  function  d(ir,  ♦,  f  )  which  may  vary  from 
point  to  point  in  the  structure,  i.e. 

d(r,e,^)  =  [d(r,  ♦  ,f)]  (r,*,$)  (2.5.1) 

Whether  or  not  r,+and  f  are  statistically  Independent 
variables  depends  on  the  method  of  manufacture.  If  the 
method  of  manufacture  is  such  that  r,  ♦,  and  f  for  each 
point  are  determined  Independently,  equation  (2.5.1)  can 
be  rewritten 

[<KT.M)](wr)  (2.5.2) 

The  notation  on  the  right  hand  side  Indicates  that  d^(r), 
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d2(e)  and  d^(f )  may  all  vary  from  point  to  point  in  the  struc^ 
ture. 

li)  Devlationa  in  Material  Properties  -  Consider  any 
material  property,  denoted  ty  A.  Both  the  mean  and  the 
variance  of  A  may  be  different  for  different  points  on 
the  structure  so  that 

d(  A)  =  [d(A)]  (r,  ♦  ,  <p  )  (2.5.3) 

Both  the  nature  of  density  function  d(A)  and  its  change 
with  r,  ♦  and  4  depend  on  the  method  of  manufacture  of 
the  material  and  the  method  of  fabrication. of  the  struc¬ 
ture.  So,  for  instance,  the  variation  in  material  prop¬ 
erties  of  a  welded  steel  structure  will  depend  on  both  the 
quality  control  of  the  steel  mill  and  the  workmanship  of 
the  individual  welds. 

iii)  Deviations  in  External  Influences  -  The  uncer¬ 
tainty  In  the  exact  nature  and  amount  of  such  external 
influences  as  earthquakes  or  blast  loading  is  very  large 
since,  contrary  to  cases  (i)  and  (ii)  above,  there  is 
little  human  control  over  them.  Deviations  in  any  specific 
influence  B  may  again  be  different  from  point  to  point  in  the 
structure,  l.e. 

d(B)  =  [d(B)]  (r,  ♦ ,  <f  )  (2.5.1l) 

iv)  Deviations  in  the  Equation  of  State  -  Any  of  the 
partial  derivatives  used  in  Section  2.3  to  specify  the 
equation  of  state  may  deviate  from  its  assumed  value. 

This  can  again  be  represented  by  a  density  function 

d(C)  =  [d(C)]  <r,  +,<f  )  (2.5.5) 

v)  Deviations  in  the  Correlation  Functions  -  The 
necessity  of  using  correlation  functions  rather  than  ratios 
to  define  the  relationship  between  model  and  prototype  can 
now  be  demonstrated.  Assume,  for  Instance,  that  the  material 
property  A  has  a  normal  distribution  in  both  model  and 
prototype.  Then  its  density  function  is 
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COMPARISON  OF  THE  DISTRIBUTION  OF  A  MODEL  PROPERTY  <*_ 

A  PROTOTYPE  PROPERTY  c\p  AND  THEIR  CORRELATION  FUNCTION  f* 
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COMPARISON  OP  THE  DISTRIBUTION  OP  A  MODEL  PROPRRTV 
A  PROTOTYPE  PROPERTY  Cj^  AND  THEIR  CORRELATION  PUNCTION  f* 
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(2.5.6) 
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On  account  of  equation  (2.5.3)  the  standard  deviation 
and  the  mean  will  be  functions  of  r,#  and  f  in  both 
model  and  prototype.  The  correlation  function 


will  thus  also  have  a  mean  and  a  standard  deviation  depend¬ 
ing  on  r,«,9so  that  f  Itself  must  in  fact  in  the  most 
general  case  be  a  function  of  r,*  and 

The  density  function  for  f^  will  not  be  geomet¬ 
rically  similar  to  that  of  A.  Assume,  for  Instance,  that 
material  property  A  has  mean  values  of  I4.  and  1  in  proto¬ 
type  and  model  respectively.  Let  the  density  functions 
have  a  triangular  graph  and  maximum  deviations  of  +  10 % 
in  the  model  and  ±  in  the  prototype,  as  shown  in  Pig. 
2.5.1.  The  density  function  for  f^  can  then  be  determined 
numerically  as  shown  in  Tables  2.5.1  and  2.5.2. 


Table  2.5.1  - 


The  table  gives  values  of 


iisisehib 

0.0312 

0.0935 

0.156 

0.219 

0.219 

0.156 

0.0935 

0.0312 

kftfl 

ESS 

3.65 

3.75 

3.85 

3.95 

4.05 

cm 

ES9 

Umm 

0.055 

0.875 

4.16 

4.29 

4.40 

4.51 

4.63 

4.75 

4.86 

4.97 

0.167 

0.925 

3.95 

4.06 

4.16 

4.27 

4.38 

4.49 

4.59 

4.70 

0.278 

0.975 

3.75 

3.85 

3.95 

4.05 

4.15 

4.26 

4.36 

4.46 

0.278 

1.025 

3.56 

3.66 

3.76 

3.86 

-3.95 

4.05 

4.15 

4.25 

O.167 

1.075 

3.40 

3.49 

3.58 

3.68 

3.77 

3.86 

3.96 

4.05 

0.055 

1.125 

3.25 

3.34 

3.43 

3.52 

3.60 

3.69 

3.78 

3.87 

d(A)  x  A  A  ia  the  probability  that  A  lies  between  (A-fcAA)  and 
(A  +  h A  A),  where  AAp  =  0.1, A  A^  =  0.05.  Subscript  m  denotes 
model,  p  denotes  prototype. 
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Table  2.5.2  -  Calculation  of  d(f~  )  for  Various  Values  of 
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0.73 2  0.262  0.124  0.035 


In  each  column, 

d({A)  w4  -  Z  el(A,)  AAp-dC^AA^-IO4, 

Corrected  values  have  been  obtained  from  the  previous  row 
by  direct  scaling  so  that 

Z  et( fa)  6^  to4  »  104 

as  required  by  theory  of  statistics.  Then 

dL(-f^)  •  *  CORRECTED  VALUE 

The  last  row  is  used  to  obtain  the  mean  for  f^\ 

yU|A  -  Z  CcL(C)aC -4  00 

The  graph  of  the  density  function  d(f^)  is  shown  in  Pig. 
2.5.1.  It  may  be  concluded  that 

a)  The  mean  of  f"1  equals  the  quotient  of  the 
means  of  the  model  and  prototype  properties. 

b)  The  maximum  deviations  of  the  correlation 
function  f^  are  considerably  larger  than  those  of  either  the 
model  or  the  prototype  property:  +  28 %  and  -  22$. 

c)  The  distribution  of  the  correlation  function 
is  skew  whereas  that  of  the  properties  themselves  was  sym¬ 
metric. 

vi)  The  Statistical  Model  Investigation  -  The  statistical 
model  investigation  proceeds  as  follows: 

a)  Density  functions  for  the  geometric  and 
material  properties,  the  external  influences  and  the  equa¬ 
tions  of  state  are  either  assumed  or  determined  experiment¬ 
ally.  They  will  depend  on  the  particular  manufacturing  pro¬ 
cess,  loading  process,  etc.  Little  work  appears  to  have 
been  done  in  this  direction  up  to  the  present. 

b)  The  density  functions  for  the  correlation 
functions  of  the  properties  in  (a)  are  determined  as  demon¬ 
strated  in  Tables  2.5.1  and  2.5.2.  The  amount  of  computa¬ 
tion  involved  is  large,  and  will  increase  if  d(Am)  and  d(Ap) 
change  from  point  to  point  in  the  structure. 

c)  The  model  output  is  measured  a  consider¬ 
able  number  of  times  so  that  its  density  function  d(Dm)  can 
be  found. 


d)  The  density  function  d(fc)  must  be  found 
before  the  density  function  d(Dp)  of  the  prototype  output 
Dp  can  be  determined.  But  this  requires  that  f^  be  known 
as  a  function  of  all  the  other  correlation  functions  of  the 
experiment,  since  D  may  be  statistically  dependent  on  all 
the  other  structural  properties,  i.e.,  variations  in  any 
of  the  structural  properties  may  cause  D  to  vary  statisti¬ 
cally.  We  come  to  the  important  conclusion  that  statisti¬ 
cal  analysis  demands  a  knowledge  of  the  general  equations 
(2.1.11)  governing  the  phenomenon:  a  esse  study  alone  is  not 
sufficient.  Methods  to  obtain  equations  (2.1.11)  have 
already  been  discussed  in  Section  2,k»  Suppose  that  it  has 
been  found  that 


x  -  xb'  Xb*rb» .  .  rb» 

tr  ts  tr  +T  **  (2.5.7) 

If  it  is  assumed  that  the  variables  on  the  right  are  statis¬ 
tically  independent,  the  density  function  for  f  can  be- 
found  from 

(2.5.8) 


If  the  right  hand  members  of  equation  (2.5.8)  vary  from 
point  to  point  on  the  structure,  d(fr)  has  to  be  obtained 
for  each  point  individually.  The  amount  of  numerical  work 
involved  makes  the  use  of  a  computer  desirable. 

e)  Once  d(f0.)  is  known,  it  is  possible  to 
determine  d(o*p)  in  the  manner  of  Tables  2.5.1  and  2.5.2, 
using  the  relationship 


dK)  *  (2.5.9) 

vii)  Conclusions  - 

a)  If  the  variables  involved  in  the  experiment 
are  statistically  independent,  the  mean  of  the  prototype 
outcome  will  have  the  same  value  as  the  outcome  when  a  "cer¬ 
tain"  scheme  of  knowledge  is  assumed. 

b)  Since  the  correlation  functions  tend  to  have 
a  broad  distribution  as  shown  in  the  example  of  Pig.  2.5.1, 


the  prototype  outcome  will  have  a  much,  larger  deviation 
than  the  model  outcome. 

c)  Conclusions  (a)  and  (b)  prove  the  impor¬ 
tance  of  taking  a  large  number  of  readings  of  the  same 
outcome  in  the  model  experiment.  Only  in  this  manner  can 
the  effect  of  the  wide  distribution  of  the  correlation 
functions  be  eliminated  through  the  determination  of  a 
reliable  mean  value. 

2.5.2  Ultimate  Load  Models.  Perry  Borges,  Ref.  11,  has 
investigated  the  effect  of  f  on  the  distribution  of  the 
prototype  results  as  predicted  from  ultimate  load  model 
studies.  His  conclusions  are 

a)  An  increase  in  dimensions  considerably 
decreases  dispersions,  regardless  of  the  failure  criterion 
adopted. 

b)  The  transformation  of  the  mean  values  depends 
largely  on  the  failure  criterion  adopted.  In  brittle  rup¬ 
ture,  the  mean  values  decrease  as  the  dimensions  increase. 

In  ductile  rupture,  the  behavior  of  the  mean  values  depends 
on  the  number  of  failure  surfaces.  Finally,  in  the  theory 
of  similitude  for  failure  by  deformation  the  mean  values 
remain  constant  with  changes  of  scale.  The  latter  result 
agrees  with  the  deduction  in  Section  2.5.1. 

Since  there  are  factors  besides  f  which 
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affect  the  prediction  of  the  prototype  results  from  the 
model  study,  this  investigation  is  by  no  means  complete. 

2.6  STATIC  AND  DYNAMIC  STUDIES  OP  MASSIVE  STRUCTURES 

In  this  section,  the  theory  derived  in  preceding 
sections  will  be  applied  specifically  to  structures  whose 
mass  plays  a  significant  part  in  their  behavior. 

2.6,1  The  Effect  of  Mass  on  Structural  Behavior. 

The  mass  of  a  structure  affects  its  behavior  in  two  different 
ways : 


a)  It  causes  body  forces  due  to  the  earth’s  attrac¬ 
tion.  The  concept  of  an  "earth  acceleration"  causing  these 
forces  is  conventional  and  frequently  convenient,  but  it 
may  give  the  erroneous  impression  that  model  studies  involv¬ 
ing  earth  attraction  are  dynamic  in  character.  It  is 
preferable  to  base  similitude  of  model  and  prototype  on 
Newton’s  Law  of  Mass  Attraction: 


where  G  =  force  of  attraction  between  masses  m  and  M 
r  =  distance  between  masses 
'f  =  gravitational  constant 

Since  in  most  cases  the  earth’s  mass  M  as  well  as y and  r, 
the  distance  from  the  object  to  the  earth’s  center,  are  the 
same  for  model  and  prototype,  the  ratio  between  the  earth 
attraction  on  them  is  entirely  determined  by  the  ratio  of 
their  total  masses.  Earth  attraction  must  therefore  be 
considered  a  static  phenomenon. 

b)  The  mass  of  a  structure  or  of  its  loads 
also  gives  rise  to  inertia  forces  once  the  structure  as 
a  whole  or  parts  of  it  are  subjected  to  accelerations. 

These  accelerations  may  be  caused  by  dynamic  surface  loads, 
support  accelerations  or  free  vibrations  of  the  structure. 
According  to  Newton’s  Second  Law 

P  =  ma  =  -  Fi  (2.6.2) 

where  F  =  force  applied  to  particle  mass  m 

a  =  acceleration  of  m  due  to  P 

P,  =  inertia  force,  equal  in  magnitude  but 
opposite  in  direction  to  P. 

Since  the  acceleration  a  is  the  second  derivative  of  dis¬ 
tance  with  time,  inertia  forces  are  a  dynamic  phenomenon. 
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2.6.2  Types  of  Model  Studies  Involving  Mass.  The 
physical  laws  outlined  in  the  preceding  section  make  the  sub¬ 
division  of  model  studies  involving  mass  into  two  categories 
natural.  Each  of  the  two  categories  will  be  treated  sepa¬ 
rately  in  more  detail  in  later  sections. 

i)  Static  Studies  Involving  Mass  -  The  body  forces 
are  caused  by  earth  attraction  only,  all  surface  loads  are 
static  and  the  supports  do  not  move.  Typical  examples  are 

a)  Massive  concrete  structures  such  as  dams. 

b)  Structures  designed  mainly  to  carry  their 
own  dead  weight,  e.g.  large  span  bridges. 

c)  Structures  subjected  to  loads  caused  by 
heavy  masses,  e.g.  grain  silos  and  coal  bunkers. 

ii)  Dynamic  Studies  Involving  Mass,  If  the  prototype 
structure  is  subjected  to  dynamic  surface  forces,  support 
accelerations  or  free  vibrations,  a  dynamic  model  study  with 
proper  regard  to  inertia  forces  has  to  be  undertaken. 

Dynamic  studies  may,  or  may  not,  irtclude  also  the  effects 

of  static  surface  loads  and  earth  attraction  acting  on  the 
prototype,  depending  on  the  relative  size  of  the  forces 
involved.  Earth  attraction  can, for  instance,  frequently 
be  disregarded.  Since  concrete,  for  example,  has  a  mass 
attraction  per  cubic  foot  roughly  equivalent  to  one  psi 
acting  on  the  surface,  the  dead  weight  stresses  of  say  12” 
thick  shells  subjected  to  blast  load  pressures  of  25  to  200 
psi  can  safely  be  neglected. 

Typical  examples  of  dynamic  studies  are 

a)  Surface  structures  such  as  dams  or  shells 
responding  to  dynamic  surface  loads. 

b)  Articulated  structures  such  as  bridges,  trans¬ 
mission  towers  or  radar  antennae  subjected  to  drag  loads 
caused  by  atomic  explosions. 

c)  Structures  supporting  heavy  masses  and  sub¬ 
jected  to  either  blast  loads  or  seismic  support  motion. 

As  in  the  case  of  model  studies  where  mass  is  not  of  impor¬ 
tance,  the  following  types  of  structural  behavior  may  be 
investigated: 


i)  Stress  or  strain  distributions  and  defor¬ 
mations  if  a  structure  is  linearly  elastic  (i.e.,  super¬ 
position  valid  if  deformations  are  small)  or  if  the  struc¬ 
ture  is  nonlinearly  elastic  or  plastic  (superposition  not 
valid) . 

ii)  Ultimate  load  bearing  capacity  of  the  struc¬ 
ture. 

iii)  Determination  of  the  loading  causing  elastic 
or  inelastic  structural  instability. 

In  addition,  however,  the  gravity  and  inertia  forces  on 
masses  supported  by  the  structure  should  be  faithfully 
reproduced  in  the  model.  Both  the  effect  of  their  inertia 
forces  on  the  overall  structural  action  and  their  destruc¬ 
tive  action  as  missiles  are  of  importance.  It  is  thus 
important  to  reproduce  the  interaction  between  the  masses 
themselves,  such  as  friction  forces,  as  well  as  the  inter¬ 
action  between  the  masses  and  the  structure. 

2.6.3  The  Laws  of  Similitude  for  Massive  Structures. 

It  will  be  assumed  in  this  section  that  the  correlation 
functions  can  be  chosen  independently  of  the  model  material, 
i.e.,  method  II  of  Section  2.I4.  is  used.  The  following 
relationships  between  the  correlation  functions  have  already 
been  derived  in  Section  2.2.8: 


Accelerations : 

11 

CB 

*-l 

v;2 

(2.6.3) 

Body  Forces: 

II 

c!J 

fgfa 

f  f  fl^ 
g  r  t 

(2.6.4) 

Surface  Forces: 

*-> 

«* 

II 

Vr  : 

V.V 

= 

Vrft2 

(2.6.5) 

These  equations  are  valid  for  both  static  and  dynamic  studies. 
Generally  f  has  to  have  a  value  considerably  less  than  unity 
to  make  a  model  study  possible.  This  can  be  achieved  by 
varying  either  tne  parameter  f  which  is  dependent  on  the 
external  influences  or  the  parameters  f ^  and  fg  which  are 
dependent  on  the  model  material. 
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The  difference  between  the  laws  of  similitude 
for  static  and  dynamic  studies  lies  in  the  freedom  that  may 
be  exercised  in  the  choice  of  the  correlation  functions. 

In  a  static  study  where  no  time  measurements  need  to  be 
carried  out  on  the  model,  equation  (2.6.3)  is  automatically 
satisfied  since  ffc  can  be  given  an  arbitrary  value.  In  a 
dynamic  model  study,  the  time  correlation  function  f^  is 
of  much  more  importance  since  the  prototype  load-time 
function  and  support  acceleration-time  function  must  be 
scaled  for  the  model  according  to  f^..  The  choice  of  f^. 
depends  largely  upon  the  nature  of  the  dynamic  load.  It 
is  convenient  to  distinguish  between  the  following  types: 

Load  Type  I  -  The  applied  surface  forces  are  known 
as  a  function  of  both  space  and  time.  If  it  is  assumed 
that  the  technical  problems  which  arise,  e.g.  in  the  model¬ 
ing  of  load  rise  times  of  about  2  ndlli seconds ,  can  be 
solved,  this  type  of  dynamic  load  imposes  no  restrictions 
on  the  correlation  functions. 

Load  Type  II  -  The  applied  surface  forces  are 
not  known  as  a  function  of  both  space  and  time,  but  are 
caused  by  phenomena  subject  to  known  physical  laws,  e.g. 
those  of  aerodynamics.  The  model  must  then  be  subjected 
to  the  actual  physical  event  that  causes  the  loading,  but 
the  variables  involved  (e.g.  overpressure)  may  be  scaled 
down  according  to  the  chosen  values  of  the  correlation 
functions.  The  laws  of  similitude  must  now  cover  not  only 
the  model  itself,  but  also  the  events  leading  to  the  load¬ 
ing.  In  many  cases  it  will  still  be  possible  to  chose 
f,.  and  ft  different  from  unity. 

Load  Type  III  -  The  dynamic  surface  forces  and 
the  physical  laws  of  the  events  by  which  they  are  caused, 
are  not  exactly  known,  as  in  some  types  of  blast  loading. 

The  model  must  then  be  subjected  to  the  actual  prototype 
loading  which  implies  that  f^  =  fr  =  1,  so  that  the  choice 
of  model  materials  now  becomes  very  limited.  There  is, 
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in  addition,  no  guarantee  that  f  is  not  also  one  of  the 
parameters  affecting  the  space-time  distribution  of  the 
loading,  so  that  in  effect  structures  subjected  to  type 
III  loads  cannot  always  be  modelled  with  confidence. 

2.6.4  Static  Model  Studies  with  Natural  Earth  Attrac 
tlon.  According  to  Newton's  Gravitational  Law,  (2.6. l), 
applied  to  the  earth  (massM)  and  model  or  prototype 
respectively,  the  body  forces  are 


Gm 

°P 

so  that 

.  °n,  . 

g_ 

ffl  _  X. 

fG 

GP 

*p  g 

Substitution  into  equation  (2.6.5)  yields 


(2.6.6) 

(2.6.7) 


f<r=  fgfr  (2.6.8) 

The  value  of  f  can  generally  be  varied  within  very  narrow 
limits  only,  as  is  shown  in  Table  2.6.1, 


Table  2,6.1  -  Typical  Average  Values  of  Density  Ratio  f  . 

O 


— 

Prototype  Material 

Model 

Steel 

Aluminum 

Concrete 

Masonry 

Material 

Spec.  Mass 

7.75 

2.65 

2.40 

2.00 

Mortar 

0.213 

0.621 

0.687 

0.825 

Gypsum 

0.296 

0.869 

0.959 

1.150 

Plexiglas 

0.165 

0.454 

0.500 

0.600 

Rubber 

1.10 

0.142 

0.415 

0.458 

0.550 

Steel 

7.75 

1.000 

2.920 

3.230 

3.875 

Aluminum 

2.65 

0.341 

1.000 

1.100 

1.325 

Ph.  Bronze 

8.90 

1.150 

3.360 

3.710 

4.450 

The  table  shows  that  the  density  ratio  is  less  than  5 
whereas  the  length  ratio  is  normally  less  than  0.125. 

It  therefore  follows  from  equation  (2.6.8)  that  the  stress 
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ratio  fy  is  less  than  1.  But  according  to  the  second  of 
equations  (2.3.10),  the  stress  ratio  for  linearly  elastic 
model  and  prototype  materials  is 


P 


Thus  the  model  material  must  be  less  rigid  than -the  proto¬ 
type  material.  It  is  suggested  that  instead  of  choosing 
f  and  then  searching  for  a  model  material  which  is  such 
that  and  f  satisfy  equation  (2.6.8),  the  inverse  pro- 

o 

cedure  is  followed,  as  demonstrated  in  the  following  example. 


Example:  A  plexiglas  model  of  a  steel  bridge  is  tested 
in  the  linear  elastic  range.  If  it  is  assumed 
that  E  =  lj.50,000  psi  for  plexiglas 

E  =  30  x  10^  psi  for  steel, 
find  the  length  correlation  function  fr  required 
so  that  static  body  forces  are  correctly  repro¬ 
duced  in  the  model. 

Solution:  Prom  equation  (2.6.9): 


f  =  o.fr-5  *  106  . 

r  30  x  106 

From  Table  2.6.1:  f  = 

g 

Prom  equation  (2.6.8):  f 


0.015 

0.165 

=  0*015  =  1_ 
0.165  11 


(Ans. ) 


This  example  demonstrates  that  unless  f  or  f  are  artifi- 

g  8 

daily  increased,  relatively  large  models  are  required  to 
reproduce  gravity  stresses  faithfully.  Where  very  flexible 
materials  are  used,  e.g.,  rubber  models  of  concrete  dams^^,, 
smaller  models  are  feasible,  but  the  high  Poisson  Ratio  of 
the  model  material  frequently  makes  the  test  results  unreli¬ 
able  (see  Section  2.3.8). 
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2.6.5  Technical  Aspects  of  Static  Studies  with  Natural 
Earth  Attraction.  If  the  model  study  Is  concerned  with  the 
ultimate  load  bearing  capacity  or  with  the  stability  of  a 
structure,  direct  observation  without  the  use  of  sensors 
and  recording  devices  for  strains  and  deflections  will  be 
sufficient,  so  that  both  the  problems  of  body  force  load¬ 
ing  and  of  model  observation  are  automatically  solved. 

If  the  stress  distribution  in  the  structure  is  to  be  deter¬ 
mined,  considerable  difficulties  both  in  the  determination 
of  the  strain  distribution  and  the  derivation  of  the  corre¬ 
sponding  stress  distribution  are  experienced. 

i)  Strain  Distribution  -  The  strains  caused  by  the  desd 
weight  of  a  model  are  extremely  small.  Consider  for  instance 
the  extremely  favorable  case  of  a  mortar  model  (E  =  75,000  psi, 
w  =  100  pcf)  of  a  500*  nigh  concrete  wall  with  f  =  0.01. 

Then  the  maximum  strain  at  the  foot  of  the  model  is  of  the 
order 


€  = 


5  x  loo 
1^4-  *  75,000 


x  10 


6  gi 


50  ^in/in. 


Not  even  the  highly  developed  electrical  strain  measuring 
devices  are  capable  of  reliable  detection  at  such  low  strain 
levels.  Another  problem  is  the  determination  of  a  zero 
reading,  I.e.,  a  strain  gage  reading  for  zero  body  forces. 

To  overcome  this  serious  problem  and  at  the  same  time  to 
increase  the  strain  levels,  tests  have  been  carried  out 
where  the  strains  ere  first  measured  with  the  model  in  an 
upright  position,  and  then  with  the  model  inverted ^ . 

In  this  manner,  strain  readings  are  doubled.  In  addition 
to  the  mechanical  problems  caused  by  the  rotation  of  large 
models  about  a  horizontal  axis,  the  method  possesses  an 
inherent  inaccuracy  where  the  stress-strain  properties  of  the 
model  are  different  in  tension  and  in  compression.  In  most 
cases,  the  strain  level  remains  too  low  for  accurate  sensing 
and  recording.  Additional  difficulties  are  encountered  If 
strains  Inside  the  model  are  to  be  measured.  Since  it  is 


53 


2.6.5  Technical  Aspects  of  Static  Studies  with  Natural 
Earth  Attraction.  If  the  model  study  Is  concerned  with  the 
ultimate  load  bearing  capacity  or  with  the  stability  of  a 
structure,  direct  observation  without  the  use  of  sensors 
and  recording  devices  for  strains  and  deflections  will  be 
sufficient,  so  that  both  the  problems  of  body  force  load¬ 
ing  and  of  model  observation  are  automatically  solved. 

If  the  stress  distribution  in  the  structure  is  to  be  deter¬ 
mined,  considerable  difficulties  both  In  the  determination 
of  the  strain  distribution  and  the  derivation  of  the  corre¬ 
sponding  stress  distribution  are  experienced. 

1)  Strain  Distribution  -  The  strains  caused  by  the  dead 
weight  of  a  model  are  extremely  small.  Consider  for  instance 
the  extremely  favorable  case  of  a  mortar  model  (E  =  75,000  psi, 
w  =  100  pcf)  of  a  500'  nigh  concrete  wall  with  f  =  0.01. 

Then  the  maximum  strain  at  the  foot  of  the  model  is  of  the 
order 

6  =  JL  £-129 - x  106  ~  50/iin/in. 

liOi  x  75,000  7 

Not  even  the  highly  developed  electrical  strain  measuring 
devices  are  capable  of  reliable  detection  at  such  low  strain 
levels.  Another  problem  is  the  determination  of  a  zero 
reading,  i.e.,  a  strain  gage  reading  for  zero  body  forces. 

To  overcome  this  serious  problem  and  at  the  same  time  to 
increase  the  strain  levels,  tests  have  been  carried  out 
where  the  strains  are  first  measured  with  the  model  in  an 
upright  position,  and  then  with  the  model  inverted 
In  this  manner,  strain  readings  are  doubled.  In  addition 
to  the  mechanical  problems  caused  by  the  rotation  of  large 
models  about  a  horizontal  axis,  the  method  possesses  an 
inherent  inaccuracy  where  the  stress-strain  properties  of  the 
model  are  different  in  tension  and  in  compression.  In  most 
cases,  the  strain  level  remains  too  low  for  accurate  sensing 
and  recording.  Additional  difficulties  are  encountered  if 
strains  inside  the  model  are  to  be  measured.  Since  it  is 
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seldom  possible  to  make  the  sensing  device  out  of  the  model 
material,  the  presence  of  the  device  will  cause  discontinu¬ 
ities  and  changes  In  the  strain  field  which  frequently  may 
lead  to  erroneous  results. 

11)  The  problems  encountered  in  deriving  the  stress 
distribution  of  a  3-dimensional  model  from  Its  strain  dis¬ 
tribution  are  of  a  technical  nature.  They  have  already 
been  discussed  in  Section  2,3.8, 

2.6.6  Static  Model  Studies  with  Earth  Attraction 
Artificially  Induced  by  Acceleration.  The  large  models 
required  by  the  method  of  Section  2.6.4  and  the  unsatisfac¬ 
torily  low  strain  levels  obtained  lead  to  the  use  of  artifi¬ 
cial  methods  of  simulating  earth  attraction.  Consider 
equation  (2.6.5)! 

f  =  f  f  f 
w  gar 

Suppose  that  in  a  particular  model  study  it  is  considered 
essential  to  reproduce  both  the  stress-strain  character¬ 
istics  and  shrinkage  and  creep  properties  of  the  prototype 
material,  and  that  this  can  be  achieved  only  by  using  the 
prototype  material  as  model  material.  Then 


f  „  =  f„ 
v  8 


and  equation  (2.6.5) 


f  n  f  _  =  1 

a  r 


=  1 

becomes 


(2.6.10) 


The  model  must  therefore  be  subjected  to  very  considerable 
accelerations,  to  which  it  will  respond  dynamically.  In 
order  to  demonstrate  the  implications  of  the  dynamic  response, 
a  typical  though  simplified  example  is  given:  consider  a 
3*0"  x  4"  x  1"  mortar  beam  (E  =  3  *  10^  psi,  w  =  100  pcf) 
which  is  simply  supported  on  a  sleigh  as  shown  in  Pig.  2.6.1. 
If  the  sleigh  is  given  a  constant  acceleration  i  along 
a  straight  line,  the  beam  is  subjected  to  constant  forces: 

P  =  *3"zL 
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MID5PAN  MOW>ENT(%  OF  "srA-TIc" MOMENT) 


FIGURE  2.6.1 


SLEIGH  MOVES  WITH  CONSTANT 
ACCELERATION  IN  THIS  DIRECTION. 


PLAN  OP  SLEIGH  CARRVINS  MORTAR  BEAM. 


o  6-i  0-4-  0-4  0-0  i-O  t-Z  t-4- 


TIME  (EEC) 

FIGURE  2.6.  2. 

THE  GRAPH  EH  OWE  the  envelopes  ONIX  OP  THE  CURVES  DESCRIBING 
VARIATION  OP  MIPSPAM  BENDING  MOMENT  WITH  TIME  OF  THE  BEAM 
IN  FIS.  2.6.1 
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at  each  of  the  supports,  where 

^  =  mass  of  beam  per  unit  length 

L  =  length  of  beam. 

The  dynamic  response  of  the  beam  is  analysed  in  Appendix  A. 
The  actual  bending  moment  at  midspan  is  calculated  and  com¬ 
pared  to  the  static  bending  moment  caused  by  a  load  ^  2  per 
unit  length  of  the  beam.  Since  ^2  is  the  body  force  per 
unit  length  which  we  actually  wish  to  obtain,  and  which  would 
have  been  obtained  if  the  beam  were  perfectly  rigid,  the 
calculations  give  the  error  caused  by  the  dynamic  response 
of  the  beam. 

The  results  are  presented  in  Pig.  2.6.2.  It 
may  be  concluded  that 

i)  The  model  is  subjected  to  stresses  far  exceeding 
those  of  the  steady  state,  represented  by  100$  in  Pig.  2.6.2. 
The  method  is  thus  not  suited  to  ultimate  strength  tests. 

ii)  Models  frequently  have  less  than  5$  of  the  critical 
damping.  No  readings  should  therefore  be  taken  before  the 
model  has  been  accelerated  for  at  least  2  seconds.  If  2 
more  seconds  are  required  to  take  readings,  and  the  model  is 
subjected  to  20  x  earth  acceleration,  the  distance  travelled 
is 

s  =  i  at2  =  5150  ft. 

Since  the  model  has  to  be  decelerated  more  gradually,  the 
total  test  run  becomes  very  long.  Linear  acceleration  may 
be  eliminated  as  a  practical  method  of  model  testing. 


Accelerations  can  also  be  achieved  by  means  of 
centrifuges.  In  addition  to  ejconomic  problems  and  techni¬ 
cal  difficulties  such  as  strain  measurement'  on  fast  moving 
models,  the  method  possesses  an  inherent  source  of  error,  which 
will  be  demonstrated  by  means  of  an  example: 

Suppose  that  the  gravity  stress  distribution  in 
a  1200’  wide  and  I4.OO  ’  high  concrete  dam  is  to  be  determined. 

A  model  material  is  chosen  so  that  f  *  \  and  f _ =  0.2. 

g  r 
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FIGURE  2.6.3 

MODEL  SUBJECTED  TO  "GRAVITY  STRESSES*  IN  CENTRIFUGE. 


FIGURE  2.6.4 


FORCES  ACTING  BETWEEN 

load  ano 

STRUCTURE. 

N  NORMAL  FORCE 

/*• 

FRICTION  COEFFICIENT 

M  MASS  OF  LOAD 

6 

body  FORCE 

at  earth  acceleration 

A 

acceleration  OF  load 

LIMITING  CONDITION  FOR 

RELATIVE  MOTION  . 

6  “Ml 
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A  10*  0"  diameter  centrifuge  is  available  and  accuracy  in 
model  manufacture  and  strain  measurements  limits  f  to  1 
in  800.  Then 


A  particle  moving  along  a  circle  with  radius  r  at  a  constant 
angular  velocity  u>  is  subjected  to  a  constant  centripetal 
acceleration 

a  =  u>^r 

Thus  if  point  A  in  figure  2.6.3  is  to  have  an  acceleration  of 

p 

160  x  32.2  ft/sec  ,  the  frequency  of  the  centrifuge  must  be 

f  =  1  /»?.*  3212~,=  5.25 

2*V  4.75 

At  points  B  and  C  on  the  center  line,  the  acceleration 
will  be  about  5%  in  error,  but  will  act  in  the  correct 
direction.  At  D,  the  component  parallel  to  OC  will  be 
5#  in  error,  and  there  is  an  additional  acceleration  equal 
to  16$  of  the  acceleration  at  A  acting  in  direction  DB. 

The  inherent  errors  are  thus  not  excessive  and  the  method 
will  be  useful  if  the  technical  difficulties  can  be  overcome. 

2.6.7  Static  Model  Studies  with  Earth  Attraction 
Artificially  Induced  by  Surface  Forces.  The  earth  attrac¬ 
tion  on  the  prototype  material  can  also  be  simulated  by  the 
application  of  surface  loads  at  discrete  locations  on  the 
model.  The  method  has,  for  instance,  been  used  in  model 
studies  of  concrete  arch  dams  by  Rocha,  Serafim  and  Ferreira . 
If  the  forces  are  applied  to  the  face  of  the  model,  the 
continuity  of  the  model  material  need  not  be  disturbed, 
but  considerable  errors  in  the  stress  distribution  in  thick 
models  may  result.  If  they  are  applied  inside  the  model, 
the  stiffness  of  the  structure  is  affected.  A  quantitative 
analysis  of  these  problems  is  given  in  Section  2.7. 
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The  effect  of  the  surface  loads  on  the  model 
can  be  represented  by  assigning  the  model  material  an 
artificial  specific  mass  g^.  Suppose 


Then 


m 


w 


m 


V 


m 


®m 


If 


then 


total  surface  load  applied  to  entire  model, 
total  weight  of  model, 
total  volume  of  model, 
earth  acceleration 


P  +  W 
m  m  _ 

V  a 
m  e 


m 


V  a 
m  e 


+ 


(2.6.11) 


artificial  specific  mass  correlation  function 
specific  gravity  of  prototype 


5" 

Sr, 


m 


+  f. 


Vv 
m  p 


g 


(2.6.12) 


It  Is  now  possible  to  select  both  the  size  of  the  model 
(i.e.,  f  )  and  the  model  material  (l.e.  f^  and  f  )  then' 
to  determine  the  required  fg  from  equation  (2.6.5)  and 
finally  to  compute  Pm  from  equation  (2.6.12).  The  rules 
according  to  which  Pm  should  be  distributed  over  the  model 
will  be  discussed  in  Section  2.7. 

Artificial  increase  of  model  specific  mass  could 
also  be  achieved  by  making  the  model  of  a  magnetic  material 
and  placing  it  in  a  magnetic  field.  The  method  has  the 
basic  limitation  that  magnetic  lines  of  force  leave  a  surface 
normally,  i.e.,  the  direction  of  the  forces  would  not  agree 
with  the  direction  of  earth  attraction.  In  addition,  the 
magnetic  force  is  not  uniformly  distributed  over  the  thickness 
of  the  model. 

2.6.8  Dynamic  Model  Studies  with  Natural  Earth  Attraction. 
A  model  study  becomes  dynamic  as  soon  as  dynamic  surface 
forces  or  support  accelerations  must  be  reproduced.  The 
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structure  may,  or  may  not,  at  the  same  time  be  subjected 
to  static  surface  loads  and  earth  attraction.  It  now 
becomes  of  importance  to  differentiate  between  loads  hav¬ 
ing  significant  mass,  and  thus  subject  to  inertia  forces 
like  the  structure  itself,  and  loads  whose  force  only  needs 
to  be  simulated.  The  relationship  between  the  correlation 
functions  and  the  importance  of  f^.  have  already  been  dis¬ 
cussed  in  Section  2.6.3.  In  this  case,  as  in  Section  2.6.I4., 


f 


ff 


so  that  fj, 


f  f 
g  r 


(2.6.13) 


In  addition,  f  =  1 

so  that  f  =  f.^  (2.6.  II4.) 

I*  b 

Once  the  model  material  has  been  chosen,  fr  and  fg  are 
known  and  f  can  be  determined  from  equation  (2.6.13). 

Then  f^  is  given  by  equation  (2. 6. 1)4.).  Since  surface 
loads  must  be  scaled  according  to  the  stress  correlation 
function  fff  ,  the  load-time  curve  for  the  model  is  now 
uniquely  determined. 

This  method  therefore  permits  the  effect  of 
both  earth  attraction  and  dynamic  surface  loads  to  be 
exactly  reproduced  in  the  model,  but  it  is  tacitly  assumed 
that  the  loading  is  of  Type  I.  If  the  prototype  is  sub¬ 
jected  to  support  accelerations,  these  must  be  applied  to 
the  model  in  the  ratio  f  =  1.  The  technical  difficulties 
are  the  same  as  in  Section  2.6,5.  Deta  recording  will, 
however,  become  more  complicated  since  strain  versus  time 
curves  must  now  be  obtained. 

2.6.9  Dynamic  Model  Studies  with  Artificially  Induced 
Earth  Attraction.  In  many  model  studies,  the  previous  method 
may  not  be  suitable  because  the  models  become  too  large 
or  because  the  loading  is  not  of  Type  I.  It  will  then  be 
necessary  to  simulate  earth  attraction  artificially.  While 
the  earth  attraction  body  forces  must  be  reproduced  artificially 
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it  should  be  noted  that  the  inertia  body  forces  will  auto¬ 
matically  satisfy  the  laws  of  similitude  as  is  demonstrated 
by  the  following  example: 

Consider  a  point  load  acting  on  a  particle 
of  unit  volume  with  specific  mass  gp,  giving  it  an  acceler¬ 
ation  a  according  to  Newton's  Second  Law: 

P 

Pp  ~  gp8p  (2.6.15) 

A  model  of  the  particle  is  made  at  scale  f  and  the  surface 
forces  are  scaled  in  ratio  f^.  .  We  wi3h  to  prove  that 
the  model  is  automatically  subjected  to  the  correct  inertia 
body  forces  but  that  the  laws  of  similitude  for  earth  attrac¬ 
tion  are  not  satisfied  unless  f  =  1,  a  case  which  we  have 
excluded  in  this  section. 

The  model  properties  are 


®m  fg%p 

Vm  =  tl  (2.6.16) 

According  to  Newton's  Second  Law 


m 


=  V. 


m 


®mam 


Substituting  from  aquation  (2,6.16): 

r»frPp  1  rrfgVm 

Using  equation  (2.6.15): 


*m  ~  f  f  aP 

r  g 

•’*  f  ^  =  frfgfa  (2.6.17) 

Equations  (2.6.17)  and  (2.6.5)  are  identical  so  that  inertia 
forces  are  always  correctly  modeled. 
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The  earth  attraction  on  the  prototype  is 


A 


P 


(2.6.18) 


where  M  and  r  for  Newton's  Gravitational  Law  have  been 
defined  previously.  The  earth  attraction  on  the  model  is 

Mfrf ffgn 

Am  =  t  (2.6.19) 


But  since  A  is  a  point  force,  it  follows  that 
m 


m 


=  fff  f  LA 

o’  r  p 


(2.6.20) 


Equations  (2.6.18)  to  (2.6.20)  yield 

f<r  =  Vg  (2.6.21) 

Equation  (2.6.21)  is  a  special  case  of  equation  (2.6.5) 
which  proves  that  the  earth  attraction  is  modelled  correctly 
only  if  f  *  1. 

If  f  ^  1,  artificial  methods  can  be  used  to 
a 

reproduce  the  earth  attraction  body  forces.  It  is  con¬ 
sidered  impractical  to  apply  both  a  blast  load  and  an 
artificial  acceleration  to  the  model,  and  the  discussion 
will  therefore  be  restricted  to  the  use  of  surface  forces 
for  the  simulation  of  earth  attraction. 

As  shown  in  Section  2.6.7,  both  the  size  of  the 

model  (i.e.,  f  )  and  the  model  material  (i.e.,  f „  end  f  ) 
r  v  g 

may  be  chosen  arbitrarily  or  in  such  a  manner  that  they 
satisfy  the  conditions  imposed  by  the  dynamic  surface  loads, 

Type  II  or  III.  Then  f  is  determined  from  equation  (2.6.5) 

O 

and  Pm  from  equation  (2.6.12).  As  shown  previously,  the  inertia 
body  forces  due  to  the  dynamic  surface  loads  are  automatically 
reproduced  correctly  in  the  model.  The  inertia  body  forces 
due  to  support  accelerations  will  be  scaled  correctly  if 
the  support  accelerations  of  the  model  are  calculated 
using  f  as  determined  from  equation  (2.6.5). 
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The  problem  of  reproduction  of  gravity  stresses 
in  models  is  thus  basically  the  same  in  static  and  in 
dynamic  studies.  One  technical  problem  is,  however,  added 
in  the  latter  case  since  the  surface  loads  simulating  earth 
acceleration  have  to  remain  constant  when  the  model  vibrates. 

2.6.10  Modelling  of  Heavy  Masses  Supported  by  Structures. 
A  great  variety  of  important  public,  industrial  and  military 
structures  support  heavy  masses  e.g.,  warehouses,  grain 
silos,  oil  storage  tanks,  etc.  Since  the  dynamic  structural 
action  in  these  cases  is  largely  determined  by  the  supported 
masses,  their  effect  may  not  be  neglected  in  a  model  study. 

In  particular,  it  is  necessary  to  take  into  account: 

a)  the  earth  attraction  on  them. 

b)  inertia  forces  due  to  dynamic  response. 

c)  friction  forces  acting  between  individual  loads 
and  between  the  loads  and  the  structure. 

d)  their  destructive  action  if  dynamic  response  of 
the  structure,  or  direct  exposure  to  dynamic  external  loads, 
turns  them  into  projectiles. 

Factors  (a)  and  (b)  may  be  studied  in  the  same  manner  as 
the  dead  weight  and  inertia  forces  of  the  structure  itself, 
as  discussed  previously.  Factors  (c)  and  (d)  are  discussed 
individually  in  the  following  two  sections. 

i)  Modelling  of  Friction  Forces  -  Consider  a  proto¬ 
type  mass  Mp  supported  by  a  structure  so  that  the  coefficient 
of  friction  between  them  is yu.p.  The>maximum  acceleration 
8p  to  which  the  prototype  structure  may  be  subjected  before 
the  mass  starts  moving  relative  to  the  structure  can  be 
deduced  from  Fig.  2,6.14. 

aP  -/V. 

where  a  =  earth  acceleration, 
e 

Similarly  for  the  model 

a  =  M  a 
m  /me 
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so  that 


a  M. 

m  _  '  m 

ap  A  (2.6.22) 

As  discussed  previously,  accelerations  in  the  model  frequently 
ere  considerably  larger  than  In  the  prototype,  but  generally 
In  these  cases,  similarly  with  regard  to  friction 
forces  can  be  achieved  by  artificial  methods  only. 

a)  Increase  the  normal  force  N.  An  additional 
force  P  normal  to  the  surface  between  model  structure  and 
load  is  provided  so  that 

N  =  P+Ve 

and 

Vm  =A»(P+Ve> 

so  that 

*P  Lv. 

Once  ffl  and  are  known,  the  required  value  of  P  can  be 
calculated  from  equation  (2.6.23).  The  force  P  can  be 
obtained  by  using  magnets  as  loads  and  making  the  model  of 
a  magnetic  substance  where  the  loads  rest  on  it.  The 
method  has  the  advantage  that,  once  the  load  starts  moving 
relative  to  the  structure,  the  full  friction  force  will 
remain  acting  between  their  surfaces  and  the  load  will  again 
come  to  rest  relative  to  the  structure  once  the  inertia 
forces  ‘become  less  than  the  maximum  available  friction  force. 
This  is  of  Importance  where  a  structure  goes  through  several 
cycles  of  vibration  before  it  reaches  maximum  stress  or  some¬ 
times  ultimate  failure  conditions. 

b)  Increase  coefficient  of  friction.  The  mass 
representing  the  prototype  load  may  also  be  glued  to  the 
model,  particularly  if  it  is  reasonable  to  assume  that  the 
load  on  the  prototype  will  not  move  relative  to  the  struc¬ 
ture,  so  that  it  is  only  necessary  to  obtain  a  sufficiently 


- 


(2.6.23) 
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large  friction  force  In  the  model.  It  must,  however,  be 
remembered  that,  once  the  bond  i-s  broken,  the  behavior 
of  model  and  prototype  loads  Is  no  longer  similar. 

These  artificial  methods  will  not  be  necessary  where  natural 
earth  attraction  Is  employed  in  the  model  study,  l.e.,  ffl  =  1. 

11)  Modeling  of  Projectile  Action  -  Once  a  load  moves 
relative  to  the  structure,  It  becomes  a  projectile  which, 
upon  impact  with  parts  of  the  structure  or  with  its  contents, 
may  cause  considerable  damage.  True  modeling  of  this 
phenomenon  is  particularly  important  since  a  general  ana¬ 
lytical  treatment  of  the  problem  is  not  available  (see  Ref.  30, 
p.  12?  et  seq.)  The  following  limitations  are  imposed  upon 
the  correlation  functions  for  the  geometry  and  material 
properties  of  the  masses  representing  the  model  loads: 

a)  Since  the  size  of  the  impact  area  may  be 

of  importance,  the  model  load  must  be  geometrically  simi¬ 
lar  to  the  prototype  load,  and  fp  for  load  and  structure 
must  be  the  same. 

b)  While  the  load  is  still  moving  with  the  struc¬ 
ture,  it  is  subjected  to  the  same  acceleration  so  that  f 

is  the  same  for  both,  and  thus  ffc  must  also  be  the  same 
(since  f  is  already  fixed). 

c)  The  ratio  between  the  prototype  and  model 

p 

kinetic  energies  ( Jmv  )  has  to  be  the  same  for  structure 
and  load  so  that  in  addition  to  the  preceding,  f  =  1. 

d)  The  percentage  of  the  kinetic  energy  which 
must  be  absorbed  by  the  structure  depends  largely  upon  the 
nature  of  the  impact  if  the  projectile  mass  is  considerably 
smaller  than  that  of  the  structure  and  the  masses  still 
adhering  to  it.  It  is  considerably  less  for  plastic  than  for 
elastic  impacts.  For  this  phenomenon  to  be  correctly  repro¬ 
duced,  it  is  essential  that  is  the  same  for  structure  and 
load. 

It  may  thus  be  concluded  that,  for  missile  action  to  be 
correctly  modeled,  the  correlation  functions  fr,  ffc,  f 


and  «fr  have  to  be  identical  for  the  structure  and  its  loads. 
In  addition,  the  requirements  of  Section  (a)  have  to  be 
met.  These  conditions  are  not,  however,  sufficient  to 
guarantee  similitude  for  impulsive  action  since  factors 
such  as  the  crystalline  structure  of  the  prototype  struc¬ 
tural  material  and  load  may  be  of  importance. 

2.7  THE  SIMULATION  OP  BODY  FORCES  BY  SURFACE  LOADS 


It  has  been  shown  in  previous  sections  that  it 
frequently  becomes  desirable  to  replace  the  body  forces 
due  to  earth  attraction  by  a  series  of  discrete  surface 
loads.  This  section  is  devoted  to  an  analysis  of  the 
errors  resulting  from  the  approximation,  and  to  methods 
that  may  be  employed  to  compensate  for,  or  at  least  to 
minimize,  the  error. 

Consider  a  structure  loaded  with  n  discrete 

surface  forces  L.  such  that  their  sum  P„.  satisfies  the 
J  m 

relationship 


* 

P  ■  /L  L  i  *  V  w  (? 

"  i”  *  *  p  %  V  (2.7.1) 

derived  in  Section  (2.6.7).  Divide  the  total  volume  Vm 


of  the  structure  into  n  subvolumes  Vj,  each  proportional 

o  ' 

Vw,  L> 


to  its  load  Lj,  so  that 


V-  «  L  TS* 


(2.7.2) 

The  subdivision  has  to  be  carried  out  in  such  a  manner  that 

the  influence  volumes  are  as  geometrically  similar  as 

possible.  Assume  that  a  typical  Influence  volume,  as 

shown  in  plan  in  Fig.  2.7.1,  is  small  enough  so  that  its 

thickness  t  can  be  considered  constant.  Let  load  L,  be 

'  J 

applied  as  a  uniform  stress  p  over  area  a^  such  that  the 
centers  of  gravity  of  Vj  and  a^  coincide  in  plan.  The 
problem  now  becomes  the  determination  of  the  stresses 
caused  in  element  v^  by  the  replacement  of  the  uniformly 
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distributed  body  force,  which  should  exist  according  to  the 
requirements  of  the  exact  theory  of  similitude,  by  the 


surface  force  L  ^ .  In  the  model,  these  stresses  will  be 
added  vectorially  to  the  stresses  caused  by  the  overall 


structural  action. 


Consider  a  section  ZZ  through  the  center  of 
gravity  of  element  Vy  taken  so  that  the  length  +  A2) 

is  a  minimum  for  the  figure.  If  shape  ABCD  i3  not  too 
irregular,  the  direction  of  ZZ  will  be  a  good  approximation 
to  the  direction  of  maximum  stiffness  of  the  element.  To 


obtain  an  upper  limit  to  the  error  stress  field,  it  will 
thus  be  safe  to  consider  a  slice  of  unit  thickness,  as  shown 
in  Pig.  2.7.2,  where  Z  is  the  smaller  of  lengths  and  tg. 
The  analysis  will  be  based  on  the  following  assumptions  and 
conditions: 


i)  Since  an  upper  limit  to  the  error  stress  field 
is  sought,  a  condition  of  plane  stress  over  EFGH  may  be 
assumed. 

ii)  The  model  material  is  homogeneous  and  linearly 
elastic  In  the  range  considered.  This  implies  either  that 
the  entire  stress-strain  diagram  is  linear,  or  that  the 
additional  stresses  caused  by  the  approximation  of  the  body 
forces  are  so  small  that  within  that  particular  range  the 
stress-strain  diagram  is  linear. 

iii)  Deformations  are  small.  This  does  not  imply  that 
the  model  as  a  whole  has  to  undergo  small  deflections,  but 
rather  that  the  relative  deflections  of  points  within 
EPGH  have  to  be  small. 


iv)  Except  for  stress  p,  faces  EH  and  PG  are  free  of 
external  loads. 


v)  Element  is  restrained  along  EP  and  GH  in  such 

a  manner  that  E  and  H  do  not  displace  and 

3u.  _  'Sv  _  A 

* 0 

along  EP  and  GH.  This  condition  satisfies  symmetry  require¬ 


ments  if  all  influence  volumes  Vj  are  identical  and  carry 
identical  loads.  P  and  G  do  displace  due  to  the  elastic 


deformation  of  lines  EF  and  GH. 
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FIGURE  2.7.1 

PLAN  OF  TYPICAL  INFLUENCE  VOLUME 


FIGURE  2.7.2 

SECTION  THROUGH  TYPICAL  INFLUENCE  VOLUME  . 
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vi)  The  nature  of  the  supports  of  the  structure  as  a 
whole  does  not  affect  the  error  stress  field  significantly 
since 

a)  the  external  loads  are  balanced  by  the  body 
forces  so  that  no  edge  reactions  ere  required. 

b)  condition  (v)  eliminates  the  effect  of  re¬ 
straints  against  edge  rotation. 

c)  small  errors  are  csused  if  edge  EP  is  not 
rigidly  restrained  at  the  support. 

Since  the  body  forces  Y  are  balanced  by  the  applied  load, 
equilibrium  requires  that  for  unit  thickness  of  element  EFGH 

2teY  =  -  2bp 


Y  = 


(2.7.3) 


Using  Timoshenko's  Notation  (12}  ,  the  equilibrium  equa¬ 
tions  are 


It  +  IVs  * 0 

*  ° 

and  the  strains  will  be  compatible  if 


The  boundary  conditions  are 


(2.7.1+) 

(2.7.5) 


v0: 

v  s*° 

-l  k*  fc  +  t 

(2.7.6) 

V*  '• 

C  fc  4-fc 

0 

(2.7.7) 

IP  \*\>b9 

P  IF  \x\  *  b 

(2.7.8) 

*■*»**  Ik  *° 

(2.7.9) 

The  solution  is  simplified  considerably  by  the  introduction 


of  the  Airy  Stress  Function^: 

M  *  •»  ♦ 

*  * 
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(2.7.10) 


The  equilibrium  equations  are  identically  satisfied  and 
the  compatibility  condition  (2.7.5)  reduces  to 

*  0 

A  solution  to  this  partial  differential  equation  is 

*  to*  ox  [c,  co *h  a^  +  ct  o'mh  cosh  oi^  -r  c4y*mh  a<^ 


Using  (2.7.10): 

<£■  co*  ax[c,a*co*h  4-C^oSinh  ay  Cjo(2»i«V»  ay  ojj  cosh  c4ct(*GM>»  ay  a^  smVi 
fry  a?  CO*  ax  [  c,  cosh  a^  +  ct  *mVt  a^  4  Caty  co*h  a^  a^j  -  ^Y 

t^«  -  a  o'ni  ax[c,a  sinhoy  qa  to*K  cuj  4  cs(co*h  ay  a^  Mnh  C4(i**V\  ay  a^  cosh  ayjj 

bu.  .  bv 


Using  (2.7.9):  T,*  •  ■== 


Using  (2.7.6): 


Using  (2.7.7) 


n'^+^a  1 

at  s  yiir 

c^  =  0 

02^  +  c-,  =  0 


C  „  at  *iwh  at 
*  *i*V»  at  t  at  cosh 


ON  x=±t 


ZiV  <*»ci 


To  satisfy  (2.7.8),  a  Fourier  Series  for  O'  is  required. 

v 

Using  the  last  three  results,  and  replacing  for  Y  from 
equation  (2.7.5) 

V&  V  +f  MjnKMk 


a^  4  ^  cosn 


<*«  4 


0  v  n«i  w  -  - 

•UT  FOR  yt,  <Ty  +  £,4*  *F*  *** 

.  c  m _ *  p  *in  at _ _ 

*  a*t  £  v»*h  at  +  at  co*h*t  +  V»d»  afcj 
All  the  coefficients  have  thus  been  determined,  and  sub¬ 
stitution  yields 

■-J  *-  ^Vtv«h  a^+  a^co*Ka^  +  (  Zcosh  ayt  ayiwha^t^  co*  ax 

(2.7.11) 
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*2- 

p  tl  l 


slnh  ay  +  ay  cosh  ay  +  dn  ay  sinh  ay 


n=l 


)*" 


cos  ax 


¥■£{ 


(2.7.12) 


ay  slnh  ay  +  dn  (slnh  ay  +  ay  cosh  ay) V  tn  sin  ax 

(2.7.13 

_ 2  sin  ab _ 

where  *n  -  at  +  at  (cosh  at  +dn  slnh  at)]  (2. 7. II4.) 


with  dn  = 


at  sinh  at 


sinh  at  +  at  cosh  at 


and  a  =  jpl  (2.7.15) 


The  results  are  In  non-dimensional  form,  depending  on  the 
ratios  t/l  and  b/t.  It  remains  to  prove  that  boundary 
condition  2.7.9  is  satisfied. 

The  displacement  in  the  y-directlon  is 

V  =  i  [(fry  -^  *x)  d^.  +f(x) 

•••&  ■  l  ->*-£)  **  ♦f'U) 


Since 


^rx  „  ^  ,  nirx 

and  contain  sin  — * — 


rx  o  y 

on  x  =  t  t  s°  th8t  Vj, 


,  the  series  vanishes 


df*  «  y  dyl  +  g’(x)=  g'(x) 

*  x  L  J  tlE  J 


=  0 


t  d  v 

Choosing  g  (x)  =  0  ensures  that  =  0  for  all  y  on 

x  =  —  C.  Since  T-  =  0,  this  also  ensures  that  ~--u- 

y  dy 

for  all  values  of  y  on  x  =  —  t.  The  condition  that  u  =  0 
on  1  t  cannot  be  met  since 


u  =  H 


(«■ 


.i 


jdbP 

•*„)  dx  +  f(y)  * - xy  +  f(y) 

y  tlE 


Choosing  f(y)  =  0  shows  that  u  is  caused  by  Poisson 
effects  only.  This  deviation  from  the  original  conditions 
will  be  neglected. 
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Examples:  Equations  (2.7.11)  to  (2.7.13)  will  now  be  used 
to  derive  the  error  stress  distribution  for 

i)  t  =  t  (shear  type  structure) 

ii)  t  *  0.1  £  (bending  type  structure) 

In  each  case,  three  alternate  widths  of  load  distribution 
are  investigated: 

b  =  0.1  l,  b  =  0.25  t  and  b  =  0.50  l 

The  case  b  =  V  is  trivial  since  the  series  terms  in  equations 
(2.7.11)  to  (2.7.13)  vanish. 

2.7.1  Error  Stress  Field  if  Height  t  Equals  Halfspan  R, 
This  case  is  typical  of  massive  structures  where  earth  attrac¬ 
tion  forces  may  be  of  great  Importance.  It  is  considered 
of  interest  to  determine  both  the  error  in  strain  readings 
on  the  surface,  and  the  distribution  of  major  principal 
stress  and  maximum  shear  stress  throughout  the  mass.  In 
order  to  facilitate  similar  calculations  for  other  ratios 
of  t:l,  typical  examples  of  the  computations  are  given  below. 
It  will  be  noted  that  the  first  four  terms  of  the  series 
solution  are  evaluated  with  slide  rule  accuracy.  Conver¬ 
gence  is  good  except  at  *  =  0,  b  =  0.1  Z,  y  =  t  where  extra 
terms  were  added. 

Table  2.7.1  -  The  information  contained  in  Table  2.7.1  and 
Table  2.7.2  is  valid  for  all  values  of  y. 

Table  2.7.2  -  The  values  of  O’.,,  O',,  andT,.,,  have  been  calcu- 

a  y  *y 

lated  for  y  =  0,  0.4t>  0.6t,  0.8t  and  t,  for  each  of  the 
three  values  of  b.  The  major  and  minor  principal  stresses 
0*^  and  CTg,  the  maximum  shear  stress  t***  and  the  angle  cc 
between  o*n  and  were  then  found  from 
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Table  2.7.2  Typical  Calculations  for  the  Error  Stress  Field 


.759  8.51x10"  2.6I4JC  10"-5  6.00x10“-^  8.51x10" 

.865  9.1i5xl0"5  5.52X10"5  9.45xlO"5  0 

.905  1.86xlO'6  1.50X10"6  -  1.31X10"6  -1.86x10"' 

.930  4.50xlO'8  4.27X10'8  0  0 


Only  the  computations  for  i^=  0.6t,  b  *  0.25  t  will  be  shown 
here.  Values  of  maximum  normal  and  shear  stresses  for 
other  values  of  y  and  b:£  can  be  read  off  Pigs.  2.7.1  to 
2.7.3  and  Tables  In  Appendix  .B.  Since 
•frinh  0-6  n*  *  cosh  0*6  wtr 


equations 


(2.7.11  to(2.7.13)  reduce  to 
s  —  21  ^  I  +  0*6ntr  +  d*  (.2.+  0*6mr)|  ft 
•  0*tS  +  £  |  l  +■  0*6  trn  + 

«  ^  w1r  +  dy,  (l  + 0*6  vnr)j 


M  S\wV>  0-6  yiV  CO*  AX 
sin  Vi  o-fcnv  cosax 
&n  Stv\Vt  0*6  via*  sin  ox 


« 

0*6viir 

— 

e*ti«iVt  0 >6nt 

Coeff 
or  SI 

o-k/p 

Icient  of 
n  ax  for 

«Vp 

Cos  ax 

Wp 

1.88 

3.77 

5.67 

7.51]- 

1.97  x  10”2 
2.05  x  10'3 
-1.90  x  10*^ 
0.0 

-0.1150 

-0.0200 

+0.0026 

0.0 

0.0850 

0.0165 

-0.0022 

0.0 

0.0800 

0.0161 

-0.0022 

0.0 

Using  the  table  for  4i**^^and  eo*2^,  we  find  the 
following  results. 
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Table  2,7.3  Typical  calculations  for  the  Error  Stres 


Coefficient  *  cos  ax 


+0.0439 

+0.0810 

+0.1060 

+0.1150 


-0.0200 

-0.0141 

0 

+0.0141 

+0.0200 

+0.0141 


-0.0141 

-0.0200 


+  0.0026 
+0.0010 
-0.0018 
-0.0024 
0 

+0.0024 

+0.0018 

-0.0010 

-0.0026 


Coefficient  x  cos  ax 


0.0165 

0.0116 

0 

-0.0116 

-0.0165 

■0.0116 

0 

0.0116 

0.0165 


-0.0022 

-0.0008 

+0.0016 

+0.0020 

0 

-0.0020 
-0.0016 
+  0.008 
+0.0022 


-0.1324 
-0.1191 
-0.0828 
-0.0322 
+0.0200 
+0.06 04 
+0.0828 
+0.0909 
+0.0924 


0.2493 

0.2391 

0.2116 

0.1728 

0.1335 

0.1040 

0.0884 

0.0841 

0.0837 


Table  2.7.U  Typical  Calculations  for  Error  Stress  Field 
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The  results  of  the  rest  of  the  computations  are  summar¬ 
ized  In  Appendix  B,  Tables  B2  to  B6.  The  major  principal  and  maximum 
shear  stresses  given  In  Tables  .B.2  to  .B.6  are  presented 
graphically  In  Figures  2.7.3,  2.7.1;  and  2.7.5.  The  variation 
of  horizontal  stress  O'*  on  the  lower  face  is  shown  In  Figure 
2.7.6,  but  the  variation  of  ffx  on  the  upper  face  (y  =  0)  Ip 
so  small  that  no  graph  has  been  drawn. 

Conclusions 

I)  The  additional  stresses  Introduced  into  the  model 
by  the  replacement  of  body  forces  by  a  surface  load  are  of 
the  same  order  of  magnitude  as  the  applied  stress  p,  rather 
than  the  total  applied  load.  If  a  1:21;  mortar  model  of  a 

21;"  thick  concrete  shell  Is  loaded  on  a  2"  x  2"  grid  (...  t- 1) , 
and  fa  =  fy  =  1,  fg  =  0.6  then  from  Eqn.  3.6.5  f g  =  24 

and  from  Eqn.  2.7.1  the  load  applied  to  each  Influence  area 
Is 

LJ  =  vj  wp  (VV 

nr  *o.k-s.wib. 

For  b  =  0,10  l:  p  =  135  psl 

b  =  0.25  Is  p  =  22  psl 

b  =  0.50  l:  p  =  5  psl 

A  great  reduction  in  the  additional  (undesirable)  stresses 
can  thus  be  achieved  by  distributing  the  applied  load  over 
a  large  percentage  of  its  Influence  area. 

II)  Strain  gage  readings  on  the  upper  face  of  the 

element  In  Fig.  2.7.2  will  not  be  affected  appreciably  by 
the  replacement  of  body  forces  bty  surface  forces  (  <.  1.3# 

of  p).  Strain  gage  readings  on  the  lower  face  are  affected 
considerably  (  O'*  <  115#  of  p).  Strain  gages  should 

thus  always  be  mounted  on  the  face  opposite  to  that  on 
which  the  load  is  applied. 
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STRESSES  IN  A  TYPICAL  INFLUENCE  VOLUME  POE  TO  THE  REPLACEMENT  OF  BOOV 
FORCES  BY  A  SURFACE  LOAD  IF  THE  WEIGHT  t  EQUALS  THE  HALFSPAN  t. 


MAXIMUM  SHEAR  STRESS  I  MAXIMUM  NORMAL  STRESS 
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EAR  STRESS  I  MAXIMUM  NORMAL  STRESS. 
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STRESSES  IN  A  TYPICAL  INFLUENCE  VOLUME  DUE  TO  TME  REPLACMENT  OF  SOPV 
FORCES  BEY  A  SURFACE  LOAD  IF  THE  HEIGHT  t  EOUALS  THE  HALFSPAN  t. 


THICKNESS  t  m  HALFSPAN 
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FIGURE  2.7. 6 

STRESSES  »N  a  TYPICAL  INFLUENCE  VOLUME  ALONG  THE  FACE  ON  NSfWlCH  the 
SURFACE  LOAD  IS  APPLIED . 
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TTPlCAL  INFLUENCE  VOLVNlE  ALONG  THE  FACE  ON  WHICH  THE  SURFACE 


TMICKNE%ft  t  3  0*1  *  HALPSPAN 


iii)  The  replacement  of  body  forces  by  surface 
loads  may  considerably  affect  the  results  of  ultimate  strength 
studies  of  mass  concrete  or  lightly  reinforced  structures. 

The  additional  tensile  stresses  may  lead  to  premature  crack 
formation  in  the  model,  or  they  may  cause  the  model  to 
crack  in  places  where  the  prototype  will  not  crack.  Figures 
2.7.3  to  2.7.5  do,  however,  show  that  the  influence  of  the 
surface  load  attenuates  rapidly  inside  the  mass,  so  that  a 
slight  strengthening  of  the  material  in  the  vicinity  of  the 
load  point  will  eliminate  the  harmful  effects.  Since  the 
stress  distribution  is  known,  the  amount  of  additional 
reinforcment  required  is  readily  calculated. 

iv)  The  effect  of  the  body  force  simulation  on 
the  buckling  behavior  is  not  studied  since  buckling  is  seldom 
a  problem  in  massive  structures. 

2.7.2  Error  Stress  Field  if  Height  t  =  0.1 

times  Half-span  L.  The  loadpoints  are  spaced 
at  twenty  times  the  thickness  of  the  element.  This  is 
typical  of  large  span  bending-type  structures.  The  horizon¬ 
tal  stresses  O' x  on  the  upper  and  lower  faces  of  the  element 
have  been  calculated  from  the  first  four  terms  of  the  series 
solution,  except  at  x  =  0,  b  =  0.1  &  where  further  terms 
were  added.  The  results  are  given  in  Table  2.7.5  and  are 
shown  in  Figures  2.7.7  and  2.7.8 

Conclusions 

i)  As  for  the  case  t  =  t,  the  additional  stresses 
introduced  into  the  model  by  the  replacement  of  body  forces 
by  a  surface  load  are  of  the  same  order  of  magnitude  as  the 
applied  stress  p,  rather  than  the  total  applied  load.  It 
is  thus  again  of  advantage  to  distribute  the  applied  load 
over  a  large  percentage  of  its  influence  area. 

ii)  Strain  gage  readings  on  both  the  upper  and  the 
lower  face  of  the  element  are  affected  by  the  replacement 
of  body  forces  by  surface  forces  (  (*■*)  max  =  30#  and  100# 
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of  the  applied  stress  p  respectively).  It  is,  however, 
possible  to  place  the  strain  gage  so  that  it  does  not 
register  the  error  stress  field.  The  exact  position 
in  this  case  depends  on  the  ratio  b:  i,  and  varies 
bet'-reen  the  quarter  and  half  points  of  Z  (see  Pig.  2.7.8). 

iii)  If  the  load  is  transferred  to  the  model  with 
a  pad,  the  stiffness  of  the  model  is  altered  locally.  This 
effect  has  not  been  taken  into  account  in  the  preceding 
analysis  since,  in  the  case  of  massive  structures,  the 
relative  stiffness  of  the  distribution  pad  can  be  made 
small. 


iv)  The  effect  of  the  body  force  simulation  on 
the  buckling  behavior  has  not  been  studied.  For  normal 
(not  massive)  structures  where  this  problem  would  be  of 
interest, it  could  be  investigated  by  determining  the  strain 
energy  locked  up  in  the  error  stress  field 


,.E.  .  /  f  (  *  A|5 


r*y 

2G 


dxd, 


This  could  then  be  compared  to  the  total  strain  energy 
required  to  buckle  the  structure.  Local  snapthrough  at 
the  loadpoint  would  also  have  to  be  investigated. 
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Table  2.7.5 

Ratio  of  horizontal  to  applied  stress  for  t  =  0.1  Jb 


X 

I 

b=0.10  l 

y  =  0 

b=0. 25  4 

0 

-0.227 

-0.285 

-0.236 

-1.000 

-0.910 

-0.660 

0.125 

-0.1^8 

-0.258 

-0. 247 

-0.535 

-0.844 

-0.752 

0.250 

-0.008 
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CHAPTER  3 


THE  TECHNICAL  ASPECTS  OF  THE  MODELING  OF  STEEL  STRUCTURES 

3.1  SIMILITUDE  REQUIREMENTS 

We  summarize  here  for  convenience  the  similitude 
requirements  for  the  static  and  dynamic  response  of  struc¬ 
tures  loaded  beyond  their  elastic  limit.  We  shall  assume 
for  simplicity  that  the  correlation  functions,  f,  do  not 
change  with  position  or  time  as  they  possibly  would  in  the 
general  case. 

3.1.1  Static  Case.  In  the  static  case  assuming  that 
we  know  our  prototype  material,  in  this  case  structural  steel, 
we  must  examine  the  possibility  of  finding  a  model  material 
with  similar  elasto-plastic  behavior.  The  correlation 
functions  obtained  from  a  dimensional  analysis  together 
with  the  compatibility  of  the  material  properties  of  the 
model  to  prototype  will  determine  the  design  of  the  model 
and  the  prediction  equations  from  which  the  prototype 
behavior  will  be  extrapolated. 

3. 1.1.1  Basic  Assumptions.  In  the  modeling 
technique  which  is  proposed  here  the  modeling  materials  are 
chosen  on  a  rational  basis  and  then  investigated  to  see 
how  accurately  they  can  predict  the  prototype  behavior  as 
compared  to  actual  tests.  The  scope  of  this  modeling  technique 
will  be  limited  to  structures  in  which  the  following  assump¬ 
tions  hold: 

(1)  The  gravity  stresses  due  to  the  dead 
weight  of  the  structure  are  not  important  in  the  study  and 
if  they  do  become  critical,  external  means  will  be  employed 
to  satisfy  similitude. 

(ii)  The  structural  members  are  essentially 
uniaxial  so  that  Poisson's  ratio  need  not  be  Identical  in 
model  and  prototype. 
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(iii)  Boundary  conditions  and  loading  will  be 
duplicated  adequately  and  are  not  the  source  of  serious 
model  errors. 

(iv)  In  the  static  case  we  assume  that  there 

are  no  important  differences  between  the  temperature  of  the 

model  and  the  prototype,  ft  *  1;  and  that  the  specific  heat 

of  the  materials  does  not  enter,  thus  f  =  1;  the  coefficient 

c 

of  linear  expansion  will  not  come  into  play  so  that  i*=  1; 
and  also,  f =  1  since  thermal  conduction  is  no  problem. 

3. 1.1. 2  Correlation  Functions.  In  the 
static  case  time  is  no  longer  a  variable  thus  we  see  that 
the  only  relationship  which  must  be  satisfied  since  acceler¬ 
ations  are  also  arbitrary  is  the  stress  ratio.  If  the 
model  material  is  known  then  the  stress  ratio  is  determined 
from  the  stress-strain  relationship  of  the  model  and  proto¬ 
type  materials.  The  length  scale  is  chosen  arbitrarily. 

If  body  forces  come  into  play,  that  is  gravity  stresses 
must  be  modeled,  then  the  relation  which  must  be  satisfied 
is  j 

*«•  *  Vr 

which  determines  the  length  scale  if  a  model  material  is 
chosen. 

3. 1.1. 3  Stress-strain  Relationship.  By 
neglecting  the  effect  of  time  and  temperature,  we  need  only 
investigate  the  stress-strain  behavior  of  the  materials. 

We  shall  consider  the  behavior  of  homologous  points  in 
model  and  prototype  which  in  most  cases  will  be  the  most 
highly  stressed  points  which  are  of  interest. 

We  assume  further  that  the  stress- 
strain  relations  for  a  tension  specimen  will  be  representa¬ 
tive  of  all  stress-strain  relations  for  compression  and  shear 
connecting  all  shear  strain  components  to  the  respective 
stress  components  in  the  general  case.  Knowing  the  stress- 
strain  relation  for  the  model  material  chosen  one  can  then 
determine  the  correlation  functions  which  make  it  compatible 
to  the  prototype  material. 
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3. 1.1. 4  The  Use  of  Phosphor  Bronze  as  a 
Modeling  Material  for  Structural  Steel,  The  first  example 
considers  the  stress-strain  relations  shown  in  fig,  3,1,1. 
This  figure  shows  the  idealized  stress-strain  relationship 
of  structural  steel  and  phosphor  bronze.  The  compatibility 
requirements  for  the  different  ranges  of  strain  follow: 


When 

0  ^  €.  ^6y 


we  have. 

CT n  =  E  £ 

m 

(3.1.2) 

0~p  =  £p  €p 

(3.1.3) 

Since  €m  =  €p  and  <5Tm  =  fff^p  3  it 

follows  from 

Equation  (2.3.10)  that 

^§m  _  f^(7p 
£>€  m  ^c^p  • 

(3.1.4) 

Em  =  f q-  Ep 

(3.1.5) 

The  materials  are  compatible. 

when 

CT  m  =  CTym 

(3.1.6) 

CT  p  =  ^yp 

(3.1.7) 

It  follows  from  equations  (3.1.2),  (3.1.3), 
that 

and  (3.1.5) 

O'  ym  Gy 

~  _  _  £ 

0  7V  &P  €  y  °* 

• 
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FIG.  3.U  IDEALIZED  STRESS-STRAIN  RELATIONS  FOR 
STRUCTURAL  STEEL  AND  PHOSPHOR  BRONZE. 
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It  follows  from  equations  (3.1.6)  and  (3.1.7)  that 


^  m  _  ^ym 

(T  P  0“  yp 


(3.1.8) 


The  materials  are  thus  compatible. 


When 


£  stm  ^  G 


0"  m  -  (y  jm 
(Tp 


=  E *m  (  £  m  -  £stm  )  (3.1.9) 

CT  yp  (3.1.10) 


Since  =  £p,  substitution  of  equation  (3.1.8) 
into  equation  (3.1.9)  yields 

v  (  Cp  -°Vp>  -  <  <eP  -  0 


(3.1.11) 


Since  £ p  £  stm,  this  requirement  cannot  be  fulfilled 
i.e.,  the  two  materials  are  incompatible  at  this  range. 
Therefore  we  arrive  at  the  conclusion  that  above  £  ^£stm 
the  model  stress-strain  curve  should  have  zero  slope. 

This  means  that  both  model  and  prototype  material  should 
strain  harden  at  the  same  strain.  If  the  plastic  deformations 
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are  not  excessive  it  nay  still  be  possible  to  use  this 
modeling  material  up  to  strains  in  the  range  of  £  ^ 
however  it  is  unlikely  that  the  full  plasticity  of  a  member 
could  be  developed  unless  plastic  strains  of  the  order  of 
magnitude  of  £stp  are  available.  We  shall  see  from  actual 
tests  that  the  plastic  strain  range  is  enough  to  develop 
full  plasticity  of  the  member. 

3. 1.1.5  The  Use  of  Ethyl  Cellulose  as  a 
Model  Material  for  Structural  Steel.  However,  if  a  material 
such  as  ethyl  cellulose  is  used  to  model  steel  the  stress- 
strain  relations  will  be  as  shown  in  fig.  3.1.2.  Looking 
at  fig.  3.1.2  and  using  the  results  of  section  3.1. 1.4 
we  can  deduce  that  for 

0  ^  e  > 


For 


^  yp  ^  ^  € atp 

there  is  no  compatibility  unless  Em  =  0.  In  other  words 

for  strains  beyond  but  less  £  .  the  flatter  the  slope 

J P  3 

°f  Em  (the  weaker  the  material  used  for  the  model)  the 

better  is  the  compatibility  of  prototype  and  model  materials 

However,  in  all  such  cases  there  will  be  discrepencies  in 

the  model,  i.e.,  the  model  will  be  "distorted”  and  the 

degree  of  distortion  will  be  a  direct  function  E  .  In  such 

m 

cases  the  model  results  cannot  be  expected  to  be  of  a  quanti 
tative  nature  if  the  model  is  strained  beyond  the  yield 
strain  of  the  prototype. 

When  ^stp  «  £  S  fcym  there  wln  b® 
a  change  in  the  value  of  fff.  It  is  now  given  by 
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STRESS 


STRAIN,  € 


FIG.  3.1.2  IDEALIZED  STRESS- STRAIN  CURVES  FOR  STRUCTURAL 
STEEL  AND  ETHYL  CELLULOSE. 


9k 


(3.1.13) 


If6ym  e  <6stm.  th®n  th® 
the  strain  hardening  modulus  of  the  steel  stress-strain 
curve  the  better  the  correlation.  It  Is  unlikely  that  the 
prototype  will  be  able  to  strain  more  than  the  early  part 
of  the  strain  hardening  range  at  the  critical  sections 
where  plastic  deformations  are  taking  place  before  the 
structure  falls  by  local  or  lateral  buckling.  Thus  €stin 
Is  probably  the  upper  limit  of  strain  to  be  considered. 

It  should  be  apparent  from  the  dis¬ 
cussion  and  fig.  3.1.2  that  the  relative  distortion  in  the 
model  will  depend  on  the  positions  of  6^  with  respect  to 
€yp.  The  closer  these  two  values  are  the  better  the  corre¬ 
lation.  Also,  the  closer  the  values  of  €  .  and  £  the 

stp  stm 

better  the  correlation. 

3.1.2  The  Dynamic  Case.  Our  main  concern  in  this 
study  is  the  dynamic  response  of  structures  on  which  the 
dynamic  loads  are  a  result  of  an  air  blast. 

The  similitude  requirements  covering  this  phenom¬ 
enon  are  derived  In  Chapter  2.  Although  similitude  does  not 
necessarily  impose  a  stress  ratio  of  unity  In  the  general 
case,  the  phenomenon  of  the  air  blast  In  order  to  be  modeled 
with  confidence  Imposes  a  stress  ratio  of  one,  l.e.,  the 
model  and  the  prototype  must  be  tested  at  the  same  pressure. 
This  means  that  In  most  cases  the  Same  material  must  be 
used  in  moa^l  and  prototype  since  the  stress-strain  relation¬ 
ship  for  the  two  materials  must  be  identical.  However, 
there  may  be  situations  where  this  restriction  is  not 
necessary.  If  the  structure  Is  made  of  open  or  frangible 
type  construction  such  as  bridges,  towers,  trusses  and 
frames  the  main  dynamic  loading  which  results  Is  of  the 
drag  type.  In  such  cases  the  complications  resulting  from 
the  reflected  overpressures  as  functions  of  the  angle  of 
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Incidence  and  the  peak  overpressure  do  not  arise.  It  will 
be  possible  then  to  test  the  model  at  a  smaller  pressure 
than  the  prototype.  This  means  that  a  weaker  material 
than  the  prototype  can  be  used  for  the  model.  It  is  for 
such  cases  that  the  modeling  technique  using  phosphor  bronze 
was  developed  in  the  study  of  steel  structures. 


From  the  basic  relation, 


f 


s 


f  f  f 
gar 


(3.1.14) 


in  the  case  of  phosphor  bronze  we  have: 


=  15  x  1  x  -  =  7.5 
? 

This  Implies  that  the  inertial  accelerations 
in  the  model  will  be  7.5  times  those  in  the  prototype. 

The  result  will  be  a  faster  rate  of  strain  of  the  model 
material. 


3.2  CHOICE  OF  MATERIALS 

The  modeling  of  steel  structures  is  a  compara¬ 
tively  new  phase  in  model  analysis.  Whereas  in  the  past 
many  elastic  model  studies  of  steel  structures  have  been 
made  these  were  done  mainly  in  plexiglass  or  some  other 
easily  shaped  plastic.  Small  scale  model  studies  in  metal 
are  therefore  relatively  scarce  in  the  literature.  The 
concept  of  loading  the  model  beyond  the  elastic  range  to 
predict  the  ultimate  behavior  of  the  prototype  and  the 
other  restrictions  imposed  by  similitude  requirements  in 
the  study  of  the  dynamic  response  of  structures  cannot  be 
attempted  without  a  clear  understanding  of  the  mechanical 
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properties  of  the  materials  involved.  In  this  section  the 
materials  of  both  the  prototype  and  the  proposed  models 
will  be  investigated. 

3.2.1  Structural  Steel.  The  mechanical  proper¬ 
ties  of  a  material  are  greatly  influenced  by  its  micro¬ 
structure  a  knowledge  of  which  is  very  essential  in  the 
choice  of  a  model  material  which  will  be  compatible  with 
the  prototype.  However,  we  shall  leave  this  realm  of 
material  behavior  to  the  specialists  in  this  newly  expand¬ 
ing  science  and  as  structural  engineers  concentrate  on  the 
mechanical  properties  which  are  of  the  most  immediate  use 
to  us. 

3. 2. 1.1  Definition  of  Structural  Steel. 

Steel  is  basically  an  alloy  of  iron  and  less  than  2$6  of 
carbon  by  weight.  Many  other  metals  can  be  alloyed  with 
it  to  give  steels  of  various  characteristics.  We  shall 
be  concerned  here  with  structural  steels  which  owe  their 
strength  and  other  properties  chiefly  to  carbon.  According 
to  the  AISC  specifications  of  November  1961,  Section  2.2 
defines  "structural  steel"  as  follows: 

"Structural  steel  shall  conform  to  one  of  the 
following  specifications,  latest  edition: 

Steel  for  Bridges  and  Buildings,  ASTM  A 7 
Structural  Steel  for  Welding,  ASTM  A373 
Structural  Steel,  ASTM  A36" 

Table  2.1  summarizes  the  properties 
of  the  structural  carbon  steels  which  we  shall  attempt  to 
model. 

Of  late  a  new  series  of  steels  which 
have  made  their  appearance  on  the  market  seem  to  show  great 
promise  in  imaginative  applications  in  structures 
resulting  in  more  economical  designs.  These  are  the  high 

strength  and  high  strength  low  alloy  steels  described  in 
table  2.2.  The  methods  which  we  shall  propose  of  modeling 
steel  structures  can  be  expanded  very  easily  to  fit  any  of 
these  steels.  The  stress-strain  curves  for  the  steels  of 
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LOW-ALLOY  STEELS  (U.  S.  Steel  Co. 
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To  one  Inch  Inclusive 

Over  one  inch,  varies  for  different  specifications 


tables  3.1  and  3.2  are  shown  in  figure  3.2.1.  These  are 
taken  from  U.S.  Steel  Co.  figures  and  indicate  the  minimum 
guaranteed  yield  points.  Note  that  the  plastic  range  of 
steel  decreases  as  the  yield  strength  increases.  In  compar¬ 
ing  tables  3.1  and  3.2  note  that  because  the  chemical  com¬ 
position  of  structural  steels  varies  with  thickness  the 
mechanical  properties  such  as  yield  and  tensile  strength 
are  about  the  same  for  all  thicknesses.  However,  in  high 
strength  steels  the  chemical  composition  is  the  same  and 
thus  the  strength  properties  are  a  function  of  the  thick¬ 
ness  as  shown  in  table  3.2. 

3. 2. 1.2  "Static"  Properties  of  Steel. 

The  most  important  mechanical  properties  of  structural  steel 
can  be  inferred  from  a  close  study  of  the  stress-strain 
curve  which  results  when  a  standard  specimen  is  pulled  in 
tension  at  a  specified  rate  of  strain.  In  fig.  3.2.2  we 
see  such  a  curve. 

(a)  shows  the  whole  stress-strain  behavior. 

(b)  a  magnification  of  its  early  portion  will  concern 
us  most  and 

(c)  the  idealization  of  the  early  portion  which  is 
often  made  for  analysis. 

The  slope  of  the  elastic  portion 
of  the  curve  is  known  as  Young's  modulus  and  is  often 
approximated  at  30  x  10^  psi.  The  most  important  pecul¬ 
iarity  of  this  stress-strain  diagram  is  the  distinct  yield 
point  which  is  almost  unique  to  this  material.  The  upper 
yield  point  £yu  is  unstable  and  the  large  plastic  deforma¬ 
tion  which  takes  place  is  at  a  constant  value,  6y^>  or  the 
lower  yield  point.  After  a  strain  of  about  1.6$  the  material 
begins  to  strain  harden  but  does  not  fracture  except  after 
considerable  elongation. 

The  uniqueness  of  the  yield  point 
in  structural  steel  can  be  explained  on  the  basis  of  its 
microstructure.  It  seems  that  carbon  atoms  strain  the  iron 
latice  which  upon  the  application  of  an  external  stress 
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TENSILE  STRESS, 


FIG. 3.2.1  IDEALIZED  STRESS  -  STRAIN  CURVES 
FOR  ASTM  STEELS 
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FIG.  3.2.2  TYPICAL  STRESS-STRAIN  CHARACTERISTICS  OF 
STRUCTURAL  STEEL. 
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field  the  strains  are  partly  relieved  and  will  cause  plastic 
yielding  at  a  constant  value  of  stress.  This  phenomenon  is 
also  observed  in  other  multiphase  systems  such  as  phosphor 
bronze  but  not  to  is  pronounced  a  degree. 

The  three  basic  methods  of  increasing 
steel  strength  are  adding  alloys,  heat  treating,  and  cold 
working.  Of  these,  alloying  is  the  basic  method  of  control¬ 
ling  the  properties  of  structural  steels.  The  yield  point 
and  the  ultimate  tensile  strength  rise  with  an  increase  in 
the  carbon  content  of  steel  but  a  loss  in  ductility  and 
weldability  accompanies  these  strength  gains.  This  is  why 
carbon  content  in  structural  steels  is  usually  limited  to 
0.28$,  the  specified  maximum  for  A373,  A36  and  Aljlj.0.  A7 
steel  is  unique  among  ASTM-apeclfled  structural  steels  in 
having  no  specified  limit  on  carbon  content.  Manganese 
and  silicon,  the  most  important  alloying  elements  apart  from 
carbon,  also  Increase  yield  and  ultimate  tensile  strength 
with  a  smaller  relative  loss  in  ductility  and  weldability 
for  a  given  percentage  of  element  added.  Copper  and  nickel 
are  the  principal  alloying  elements  used  to  improve  the 
corrosion  resistance  of  high  strength  steels.  It  should 
be  pointed  out  that  all  the  high-strength  structural  steels 
resist  atmospheric  corrosion  better  than  carbon  steels. 

3. 2, 1.3  Dynamic  Properties  of  Steel. 
Experimental  evidence  Indicates  that  the  strength  behavior 
of  most  materials  depends  among  other  things  on  the  rate 
of  strain  during  testing.  In  steel  these  changes  are  shown 

i 

in  fig.  3.2.3  for  the  case  of  ASTM  A-7  structural  steel. 

The  effects  of  increasing  rate  of 
strain  can  be  suxanarlzed  as  follows,  looking  at  fig.  3.2.3: 

1.  The  yield  stress  Increases  to  some  dynamic  value,  . 

2.  The  yield  point  strain  €y  Increases. 

3.  E,  the  modulus  of  elasticity,  in  the  elastic  range 
remains  constant. 
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DYNAMIC  YIELD  STRESS 


STRAIN,  e  (%) 


FIG.  3.2.3  EFFECT  OF  RATE  OF  STRAIN  ON  STRESS-STRAIN 
CURVE  FOR  STRUCTURAL  STEEL  (REF.  30) 


TIME  TO  REACH  YIELD  POINT,  typ  (SEC) 


FIG.  3.2.5  EFFECT  OF  RATE  OF  STRAIN  ON  YIELD  STRESS 
(REF.  30) 

104 


ij..  The  strain  at  which  strain  hardening  begins, 
also  increases. 

5.  The  ultimate  strength  Increases  slightly. 

The  most  important  effect  which  will  influence  the  design 
of  steel  structures  to  resist  dynamic  loads,  is  the  increase 
in  the  yield  stress.  In  fig.  3.2.k  the  percentage  increase 
is  yield  stress  is  given  as  a  function  of  the  rate  of  strain 
for  two  steels  of  different  static  yield  stress.  It  is 
evident  from  the  figure  that  the  increase  in  dynamic  yield 
point  is  greater  for  steels  with  lower  static  yield  points 
as  is  shown  by  the  higher  slopes  of  curve  A.  Figure  3.2.5 
8 hows  the  dynamic  yield  stress  as  a  function  of  the  time 
required  to  reach  that  value  of  stress  (t  )  for  ASTM  A7 

v  r 

steel.  From  this  curve  values  of  design  yield  stress  could 
be  found  knowing  the  time  to  reach  yield  stress  in  a  particu¬ 
lar  structure. 

3.2.2  The  Copper  Base  Alloys  as  Model  Materials. 

Our  first  investigation  in  the  possibility  of  finding  an 
appropriate  model  material  which  could  be  used  to  model  steel 
structures  is  the  family  of  metals  known  as  the  copper 
base  alloys.  The  similitude  requirements  for  the  case  of 
blast  loads  acting  on  structures  indicate  that  if  the  sur¬ 
face  loads  must  be  kept  the  same,  which  may  be  the  case  in 
some  structures  such  as  shells  and  other  non-uniform  sur¬ 
faces,  then  the  same  material  must  be  used  for  model  and 
prototype  since  strain  must  be  equal  also.  However,  in 
many  structures  made  of  thin  exposed  members  or  frangible 
sidings  it  is  the  drag  forces  which  are  important  and 
therefore  in  testing  models  of  such  structures  it  is  possible 
to  use  lower  pressures  than  in  the  prototype.  This  makes  it 
possible  to  use  a  large  number  of  materials  which  have  less 
strength  than  steel  but  still  compatible  to  its  stress- 
strain  diagram.  It  is  for  this  class  of  problems  of  welded 
construction  in  steel  for  which  the  modeling  technique 
described  in  subsequent  sections  was  developed. 
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FIG.  3.2.4  INCREASE  OF  LOWER  YIELD  POINT  OF  STEEL  WITH  STRAIN  RATE.  ( REF.42 ) 


3. 2. 2.1  Properties.  Alloys  having  copper  as 
a  base  are  on  the  whole  the  most  readily  machinable  metals 
and  the  ease  with  which  they  can  be  milled  Into  various 
shapes  may  be  of  great  advantage  In  model  making.  In  addi¬ 
tion  the  ease  of  jointing  them  by  the  relatively  simple 
process  of  silver  soldering  adds  greatly  to  their  favor  In 
cases  where  welded  structures  are  to  be  Investigated  or 
where  a  large  number  of  parts  must  be  assembled.  An  exam¬ 
ination  of  the  stress-strain  properties  of  most  of  these 
alloys  Indicates  that  a  distinct  yield  point  as  In  the 
case  of  structural  steel  does  not  exist.  However,  in  two 
commerlcal  alloys  namely  phosphor  bronze  and  silver  nickel 
a  definite  yield  point  can  be  established  after  annealing. 

The  important  difference,  however, ^la  that  the  plastic  yield¬ 
ing  in  these  materials  is  not  as  pronounced  as  in  steel. 

3. 2. 2. 2.  Phosphor  Bronze.  Of  the  two  commer¬ 
cially  available  copper  base  alloys  that  could  be  investi¬ 
gated  as  possible  substitutes  to  steel  in  model  making, 
phosphor  bronze  is  chosen  on  the  basis  of  lower  cost  and 
better  machlnabllity.  The  "free  cutting"  phosphor  bronze 
which  la  Investigated  contains  1$  lead  (added  to  Improve 
machlnabllity),  l\.%  tin,  1$  zinc  and  68£  copper. 

To  obtain  a  definite  yield  point 
In  phosphor  bronze  a  cold  drawn  specimen  must  be  carefully 
annealed  to  control  the  crystal  size.  The  percent  cold 
work  or  the  amount  of  reduction  In  cross-sectional  area 
during  the  rolling  process  must  be  between  5 0 %  and  7 5%. 

In  the  case  of  $0%  cold  worked  samples  stress-strain  curves 
from  strips  Indicate  an  optimum  annealing  temperature  of 
675°C  for  one  half  hour.  The  effect  of  cold  work  on  the  yield 
stress  of  phosphor  bronze  Is  shown  in  fig.  3.2.6  and  the 
effect  of  annealing  temperature  In  fig.  3*2.7. 

These  curves  were  established  by 
annealing  l/l;  Inch  rods  of  two  hardnesses  (50#  and  75%) 
at  various  temperatures  for  one  half  hour.  The  composition 
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FIG.  3.2.6  EFFECT  OF  COLD  WORK  ON  THE  YIELD  STRESS  OF 
PHOSPHOR  BRONZE. 


FIG.  3.2.8  DIMENSIONS  OF  PHOSPHOR  BRONZE  STRIPS 
USED  AS  TENSILE  SPECIMENS, 
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of  the  phosphor  bronze  was  that  described  above.  The 
stresses  quoted  in  fig.  3.2.7  as  ’’yield  stresses"  are 
actually  the  stresses  at  which  the  load  was  observed  to  re¬ 
main  approximately  constant  during  the  test.  It  is  inter¬ 
esting  to  note  that  the  annealing  has  a  greater  influence 
in  lowering  the  yield  stress  in  samples  of  larger  cold  work. 
This  may  be  due  to  the  fact  that  in  this  case  the  crystals 
are  left  more  distorted  and  strained  by  the  rolling  proc¬ 
esses. 

3. 2. 2. 2.1  Static  Stress-Strain  Relations. 
The  stress-strain  properties  of  phosphor  bronze  are  not  readily 
available  in  the  literature.  It  becomes  necessary  then 
to  test  actual  strips  the  size  of  the  flanges  of  the  wide 
flange  model  beams  which  were  machined  for  the  various  tests. 
Fig.  3.2.8  shows  the  dimensions  of  the  test  strips  and 
the  precautions  taken  at  the  end  to  prevent  failure  at 
the  ends  by  the  jaw  grips  of  the  testing  machine.  The 
tests  are  made  on  the  constant  rate  of  strain  "Ingtrom" 
type  of  machine  of  the  Plastics  Research  Lab.  at  M.I.T. 

A  gage  length  of  2  inches  was  used  and  externally  controlled 
to  define  the  points  on  the  stress-strain  diagram.  This  is 
a  definite  improvement  over  other  constant  rate  of  strain 
machines  where  strain  is  measured  from  the  jaw  travel  of 
the  machine.  Two  important  errors  arise  from  tests  of  this 
type.  First,  there  may  be  slippage  in  the  Jaws  which 
cannot  be  accounted  for  on  the  usual  plot  which  the  machine 
gives  as  output  and,  second,  if  the  whole  length  of  the  spec¬ 
imen  between  the  two  jaws  is  used  as  a  gage  length  there 
will  be  errors  because  of  the  stress  distortions  at  the  ends 
of  the  specimen.  A  uniform  stress  field  will  exist  only 
in  the  centeral  portion  of  the  specimen  which  is  the  region 
far  enough  from  the  ends  to  allow  this  uniform  field  to 
develop. 

The  material  from  which  the 
strips  are  cut  has  a  cold  work  hardness  of  $0%.  It  is  then 
annealed  for  one  half  hour  at  various  temperatures  and  then 
tested.  The  results  of  three  tests  at  a  strain  rate  of 
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1  x  10~*  in/in/sec  are  shown  in  figures  3.2.9,  3.2.10,  and 
3.7.11.  These  results  indicate  that  the  optimum  annealing 
temperature  is  675°C  for  the  time  interval  chosen.  This 
is  therefore  adopted  for  all  the  model  specimens. 

A  comparison  of  the  stress- 

strain  properties  of  phosphor  bronze  with  ASTM  A7  steel  is 
made  in  fig.  3.2.12.  We  see  from  the  two  curves  that 
yielding  takes  place  at  about  the  same  strain.  This  is  very 
important  in  model  studies  since  a  strain  ratio  of  unity 
will  not  result  in  distorted  models.  However,  we  notice  that 
although  phosphor  bronze  shows  some  plastic  yield  at  constant 
stress  it  does  not  possess  the  uniquely  large  plastic  range 
of  structural  steel.  Strain  hardening  begins  much  earlier 
in  phosphor  bronze  than  in  steel  so  that  at  very  large 
strains  the  two  materials  are  not  fully  compatible.  More 
will  be  said  about  this  point  when  actual  test  results  are 
discussed. 

3. 2. 2. 2. 2  Dynamic  Stress-strain  Relations. 
Strips  of  the  same  dimensions  as  in  fig.  3.2.8  are  strained 
at  various  rates  of  strain  to  see  how  the  yield  point  is  af¬ 
fected.  These  specimens  are  also  annealed  at  675°C  for  one 
half  hour  and  then  pulled  to  failure  at  the  rates  indicated 
in  fig.  3.2.13,  The  yield  point  definitely  increases  with 
increasing  rate  of -strain  in  this  material.  Pig.  3. 2. II4. 
shows  the  effect  of  the  strain  rate  on  the  yield  point 
stress  of  phosphor  bronze.  Knowing  the  rate  of  strain  of 
the  model  we  can  estimate  the  yield  stress  of  the  criti¬ 
cally  stressed  points  where  plastici  deformations  will 
take  place.  The  comparison  is  made  in  fig.  3.2.14  between 
structural  steel  and  phosphor  bronze  yield  stresses  as 
affected  by  strain  rate.  It  can  be  seen  that  in  steel  the 
effect  is  more  pronounced,  but  since  the  rates  of  strain 
used  are  of  a  different  range  no  definite  conclusions  can 
be  drawn. 

3.2.3  Using  Plastics  to  Model  Steel  Structures. 

The  possibility  of  using  plastics  as  a  model  material  for 
the  structural  elements  of  steel  is  also  Investigated.  It 
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FIG.3.2.9  STRESS  vs  STRAIN  CURVE  FOR  PHOSPHOR  BRONZE  STRIP 
ANNEALED  AT  650 °C 
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FIG.3.2.10  STRESS  vs  STRAIN  CURVE  FOR  PHOSPHOR  BRONZE  STRIP 
ANNEALED  AT  675 °C 
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FIG.3.2.II.  STRESS  vs  STRAIN  FOR  PHOSPHOR  BRONZE  STRIP 
ANNEALING  TEMPERATURE  AT  7I2°C 


TENSILE  STRAIN- 


TENSILE  STRAIN  x  I0"2  (in/in) 

FIG. 3.2.12  COMPARATIVE  STRESS'  STRAIN  CURVES 
FOR  STRUCTURAL  STEEL  AND 
PHOSPHOR  BRONZE 
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FIG.  3.2.13  EFFECT  OF  RAPID  STRAINING  OF  PHOSPHOR  BRONZE  STRIPS 
ANNEALED  AT  675 °C  FOR  30  MIN. 
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is  observed  that  most  of  the  plastics  are  characterised  by 
a  non-linear  stress-strain  relationship  throughout  the  entire 
range  of  loading,  low  initial  modulus  of  elasticity,  lack 
of  ductility  and  a  marked  degree  of  creep.  Most  plastics, 
in  general,  provide  a  considerable  contrast  to  the  mechani¬ 
cal  properties  of  structural  steel.  However,  unlike  other 
plastics,  ethyl  cellulose,  a  thermo-plastic  material,  has 
a  remarkably  similar  stress-strain  curve  to  that  of  struc¬ 
tural  steel  (ASTM  A7) . 


3.2. 3.1  Properties  of  Ethyl  Cellulose. 

Ethyl  Cellulose,  an  ethyl  ether  of  cellulose,  is  a  thermo¬ 
plastic  material.  Its  chief  characteristics  are  high 
impact  strength,  toughness  and  a  considerable  degree  of 
ductility.  Fluctuations  in  its  dimensions  due  to  normal 
changes  in  temperatures  and  effect  of  humidity  are  quite 
small.  In  spite  of  its  application  in  various  electric 
appliances,  radio  housings,  tooth-brushes  and  other  indus¬ 
tries,  its  manufacturers  are  relatively  very  few. 

Tests  on  standard  specimens  of 
ethyl  cellulose  subjected  to  standard  rates  of  strain 
gave  the  average  stress-strain  curve,  as  shown  in  fig.  3.2.15. 
On  idealising  such  a  stress-strain  relationship,  the  salient 
points  are: 


Modulus  of  Elasticity 
Yield  Stress 
Strain  at  Yield  Stress 
Ultimate  Strain 


=  0.184  x  106  psi 

=  4.15  *  io^  psi 

*  0.0225  in/ in 

=  0.17 2  in/in 


The  specific  gravity  of  ethyl 
cellulose  is  1.10.  A  comparison  of  the  idealized  stress- 
strain  curve  of  ethyl  cellulose  is  made  to  that  of  steel 
in  fig.  3.2.16. 


3.3  MODEL  MANUFACTURE 


The  choice  of  the  materials  from  which  the  model 
is  made  according  to  similitude  requirements  is  the  most 
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FIG. 3.2.15  STRESS  "STRAIN  CURVE  OF  ETHYL  CELLULOSE 


FIG  3.2.16  COMPARISON  OF  STRESS- STRAIN  CURVES  OF  STRUCTURAL  STEEL  AND  ETHYL 
CELLULOSE. 


important  step  in  the  structural  model  analysis.  The 
second  most  important  aspect  is  the  actual  process  of  making 
the  model.  During  this  step  many  considerations  must  be 
evaluated  before  an  efficient  modeling  technique  can  evolve. 

In  this  section  some  of  the  considerations  on  which  the 
decisions  of  the  analyst  depend  will  be  examined.  Also 
the  details  of  the  modeling  technique  for  steel  structures 
will  be  described. 

3.3.1  Some  Considerations  in  Model  Making,  When 
small  scale  models  a he  to  be  built  the  accuracy  of  the 
dimensions  must  be  very  high.  The  difficulty  of  achieving 
this  increases  as  the  scale  of  length  decreases.  At  some 
point  it  may  become  very  hard  or  even  impossible  to  get 
good  results  from  a  small  scale  model.  The  accuracy  to  be 
achieved  in  the  model  also  depends  on  the  results  required 
from  the  model  and  the  time  and  funds  available  for  the 
model  study.  If  the  model  is  to  be  used  to  predict  the 
ultimate  behavior  of  the  prototype,  then  boundary  conditions 
and  other  details,  such  as  joints,  also  become  important. 

If  time  and  money  were  not  limited  then  the  model  could  be 
refined  to  any  degree  with  increasing  chances  of  giving 
more  accurate  results  provided  of  course  that  similitude 
requirements  are  satisfied. 

It  may  also  be  possible  to  use  many  alternate 
materials  for  the  model  especially  in  cases  where  elastic 
behavior  is  to  be  considered.  Some  materials  are  much  easier 
to  work  with  than  others  so  a  careful  consideration  of  all 
available  methods  may  save  considerable  amount  of  time  and 
effort.  Since  in  general  model  studies  are  on  the  high  side 
with  respect  to  time  consumption,  any  effort  on  the  part 
of  the  analyst  to  minimize  this  is  highly  advantageous. 

3.3.2  Specific  Examples  of  Modeling.  We  shall  consider 
here  some  examples  of  modeling  of  steel  structures  in  which 
plastic  behavior  during  the  dynamic  response  is  of  primary 
concern.  Tests  on  actual  full  scale  structures  and  struc¬ 
tural  elements  afford  an  excellent  opportunity  for  comparisons 
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with  extrapolations  from  model  studies.  So  far  little  has 
been  done  in  the  study  of  the  dynamic  response  of  actual 
full  scale  structures  or  structural  elements  although  some 

/  pA  ) 

model  studies  have  been  reported;  '  The  case  of  static 
loads  has  received  more  attention  and  a  host  of  data  from 
actual  full  scale  structures  and  structural  elements  in 
steel  have  been  reported  by  many  investigators  especially 
the  Lehigh  University  group. (22,25,37,41 Before  a 
modeling  technique  can  be  generalized  as  an  efficient 
solution  to  a  class  of  structural  problems,  correlation  tests 
must  be  made  to  prove  its  accuracy  in  predicting  the  behavior 
of  the  prototype.  Comparisons  are  thus  made  with  actual 
tests  on  prototype,  not  with  analytical  solutions,  since 
the  aim  of  the  model  study  is  to  bypass  these.  At  any  rate 
the  behavior  of  the  real  structure  will  not  be  the  same 
as  the  analytical  solution  would  indicate  and  will  be  closer 
to  the  model  prediction  depending  on  the  relative  accuracy 
of  the  basic  assumptions  made  in  each  case. 

3. 3. 2.1  Phosphor  Bronze  Models.  A  decisive 
advantage  of  using  this  material  is  its  ease  of  handling. 
Machining  operations,  fitting  and  especially  jointing  are 
easily  performed  by  the  inexperienced  analyst.  It  should 
be  remembered,  however,  that  a  modeling  technique  is  by 
nature  specific  in  that  one  attempts  to  solve  the  immediate 
problems  of  the  structure  to  be  modeled.  In  some  cases 
the  method  may  be  easily  applied.  To  other  problems  it 
may  not.  Keeping  this  in  mind  then,  we  realize  that  the 
phosphor  bronze  modeling  technique  of  making  models  of 
steel  structures  is  not  claimed  to  be  the  best  nor  the  most 
general . 

3. 3. 2. 1.1  Machining  the  Shapes.  Assuming 
that  most  structures  will  be  designed  from  existing  rolled 
sections  and  not  be  predominantly  built  up,  although  occasion 
may  require  some  such  sections,  we  shall  try  to  duplicate 
through  machining  the  rolled  sections  in  smaller  scales. 

For  convenience  one  section  is  chosen  for  duplication  of 
some  of  the  Lehigh  tests.  The  section  used  mostly  there 
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la  the  8WPI4.O  which  Is  shown  in  fig.  3.3.1(a).  A  scale  of 
1:15  is  chosen  as  the  basis  of  being  the  smallest  practical 
scale  on  which  the  technique  would  be  evolved  and  also 
a  scale  in  which  the  joints  and  boundary  conditions  could 
be  duplicated  faithfully.  In  order  to  facilitate  machining 
and  to  cut  down  on  the  cost  the  exact  dimensions  to  1:15 
scale  were  modified  as  shown  in  fig.  3.3.1(c),  assuming 
that  corrections  could  be  applied  to  the  results  to  over¬ 
come  the  geometric  dlscrepencies.  Machining  of  the  sections 
was  done  in  30  inch  pieces  from  &  inch  square  stock  by  the 
Chicago  firm  of  Milled  Shapes  Inc.  The  tolerances  are  less 
than  +.005M  on  any  dimension.  Machining  always  tends  to 
warp  thin  sections  and  the  machined  phosphor  bronze  pieces 
as  received  showed  some  signs  of  warpage  but  this  was  not 
excessive.  Another  consideration  in  the  choice  of  the 
length  scale  of  1:15  was  the  maximum  capacity  of  the  dynamic 
loading  machine  of  4  kips  and  the  size  of  portal  frame  which 
could  be  handled  by  the  machine  dimensions  with  comparative 
ease. 

3.3. 2.1. 2  Annealing  of  Phosphor  Bronze. 

The  annealing  of  the  machined  shapes  is  carried  out  before 
the  components  are  assembled.  The  optimum  annealing  temper¬ 
ature  of  this  material  is  found  to  be  675°C  for  one  half 
hour,  from  the  stress-strain  curves  on  thin  strips.  Since 
the  annealing  temperature  and  the  time  of  annealing  strongly 
influence  the  yield  point  stress  in  phosphor  bronze,  it  is 
necessary  to  find  a  silver  solder  with  a  lower  flowing 
temperature  than  675°C  to  prevent  double  annealing  upon 
assembly  of  the  parts.  During  annealing  there  may  be 
variations  In  the  temperature  and  this  must  be  guarded 
against  since  only  a  short  time  interval  is  used. 

3.3. 2.1. 3  Silver  Soldering.  The  biggest 
single  advantage  of  using  phosphor  bronze  in  model  making 
is  the  ease  with  which  strong  accurate  joints  can  be  made. 

It  is  this  advantage  that  prompted  the  use  of  this  material 
in  the  small  scale  modeling  of  steel  structures  although 
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(a)  8  WF  40  (NOMINAL) 


(c)  APPROXIMATE  SECTION 
(PHOSPHOR  BRONZE) 


(b)  8  WF  40  (1/15  SCALE) 


FIG.  3.3.1  THE  SECTIONS  OF  THE  PROTOTYPE  AND  MODEL 
STRUCTURES  INVESTIGATED. 
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as  seen  from  the  stress-strain  curves  the  compatibility  of 
the  two  materials  does  not  extend  over  the  whole  strain 
range.  If  the  method  proves  satisfactory  then  small  scale 
model  studies  for  shock  tube  tests  where  size  is  an  impor¬ 
tant  consideration  could  be  made  of  structures  in  the 
1/15  -  l/30  scale  range  very  easily.  The  problems  presented 
by  welding  of  steel  sections  at  these  scales  are  beyond  the 
capabilities  of  the  average  models  laboratory.  Special 
equipment  must  be  available  which  may  hinder  the  making  of 
small  scale  models  of  welded  steel  structures.  In  fact 
is  is  highly  Improbable  that  a  normal  size  weld  as  used 
in  steel  structures  could  be  modeled  at  such  small  scales. 

As  shown  by  fig.  3.3.2  the 

process  of  silver  soldering  is  extremely  simple  with  a  minimum 
of  equipment  necessary.  All  that  is  required  is  a  gas 
torch,  silver  solder  wire  and  a  borax  flux.  The  technique 
can  be  very  easily  adopted  by  technicians  and  the  joints 
can  be  done  with  excellent  control  as  to  thickness  of  weld. 

The  accuracy  to  be  maintained  in  silver  soldering  is  high 
even  for  a  scale  of  1:1$  as  shown  by  fig.  3.3.3  where  the 
stiffeners  which  go  into  the  wide  flange  beams  are  shown 
before  its  assembly  and  below  it  is  a  finished  beam  ready 
for  Instrumentation  and  subsequent  testing. 

It  is  important  in  the  assembly 
to  make  the  joints  as  quickly  as  possible  so  that  over-* 
heating  of  the  member  in  the  vicinity  of  the  weld  will 
be  avoided  and  the  material  will  not  be  weakened.  The 
annealing  temperature  of  the  models  is  bl$°C  and  the  flow¬ 
ing  temperature  of  the  silver  solder  used  was  628°C  so  that 
we  see  that  if  excessive  overheating  occurs  it  could  affect 
the  strength  of  the  model  at  the  local  point  of  soldering. 

3. 3. 2. 2  Fabrication  of  Ethyl  Cellulose  Models. 

It  is,  a  simple  process  to  bond  ethyl  cellulose  to  Itself. 

The  following  considerations  are  necessary  for  a  good  bond: 

1)  Design  and  Fabrication 

ii)  Selection  of  Adhesive 

ill) Application  and  Assembly  Technique 
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PIG.  3.3.2  EQUIPMENT  NEEDED  FOR  SILVER  SOLDERING 
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PIG.  3.3.3  PARTS  WHICH  GO  INTO  THE  MAKING  OF  A  SIMPLE 
WIDE  FLANGE  BEAM  SCALE  1/25. 
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i)  Design  and  Fabrication:  Strength  and  ductility 
are  important  factors  in  the  design  of  joints.  The  tensile 
strength  of  adhesive  bond  is  slightly  lower  than  that  of 
ethyl  cellulose.  Therefore,  it  is  necessary  to  Increase 
the  bonding  area  in  order  that  the  joint  will  be  as  strong 
as  the  molded  part  of  the  material.  This  can  be  done  either 
by  thickening  the  wall  section  around  a  butt  Joint  or  by 
using  a  special  type  of  joint.  Surface  contamination  or 
poor  contact  of  mating  surfaces  should  be  avoided  in  order 
to  obtain  good  bond. 

ii)  Selection  of  Adhesive:  There  are  three  available 
adhesives  especially  suited  for  Ethyl  Cellulose.  They 

are  known  commercially  as  P-34A  and  P-35*»  These  are 

all  the  solvents  for  the  material.  The  first  two  adhesives 
have  the  composition  of  Alcohol  and  Toulene  in  different 
proportions.  The  last  one  is  essentially  Ethylene  Dichloride. 

ill)  Application  and  Assembly  Techniques:  Stronger 
bonds  result  when  the  adhesive  is  applied  to  both  pieces  of 
an  assembly  rather  than  to  only  one  of  the  mating  surfaces. 
This  can  be  explained  by  considering  that,  when  only  one 
piece  is  coated,  the  solvent  first  softens  the  coated  plastic 
and  then  begins  to  thicken  before  the  non-coated  piece  is 
pressed  against  it.  Thus  the  solvent  does  not  have  as  great 
a  chance  to  soften  the  second  piece  and  usually  a  weaker 
bond  results. 

The  adhesive  itself  may  be  applied  by 
various  methods,  such  as  saturated  felt  pad,  flow  gun, 
dripping  or  brush.  High  vapor  concentrations  of  the  adhesives 
should  not  be  inhaled.  Good  ventilation  at  the  place  of 
fabrication  is  desirable,  therefore. 

3.I4.  TESTING  APPARATUS  AND  INSTRUMENTATION 

To  test  the  small  scale  models  made  from  phosphor  bronze 
* 

Dealer  -  Dow  Chemical  Co.,  Midland,  Michigan. 
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special  loading  and  supporting  devioea  are  needed.  The  static 
tests  are  performed  in  the  Models  Laboratory  of  the  Civil 
Engineering  Department  of  M.I.T.  where  a  special  lever 
mechanism  is  available.  The  dynamic  tests  are  performed 
with  the  aid  of  the  impulsive  loading  machine  of  the  Soil 
Mechanics  Division  of  the  same  department.  A  description  of 
this  apparatus  and  the  instrumentation  used  in  the  model 
tests  is  given  below. 

3.4.1  Loading  Frame  for  Static  Tests.  The  loading 
frame  for  carrying  out  static  tests  shown  in  fig.  3.4.1 
has  a  maximum  capacity  of  5  kips  concentrated  load.  At 
the  lower  bound  any  small  load  of  a  few  pounds  only  could  be 
applied  since  the  friction  in  the  system  was  found  to  be 
very  small.  The  frame  is  basically  a  beam  pivoted  at  one 
end  through  a  roller  bearing  passing  through  its  center  of 
gravity  as  shown  in  the  detail  of  fig.  3.4.1.  The  load  is 
applied  in  increments  at  the  furthest  extremity  from  the 
fulcrum  and  the  model  is  positioned  at  any  point  in  between 
the  fulcrum  and  the  point  of  load  application.  By  varying 
the  position  of  the  model  from  the  fulcrum  a  variety  of 
multiplication  factors  can  be  achieved  for  the  applied  load. 

As  point  B  in  fig.  3.4*1  approaches  A  the  multiplication 
factor  increases  with  a  practical  maximum  of  about  10  for 
this  setup. 

Since  most  of  the  structures  to  be  tested  are 
of  a  planar  nature  a  rectangular  box  made  of  aluminum  plates 
is  used  to  support  the  small  scale  models  in  one  plane. 

Fig.  3.4*2  shows  the  dimensions  of  the  box  and  the  method 
of  supporting  of  portal  frames  which  are  tested  to  ultimate 
failure.  Thin  brass  wires  of  1/16"  diameter  are  used  for 
lateral  supports  of  beams  and  frames  to  allow  the  model  to 
reach  its  full  plastic  moment.  It  should  be  noted  that  the 
prototype  structures  which  were  modeled  are  also  supported 
laterally  but  of  course  the  details  vary  in  model  and 
prototype  since  size  limitations  do  not  allow  in  many  cases 
the  reproduction  of  details  in  exactly  the  same  way.  In 
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A.  ROLLER  BEARING  PIVOT  THROUGH  WEB  OF  BEAM  (SECTION  l-l). 

B.  POINT  OF  LOAD  APPLICATION  TO  MODEL. 

C.  ROLLER  BEARING  GUIDES  (SECTION  2-2). 

D.  POINT  OF  LOAD  APPLICATION  AT  END  OF  BEAM. 

E.  COUNTER  BALANCING  WEIGHTS  CAUSING  eg.  PASS  THROUGH  POINT  fi 

F.  ALUMINUM  SUPPORTING  BOX  FOR  MODELS 
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SECTION  l-l 


SYMMETRIC  ABOUT  £ 


FIG.  3.4.1  SCHEMATIC  DRAWING  OF  STATIC  LOADING  FRAME. 
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Fig.  3, 4.3  a  simply  supported  model  beam  Is  being  tested 
In  the  above  mentioned  box.  Note  the  lateral  supports 
given  by  the  brass  wires. 

3. 4. 1.1  Method  of  Load  Application  on  Models. 

Loads  on  beams  and  frames  were  applied  through  small  brass 
bearing  plates  placed  on  the  top  flange.  At  points  of 
load  application  stiffeners  were  provided.  Fig.  3.4*3 
shows  a  close  up  of  the  loading  and  supporting  arrangement 
for  a  simply  supported  beam.  The  concentrated  single 
load  of  the  loading  frame  is  distributed  to  two  or  more 
points  by  small  stiff  beams  as  shown  in  figs.  3.4*4  and  3.4*5. 
The  method  of  load  application  is  the  same  for  some  of  the 
prototype  tests.  However,  in  cases  where  the  load  was 
applied  through  the  web  of  the  beam  an  exact  reproduction 
on  the  model  was  not  attempted  and  all  loads  are  applied 
through  bearing  plates  on  the  top  flange.  Since  the  load 
increments  are  applied  by  hand  at  the  end  of  the  beam  the 
load  variation  with  time  is  a  step  variation  rather  than 
a  smooth  one.  Unloading  of  the  specimen  during  testing  is 
avoided  as  much  as  possible  especially  after  the  model 
passes  its  elastic  limit. 

3.4*1* 2  Deflection  Measurements.  All  deflection 
measurements  for  the  static  tests  are  taken  with  Ames 
dial  gages  with  an  accuracy  of  0.001  inch.  Usually  the 
mid-span  deflection  is  of  particular  interest  in  the  test 
but  in  some  cases  curvatures  are  computed  from  deflection 
measurements.  Figs.  3.47  and  3.43  show  the  mode  of  deflec¬ 
tion  measurements  for  different  model  tests.  Since  the 
scale  of  these  tests  is  1:15  the  deflection  readings  on 
the  model  are  1/15  of  those  of  the  prototype.  The  gages 
used  are  of  sufficient  accuracy  to  allow  good  extrapolation 
to  the  prototype. 

3.4.1. 3  Strain  Measurements.  Strain  measurements 
are  made  by  means  of  SR-4  foil  type  strain  gages  with  an 
epoxy  backing.  Readings  are  taken  with  a  Baldwin-Llma- 
Hamilton  Strain  Indicator.  A  typical  test  setup  is  shown 
in  Fig.  3.4.6. 
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PIG.  3.4.3  CLOSEUP  OF  FRAME  SHOWING  MODEL  BEAM  READY  FOR  TEST 


FIG.  3.4.4  CLOSEUP  VIEW  OF  THE  LOAD  APPLICATION  AND  THE 
SUPPORT  OF  A  SIMPLY  SUPPORTED  MODEL  BEAM 


PIG.  3.4-5  PORTAL  FRAME  IN  POSITION  TO  BE  TESTED.  NOTE 
THE  LOAD  APPLICATION  AND  SUPPORT  DETAILS. 


„  'X 


PIG.  3.4.6  CLOSEUP  OP  SUPPORTING  BOX  SHOWING  THE  LATERAL 
SUPPORTS  FOR  A  PCRTAL  FRAME,  THE  METHOD  OP 
LOAD  APPLICATION  AND  THE  STRAIN  INDICATOR. 
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FIG.  3.4.7  CLOSEUP  OF  LOADING  FRAME  SHOWING  A  PORTAL 

FRAME  READY  FOR  TEST.  NOTE  LATERAL  SUPPORTS. 


FIG.  3.4.8  DEFLECTION  MEASUREMENTS  OF  JOINT  TEST  WITH 
AN  AMES  DIAL. 
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The  foil  gages  used  varied  in  length 
from  3/16”  to  1"  and  are  bonded  with  EFY-150  cement 
according  to  the  recommended  procedure  of  the  manufacturer. 
This  cement  is  chosen  for  its  relatively  fast  cure.  However, 
from  observations  it  appears  that  some  of  the  gages  began 
to  creep  at  strains  much  less  than  the  maximum  range  of 
about  1 %.  This  is  thought  to  be  caused  by  failure  of  the 
bond  between  the  gage  and  the  model  material  and  measurements 
on  such  gages  should  be  taken  only  after  the  creeping  has 
ceased.  Whereas  the  metal  material  responds  quickly  to  an 
increment  of  stress,  the  same  is  not  true  of  the  epoxy 
which  is  a  thermo  setting  plastic.  When  enough  time  is 
allowed  for  creep  to  occur  the  strain  measurements  give 
better  results. 

The  gages  are  temperature  compensated 
for  1018  Steel  but  since  these  are  used  with  phosphor  bronze 
with  a  different  coefficient  of  expansion  a  dummy  gage  is 
used  in  all  the  tests. 

3.4*2  The  Dynamic  Loading  Machine.  The  machine  used 
to  test  the  model  structures  is  one  built  for  high  speed 
triaxial  testa  of  soil  samples.  Pig.  3*4*9  shows  a  drawing 
of  the  machine  and  Pig.  3.4*10  shows  its  mode  of  operation 
in  a  schematic  form.  Since  the  system  had  no  feedback  control 
on  the  load  rate,  the  obtained  load  time  curve  is  the  result 
of  interaction  between  the  specimens  and  the  machine  and 
as  such  the  applied  load  is  an  impact  rather  than  a  predeter¬ 
mined  impulse. 

Loads  up  to  5  kips  are  obtained  with  this 
machine  and  load  rise  times  ranging  from  3  milliseconds  up. 

3.4, 2.1  Support  Frame.  A  simple  braced  frame 
was  built  to  hold  the  specimens  to  be  tested.  Fig.  3.4*11 
shows  the  frame  with  a  model  phosphor  bronze  portal  frame 
ready  to  be  tested.  The  support  ^rame  is  clamped  to  the 
platform  of  the  machine  during  the  test. 
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FIG.  3.4.9  THE  DYNAMIC  LOADING  MACHINE  USED  IN  MODEL 
TESTING. 
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FIG.  3.4.10  BASIC  UNITS  OF  MODEL  LOADING  MACHINE 
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3.U.2.2  The  Load  Cell,  The  load  cell  is  attached 
to  the  ram  and  the  load  applied  to  the  model  through  it.  It 
conaiated  of  a  J  inch  diameter  aluminum  shaft,  2  inches 
long  onto  which  are  bonded  two  strain  gages  (  Baldwin-Lima- 
Hamilton  Type  C?) .  The  gages  are  laid  along  diametrically 
opposed  generators  of  the  cylindrical  shaft  and  wired  to 
form  two  opposite  arms  of  a  Wheatstone  bridge.  The  other 
two  arms  of  the  bridge  consisted  of  two  dummy  gages  attached 
inside  an  unloaded  aluminum  box.  This  system  caused  the  por¬ 
tions  of  strains  due  to  bending  to  cancel  each  other  out  and 
the  bridge  output  voltage  to  be  proportional  to  the  axial 
load  in  the  load  cell.  The  circuit  is  excited  by  a  22J 
volt  battery  and  the  output  signal  diaplayed  on  a  cathode 
ray  oscillograph. 

3.I4..2.3  Recording  Equipment.  The  oscillograph 
as  shown  in  Pig.  3.U.12  is  a  Tektronics  Type  502  which  had 
a  calibrated  time  sweep  and  could  be  triggered  by  the  load 
signal  Itself.  Permanent  records  are  obtained  by  using 
Dumont  Type  287  oscillograph  record  cameras  operating  with 
Polaroid-land  film. 

A  Twin  Vise  Model  60-1300  Sanborn  Recorder 
is  used  during  calibration  of  the  load  cell. 

The  load  cell  is  calibrated  in  a  standard 
beam  type  testing  machine,  the  output  signal  being  recorded 
on  a  Sanborn  unit.  The  signal  due  to  a  known  load  is  com¬ 
pared  to  the  signal  due  to  a  fixed  unbalance  in  a  Wheatstone 
bridge.  After  each  test  the  signal  due  to  this  fixed  unbal¬ 
ance  would  be  recorded  giving  a  load  scale  for  the  data 
obtained. 

The  deflections  are  measured  with  a  Linear 
Variable  Differential  Transformer  (L.V.D. T),  the  signal  is 
converted  to  D.C.  and  displayed  on  the  Oscillograph.  Thus 
two  signals  are  recorded  with  time  on  the  same  picture.  For 
any  time  during  the  response  the  load  deflection  relation 
would  be  obtained. 


138 


FIG.  3.4.11  SUPPORTING  FRAME  FOR  MODEL  FRAMES.  NOTE  THE 
L.V.D.T.  FOR  THE  DEFLECTION  MEASUREMENTS. 


FIG.  3.4.12  SETUP  SHOWING  THE  DYNAMIC  LOADING  MACHINE 
AT  THE' LEFT  AND  INSTRUMENTATION. 
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The  L.V.D.T.  is  calibrated  with  an  Ames 
dial  and  a  plot  of  signal  v.s.  deflection  is  obtained  for 
every  test.  From  this  plot  the  deflections  for  the  response 
signal  recorded  by  the  camera  could  be  obtained. 


CHAPTER  4 


PHOSPHOR  BRONZE  MODEL  STUDIES  OF  WELDED  STEEL  STRUCTURES 


Tests  were  made  on  wide  flange  beams,  portal  frames, 
and  connections  for  the  static  loads,  and  portal  frames 
in  the  case  of  dynamic  loads.  In  all  these  tests  a  scale 
of  1:15  was  used  for  convenience  and  the  section  used 
closely  approximates  an  8WF40  rolled  section.  In  the 
static  tests  comparisons  are  made  between  the  model  predic¬ 
tions  and  the  actual  prototype  tests  but  in  the  dynamic 
case  no  prototype  results  were  available. 

4.1  STATIC  TESTS  ON  MODELS  BEYOND  THE  ELASTIC  LIMIT 

4.1.1  Wide  Flange  Beams.  A  series  of  seven  tests 
were  performed  on  wide  flange  beams  with  simple  supports 
and  fixed  ends  from  which  load  deflection  and  moment  curva¬ 
ture  relations  were  obtained.  The  test  results  of  the 
model  beams  are  summarized  in  table  4.1.  Using  the  full 
scale  results  of  tests  performed  at  Lehigh  University  the 
predicted  behavior  of  the  prototype  from  the  model  test  is 
compared. 

4.1. 1.1  Prototype  Tests.  The  simply  supported 
wide  flange  beams  are  of  two  different  types.  The  differ¬ 
ence  is  in  both  length  and  mode  of  loading  as  shown  in 
Pigs.  4.1.1  and  4.1.2. 

The  "pilot"  tests  were  made  on  an  8WF40 
beam  of  12-ft.  span,  simply  supported  and  loaded  at  the 
third  points.  The  load  was  applied  to  the  beam  through 
bearing  blocks  resting  on  the  top  flange  and  the  web  was 
stiffened  under  the  load  with  3/8  in.  plates  welded  per¬ 
pendicular  to  the  web.  Prom  this  test  it  was  observed  that 
a  non-linear  strain  variation  exists  at  the  mid-span  after 
the  elastic  limit  was  exceeded.  This  caused  the  neutral 


TABLE  U.l  SUMMARY  OF  BEAM  TESTS 
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This  Is  the  value  for  the  as  delivered  beam.  The  value  for  the  annealed  beam  is 
U9  kips. 

This  is  the  value  for  support  condition  (a)  as  shown  in  Fig.  (14.1.12).  The  value 
for  support  condition  (b)  is  108  kips. 
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PIG.  4,1.1  PROTOTYPE  PILOT  TEST  OP  8WF40  SIMPLY  SUPPORTED 
BEAMS  LOADED  ON  THE  FLANGE  AS  SHOWN  PERFORMED 
AT  LEHIGH  (REF.  25  ). 
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PIG.  4.1.2  PROTOTYPE  8WFI4.O  SIMPLY  SUPPORTED  BEAMS  WITH  LOADS 
APPLIED  THROUGH  THE  WEB  AS  SHOWN  ABOVE  TESTED 
AT  LEHIGH  (REP.  25  ). 
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axis  to  move  towards  the  tension  flange  but  at  higher  per 
cent  strains  it  moved  back  towards  the  centroidal  axis. 

The  pilot  test  developed  less  plastic  strength  than  that 
calculated. 

The  "regular”  tests  were  made  on  8WFij.O 
beams  of  ll(.-ft.  span  simply  supported  and  loaded  at  the 
third  points.  One  test  was  done  on  the  beam  in  the  as 
delivered  condition  and  the  other  was  stress-relief  an¬ 
nealed  prior  to  testing.  The  method  of  load  application 
was  made  to  simulate  a  beam-to-girder  connection  by  bring¬ 
ing  the  load  directly  onto  the  web.  This  method  of  load 
application  could  not  be  duplicated  in  the  model  tests 
without  undue  loss  of  time.  Thus  the  method  of  load  appli¬ 
cation  in  all  the  model  tests  was  done  through  bearing 
blocks  on  the  top  flange  as  shown  in  Fig.  I4..I.3.  The  beams 
were  laterally  braced  to  prevent  lateral  bucking  upon  passing 
of  the  elastic  limit. 

From  the  results  of  these  simply  supported 
beams  the  agreement  between  the  theoretical  and  experimental 
results  was  better  for  the  annealed  beam  in  the  "regular" 
test.  The  "pilot"  beams  were  not  annealed  nor  laterally 
supported  and  also  loaded  on  the  top  flange  through  bearing 
blocks.  These  may  be  the  causes  of  the  lower  than  predicted 
moment  of  the  pilot  tests.  Annealing  although  lowering  the 
residual  stresses  considerably  also  lowered  the  yield  stress 
and  thus  the  strength  of  the  annealed  specimens  was  less. 

It  is  concluded  from  the  test  results  that  the  moment  curva¬ 
ture  curves  for  structural  steel  beams ,  laterally  supported, 
can  be  predicted  from  the  stress-strain  tensile  test  curves 
and  the  usual  theory  of  plastic  bending  based  on  the  assump¬ 
tions  of  a  linear  strain  variation  across  the  section  and  a 
uniform  distribution  of  yield. 

The  prototype  for  the  continuous  beams 
were  also  of  ll*-ffc.  span.  However,  the  "fixity"  in  these  tests 
results  from  the  2  side  spans  of  7-ft.  each  of  a  3  span 
continuous  beam.  The  central  span  of  ll|.-ft  is  considered 
fixed  ended.  Two  types  of  tests  are  reported  from  the 


Lehigh  series  on  continuous  beams. One  has  the  reaction 
applied  to  the  flange  through  a  bearing  plate  and  the  other 
the  reaction  is  transmitted  to  the  web  by  a  beam  to  column 
simulated  connection.  The  test  results  indicate  that  the 
latter  support  detail  has  more  stiffness  and  thus  the  beam 
mobilized  higher  loads. 

4. 1.1. 2  The  Model  Tests.  The  models  of  the  simply 
supported  beams  were  of  1:15>  scale  and  are  soown  in  Pigs. 

4.1.3  and  4.1.4*  The  cross-sectional  geometry  was  only 
approximately  scaled,  however,  due  to  machining  difficulties 
(see  Appendix  C).  The  phosphor  bronze  sections  machined 
from  square  stock  were  first  annealed  as  described  earlier 
and  then  assembled.  Pig.  4.1.5  3hows  a  beam  before  testing 
and  below  it  a  failed  simply  supported  beam,  failure  due  to 
excessive  deflection.  In  beam  tests  B1  to  b5  only  centerline 
deflections  were  recorded  with  load  increase.  Beams  B1A  and 
BIB  were  instrumented  with  strain  gages  at  the  mid-span 
section  and  also  in  BIB  deflections  were  measured  at  2  other 
points  so  that  moment  curvature  relations  could  be  established 
from  the  model  test.  The  position  of  the  strain  gages  is  shown 
in  Fig.  4.1.3.  The  beams  B1A  and  BIB  are  shown  in  Pig.  4.1.6 
after  failure.  Beam  B2  and  B3  are  shown  in  Pig.  4.1.7. 

For  the  fixed  ended  beams  B4  and  B5  the 
fixity  was  achieved  in  the  model  by  silver  soldering  the 
ends  to  a  tube  in  the  case  of  B4  and  pouring  molten  lead 
in  a  similar  tube  in  the  case  of  beam  B$.  The  ends  were  then 
fastened  to  the  supports  with  two  screws  on  each  support. 

The  detail  of  the  model  "fixed"  support  is  shown  in  Fig.  4.1.12. 
The  fixed  ended  model  beams  are  shown  after  failure  in 
Fig.  4.1.8. 

4. 1.1. 3  Test  Results.  The  results  of  the  model 
beams  representing  the  "pilot"  test  of  the  prototype  are 
given  in  Tables  4-2,  4.3,  and  4*4.  These  are  for  model  beams 
Bl,  B1A  and  BIB  respectively.  The  explanations  following  each 
of  these  tables  gives  the  procedure  followed  in  obtaining  the 
numerical  values  in  each  column.  Tables  4.5  and  4.6  give  the 
test  results  of  model  beams  B2  and  B3  and  Tables  4.7  and  4.8 
give  the  results  on  the"fixed"  end  model  beams  B4  and  B£. 


* -  9.6  - 

knife  edge 


FIG.  4.1.3  MODEL  BEAM  OF  THE  "PILOT"  PROTOTYPE  TEST 
AT  1:15  SCALE. 


FIG.  4.1.4  MODEL  BEAM  OF  THE"  REGULAR"  PROTOTYPE  TEST 


AT  1:15  SCALE. 


PIG.  4.1.5  MODEL  BEAM  B1  AFTER  FAILURE. 

MATERIAL,  PHOSPHOR  BRONZE. 


SCALE  lilSt 


I 


BIB 


FIG. 


k.1.6  MODEL  BEAMS  B1A  AND  BIB  AFTER  FAILURE.  NOTE 
THE  STRAIN  GAGES  AT  THE  CRITICAL  SECTION. 
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S33SBBC 


PIG.  4.1.7  SIMPLY  SUPPORTED  MODEL  BEAM  l:l£  SCALE 
AFTER  FAILURE. 


FIG.  4.1.8  FIXED  ENDED  MODEL  BEAM  AFTER  FAILURE. 
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Comparison  of  these  results  with  the  prototype 
tests  is  given  in  Fig.  4.1.9  for  the  "pilot"  test  or  the 
8WF40  simply  supported  beam  of  12-ft.  span.  This  figure  shows 
the  M-^S  curves  obtained  from  strain  gage  readings  and  also 
deflection  readings  in  the  case  of  beam  BIB. 

Fig.  4- 1.10  is  a  load  deflection  curve  for  the 
predicted  behavior  of  the  prototype  from  3  model  tests.  This 
shows  the  method  that  could  be  used  to  predict  ultimate 
failure  loads  in  more  complicated  structures  from  a  few  model 
tests.  From  this  figure  the  predicted  failure  load  is 
conservatively  picked  at  56  kips.  Backf iguring  the  ultimate 
moment  of  Fig.  4-1.9  (the  load  deflection  curve  is  not  given 
by  LUxion  et  al.)^^  we  get  a  value  of  52  kips  for  the 
ultimate  load. 

Fig.  4-1-11  shows  the  load  deflection  curves 
from  model  beams  B2  and  B3  together  with  the  prototype 
results . 

Fig.  4-1-12  shows  the  load  deflection  pre¬ 
diction  curves  from  model  tests  B4  and  B5  as  compared  to 
two  prototype  tests  of  different  boundary  conditions. 

4. 1.1. 4  Conclusions  from  the  Model  Comparisons. 

From  the  few  tests  which  are  reported  herein  a  number  of 
conclusions  could  be  drawn  concerning  the  correlation  of 
model  predictions  and  actual  prototype  tests. 

1.  It  is  obvious  from  all  the  load  deflection  and 
the  moment  curvature  plots  that  the  model  material  strain 
hardens  at  earlier  values  of  strain  than  steel.  This  was 
anticipated  and  the  tests  verify  it.  However,  as  shown 

by  Figs.  4-1-11  and  4-1-12  the  ultimate  loads  on  these  beams 
could  have  been  predicted  within  the  accuracy  required  for 
design  purposes. 

2.  The  results  of  the  two  simply  supported  model 
beams  B2  and  B3  which  modeled  the  "regular"  tests  of  14-ft. 
spans  indicate  that  the  model  results  approach  closer  the 
annealed  8WF40  beam  rather  than  that  as  delivered  as  shown 
by  the  values  of  Table  ,4-1-  The  failure  criterion  chosen 
for  the  beams  is  of  excessive  deflections  due  to  bending 
deformations. 
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MOMENT  IN  KIP  -  INCHES 


FIG. 4.1.9  SIMPLY  SUPPORTED  BEAMS  M"$  RELATIONS 
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TOTAL  LOAD, 2P,  (KIPS) 


MIDSPAN  DEFLECTION  IN  INCHES 


FIG.4.UI  LOAD  DEFLECTION  CURVES  FOR  SIMPLY 
SUPPORTED  WIDE  FLANGE  BEAMS 
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TOTAL  LOAD, 2 P,  (KIPS) 


FIG. 4.1.12  LOAD  DEFLECTION  CURVES  FOR  FIXED 
ENDED  WIDE  FLANGE  BEAMS 
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TABLE  4.2  TEST  RESULTS  FOR  MODEL  BEAM  B1 


r  . .  —  — Hr 

; 

1 

j  Load  at  j 

j  End  of  Beam 
(Lbs.) 

Mid-Span 

Deflection 

( Inches) 

Total  Load 

2P 

on  Model 

(Lbs. ) 

Predicted* 

Prototype 

Load 

(Kips) 

Predicted 

Prototype 

Deflection 

( inches) 

1 

.004 

5.92 

3.88 

.06 

2 

.006 

11.84 

7.74 

.09 

3 

.008 

17.76 

11.6 

.12 

l 

4 

.01 

23.68 

15.5 

.15 

5 

L 

.013 

29.60 

19.3 

.195 

O 

7 

.019 

41.144- 

27.0 

.285 

8 

.0215 

47.36 

30.9 

.324 

9 

.024 

53.28 

34.8 

.36 

10 

.028 

59.20 

37.4 

.42 

11 

.031 

65.12 

42.8 

.465 

12 

.035 

71.04 

46.4 

.525 

13 

.0395 

76.96 

50.2 

.593 

14 

.046 

82.88 

54.0 

.69 

15 

.069 

88.80 

57.7 

1.035 

16 

.148 

94.72 

61.8 

2.22 

17 

.267 

100.64 

65.5 

4.00 

18 

.329 

106.56 

69.4 

4.94 
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Explanation  of  Table  4.2. 

Col.  1  gives  the  load  increments  on  the  beam  of  the  load¬ 
ing  frame. 

Col,  2  is  the  midspan  deflection  read  on  the  Ames  dial. 

Col.  3  gives  the  total  load  on  the  model  beam  or  the 
sum  of  the  2  third  point  concentrated  loads.  Col.  3  gives  the 
values  of  Col.  1  times  5.92,  the  ratio  of  lengths  on  the 
beam  or  the  load  multiplication  factor  of  the  loading  frame. 

Col.  4  is  the  predicted  load  of  the  prototype.  It  is 
obtained  from  the  following  relation: 

P„  =  x  f  "2  x  f  "1  x  (G.P.)  (4.1.D 

p  m  r  9 


where 

f  is  the  length  scale 

fg.  is  the  ratio  of  stresses 

G.P.  is  a  geometric  factor  to  account  for  the 
difference  in  moment  of  inertia  as  explained  in  Appendix  c, 

fr'2  x  f^'1  x  G.P.  =  (15)2  x  2  x  1.45  =  652 


P„  =  .65 2  Pm  (kips) 

p  m 


(4.1.2) 


Cel.  5  gives  the  predicted  defection  at  mid  span  of  the 
prototype.  The  values  are  obtained  from  the  expression: 


U  f 
m  r 


-1 


(4.1.3) 


where 

Up  is  a  vector  displacement  in  the  prototype. 

U  is  a  vector  displacement  at  the  same  point  and 
direction  in  the  model. 

r  is  the  length  scale  which  is  l/l5  in  this  case. 
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TABLE  I4..3  TEST  RESULTS  OP  MODEL  BEAM  BlA 
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TABLE  4.3  (CONT. )  TEST  RESULTS  OP  MODEL  BEAM  B1A 
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TABLE  4.3  (CONT. ).  TEST  RESULTS  OP  MODEL  BEAM  B1A 
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Col,  1  -  Col,  5  are  the  same  as  In  Table  4.2. 

Col,  6  gives  the  midspan  moment  in  the  model.  This  is 
obtained  from  £  x  Col.  2  x  L/3,  where  L/3  =  3.2”. 

Col.  7  gives  the  midspan  moment  in  the  prototype.  The 
expression  for  this  value  is: 

Mp  =  Mm  *  fr"3  *  V1  *  G*F*  (4.1.4) 

where,  f  =  the  length  scale 

f  =  the  ratio  of  stresses 

G.P.  =  the  correction  factor  due  to  discrepencies 
in  the  moment  of  inertia  as  given  in 
Appendix  C 

Mp  =  Prototype  moment  at  midspan. 

M  =  Model  moment  at  midspsn 
m 

fr‘3  x  f^'1  x  G.P.  =  (1$) 3  x  2  x  1.45  *  976 


therefore, 

Mp  =  .976  Mn  (Kip-in)  (4.3.5) 

Col.  8  -  Col.  11  give  the  strain  as  read  on  the  strain 
indicator  from  one  SR4  strain  gage. 

Col.  12  gives  the  curvature  or  the  rotation  per  unit 
length  at  midspen.  The  expression  used  to  calculate  this 
as  derived  in  Appendix  D  from  the  strain  measurements: 


(4.1.6) 


Col.  13  gives  the  predicted  curvature  or  rotation  per 
unit  length  of  the  prototype  at  midspan.  It  is  obtained 
from 

<1/R)p  =  <l/H>B*fr  (4.1.7) 

where  the  subscripts  indicate  prototype  and  model  respectively 
and  f  is  the  length  scale. 
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TABLE  4.4.  TEST  RESULTS  OP  .MODEL  BEAM  BIB 
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Explanation  of  Table  lull.. 

Col.  1  -  Col.  13  are  obtained  In  the  same  manner  as 
indicated  for  Table  1^.3. 

Col.  14  is  the  curvature  or  the  rotation  per  unit 
length  at  midspan  as- obtained  from  the  3  measured  deflec¬ 
tions.  Since  the  moment  is  constant  over  the  middle  third 
of  the  beam,  the  curvature  is  constant  and  has  the  shape 
of  a  circular  arc.  The  rate  of  change  of  slope  at  the 
center  line  is  then  given  by  the  following  expression  as 
derived  in  the  Appendix  E 


(Axp  (I4..I.8) 

=  the  deflection  at  left  third  point 

=  the  deflection  at  the  center  line 
=  the  deflection  at  the  right  third  point 

*  One  sixth  of  span  of  the  beam. 


where. 
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£c 

£r 

Ax 
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TABLE  4.5.  TEST  RESULTS  OP  MODEL  BEAM  B2 


1 

2 

3 

— 

k 

£ 

Load  at  End 

of  Beam 

(Lbs.) 

Midspan 

Deflection 

( Inches) 

Total 

Load  on 

Model 

(Lbs.) 

Predicted 

Total  Load 

on  Prototype 
(Kips) 

Predicted 

Prototype 

Midspan 

Deflection 

( Inches) 

1 

.0115 

3.88 

.172 

2 

.017 

7.74 

.255 

3 

.022 

17.76 

11.6 

.33 

4 

.027 

23.68 

15.5 

.405 

5 

.0315 

29.60 

19.3 

.472 

6 

.036 

35.52 

23.2 

.54 

7 

.042 

41.44 

27.0 

.63 

8 

.046 

47.36 

30.9 

.69 

9 

.052 

53.28 

34.8 

.78 

10 

.0575 

59.20 

37.4 

.862 

11 

.067 

65.12 

42.8 

1.005 

12 

.085 

71.04 

46.4 

1.275 

13 

.170 

76.96 

50.2 

2.55 

14 

.265 

82.88 

54.o 

3.98 

15 

.3515 

88.80 

57.7 

5.27 

16 

.464 

94.72 

61.8 

- k 

6.96 
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TABLE  I4..6.  TEST  RESULTS  OP  MODEL  BEAM  B3 


1 

1 - r 

2  1 

1 

3 

k 

5 

i 

Predicted 

Load  at  End 

Midspan  j 

Total 

Predicted 

Prototype 

of  Beam 

Deflection  ! 

Load  on 

Total  Load 

Midspan 

1 

1 

Model 

on  Prototype 

Deflection 

(Lbs. ) 

(Inches) 

1 

(Lbs.) 

(Kips) 

( Inches) 

1 

.005 

5.92 

3.88 

.075 

2 

.01 

11.81+ 

7.71+ 

.15 

3 

.015 

17.76 

11.6  j 

.225 

1+ 

.019 

23.68 

15.5 

.287 

5 

.0225 

29.60 

19.3 

.336 

6 

.027 

35.52 

23.2 

.1+06 

7 

.0315 

1+1.14+ 

27.0 

.1+73 

8 

.036 

1+7.36 

30.9 

.51+ 

9 

.01+1 

53.28 

31+.8 

.615 

10 

.01+6 

59.20  ! 

37.1+ 

.68 

11 

.052  , 

65.12  ; 

1+2.8 

.78 

12 

.063 

71.01+ 

1+6.1+ 

.95 

13 

.107 

76.96  ' 

50.2 

1.61 

11+ 

.255 

82.88  ! 

51+.0 

3.38 

15 

.1+03 

88.80 

57.7 

6.06 

16 

.592 

91+.72 

61.8 

8.90 
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TABLE  4.7  TEST  RESULTS  OP  MODEL  BEAM  B4 


i 

1 - 

1 

2 

3 

|  k 

1 

5 

Load  at  End 
of  Beam 

Midspan 

Deflection 

Total  Load 
on  Model 

Predicted 
Total  Load 
On 

Predicted 

Midspan 

Deflection 

(Lbs.) 

( Inches ) 

(Lbs.) 

Prototype 
i  (Kips) 

Prototype 
( Inches) 

1 

.002 

5.92 

3.88 

.03 

2 

.003 

11.84 

I  7.74 

.045 

3 

.0045 

17.76 

11.6 

.0676 

4 

.006 

23.68 

15.5 

.09 

5 

.008 

29.60 

!  19.3 

.12 

6 

.0095 

35.52 

!  23.2 

.1425 

7 

.0115 

41.44 

|  27.0 

.1725 

e 

.013 

47.36 

30.9 

.195 

9 

.015 

53.28 

34.8 

.225 

10 

.017 

39.20 

37.4 

.255 

11 

.019 

65.12 

42.8 

.285 

12 

.021 

71.  Oli 

46.4 

.315 

13 

.023 

76.96 

50.2 

.345 

14 

.025 

82.88 

54.0 

.375 

15 

.029 

88.80 

57.7 

.436 

16 

.031 

94.72 

61.8 

.465 

17 

.033  ! 

100.64 

65.5 

.496 

18 

.036 

106.56 

69.6 

.54 

19 

0 

0 

• 

112.  48 

74.0 

.60 

20 

.044 

118.40 

77.1 

.66 

21 

.049 

121|.32 

81.2 

.736 

22 

.057 

130.24 

85.0 

.856 

23 

.065 

136.16 

89.0 

.975 

24 

.085 

142.08 

92.8 

1.275 

25 

.117 

148.00 

96.6 

1.755 

26 

.163 

153.92 

100.0 

2.44 

27 

.206 

159.84 

104.0 

3.09 

28 

.272 

165.76 

107.8 

4.08 

29 

.363 

171.68 

111.8 

5.45 
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TABLE  I4..8  TEST  RESULTS  OP  MODEL  BEAM  B5 


16? 


p 


p 


<L 


MIDSPAN  DEFLECTION  IN  INCHES  — •> 


FIG. 4.1.13  CRITERION  FOR  ULTIMATE  LOAD  PREDICTION 


3.  The  agreement  In  the  case  of  the  two  ''fixed”  ended 
beams  is  closer  to  the  prototype  test  with  support  conditions 
(a)  as  shown  in  Pig.  4., 1.12.  Since  the  model  boundary  con¬ 
ditions  in  the  model  are  simulated  in  a  different  way  it 

is  hard  to  draw  any  conclusions  as  to  which  type  of  boundary 
conditions  the  model  beams  approached.  However,  it  should 
be  noted  that  there  was  some  rotation  at  the  supports  of  the 
model  and  the  extra  heat  applied  in  the  manufacture  of  the 
model  at  the  two  ends  may  have  weakened  the  material  at 
those  two  points  considerably  more  than  the  yield  stress 
of  18,000  psi  of  the  regular  annealing  process. 

4.  The  criterion  used  in  picking  the  ultimate  loads 
from  the  load  deflection  curve  will  be  illustrated  with 
the  idealized  bi-linear  load  deflection  curve  of  beam  B2 
shown  in  Pig.  4.1.13.  At  the  intersection  of  the  two  linear 
approximations  to  the  predicted  load  deflection  curve  we 
read  a  total  load  of  4&  Kips  and  a  midspan  deflection  of 

1  inch.  This  is  actually  the  minimum  value  of  a  so-called 
"ultimate  load"  since  at'r  larger  deflection  the  correlation 
between  the  model  test  and  the  prototype  test  will  be  better. 
In  fact  at  a  midspan  deflection  of  2.25  inches  the  total 
ultimate  load  is  the  same  for  both  or  48.5  Kips.  Thus  if 
the  maximum  allowed  deflection  of  the  prototype  structure  is 
known  (which  is  the  case  in  most  structural  problems)  then 
it  is  a  simple  matter  to  pick  the  ultimate  load  from  the  pre¬ 
dicted  load  deflection  curve  obtained  in  a  model  study. 

4.1.2  Portal  Frames.  In  order  to  test  the  modeling 
technique  better  in  cases  of  real  structures  with  more  elab¬ 
orate  boundary  conditions,  it  was  decided  to  make  model 

tests  of  singly  redundant  portal  frames  as  tested  at  full 

X  37) 

scale  at  Lehigh  University  .  The  results  from  the  model 
test  were  then  used  to  predict  the  behavior  of  the  prototype 
and  a  comparison  was  made  between  this  prediction  and  the 
actual  prototype  results. 

4. 1.2.1.  Prototype  Tests.  The  dimensions  of  the 
prototype  and  the  loading  set  up  are  shown  in  Fig.  4.1.14  • 
The  rolled  section  used  for  beams  and  columns  are  the  8WF40 
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PIG.  4..I.U4.  TEST  SET  UP  FOR  THE  PROTOTYPE  TEST  ON  PORTAL 
FRAMES  (REP.  37  ). 


FIG.  i4-.l-.i5  HORIZONTAL  REACTION  ASSEMBLY 


section  used  before  in  the  wide  flange  beams.  The  loads 
applied  were  two  vertical  concentrated  loads  at  the  3/8 
points.  Side-sway  and  lateral  buckling  were  prevented  In 
these  tests.  The  boundary  conditions  used  are  shown  in 
Fig.  4*1*14*  They  simulate  pin  ended  supports  at  the  base 
of  the  columns.  The  detail  of  the  horizontal  reaction  meas- 
ureing  assembly  Is  shown  in  Fig.  4.1.15.  Side-sway  was 
prevented  by  a  longitudinal  support  on  one  of  the  knees. 

Lateral  supports  were  provided  at  the  corner  connections  and 
at  4  points  along  the  length  of  the  beam  where  severe  plastic 
straining  was  expected.  The  loads  were  transferred  to  the 
web  of  the  beam  as  in  the  case  of  the  beam  tests  described 
in  the  previous  section. 

4. 1.2. 2  The  Model  Tests.  The  model  frame  was 
fabricated  from  the  phosphor  bronze  section  shown  in  Fig.  4* 1.16 
and  was  at  a  scale  of  1:15  of  the  prototype.  Lateral 
supports  were  provided  by  means  of  thin  brass  wires  attached 
to  the  stiffeners  of  the  beam  and  knees  as  shown  in  Fig.  4*1*17. 
The  joints  used  are  a  modified  version  of  the  prototype. 

This  and  other  types  of  connections  will  be  discussed  in 
a  section  that  follows.  The  horizontal  reaction  was 
measured  by  reading  the  strain  on  a  stiff  steel  bar  1"  x  .03" 
in  cross  section.  The  bar  together  with  other  details  of 
the  test  can  be  seen  in  Fig.  4*1*17.  The  position  of  the 
strain  gages  is  shown  in  Fig.  4*1*16  together  with  the 
horizontal  reaction  measuring  bar. 

Loads  were  applied  to  the  frame  by  means 
of  bearing  plates  at  the  top  flange  as  seen  in  Fig.  4*1*17. 
Side-sway  was  prevented  by  attaching  a  thin  brass  wire  to 
the  right  hand  side.  However,  this  proved  to  be  Inadequate 
under  very  large  deflections  and  the  ultimate  failure  of  the 
frame,  after  considerable  plastlficatlon  took  place,  was  a 
side-sway  failure  to  the  right.  This  is  shown  very  clearly 
in  Figs.  4*1*19  and  4*1*20  of  the  failed  specimen. 

Deflections  were  measured  at  midspan  by 
means  of  an  Ames  dial  as  shown  in  Fig.  4*1*18. 
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SECTION  A-A 


FIG.  4.1.16  MODEL  FRAME,  1:15  SCALE,  SHOWING  HORIZONTAL 
REACTION  MEASURING  DEVICE  AND  POSITION  OF 
STRAIN  GAGES. 
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FIG.  4.1.17  POSITIONING  OF  THE  MODEL  PORTAL  FRAME  IN 
SUPPORTING  BOX,  SHOWING  LOADING  METHOD. 


FIG.  4.1.18  TEST  SET  UP  FOR  MODEL  FRAME 
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PIG.  4.1.19  MODEL  FRAME  AFTER  TEST. 


FIG.  4.1-20  MODEL  FRAME  SHOWING  A  LATERAL  SWAY  FAILURE 
INITIATED  BY  LOCAL  BUCKLING  OF  THE  TOP 
FLANGE  AT  LARGE  DEFORMATION  DUE  TO  BENDING. 
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4. 1.2.3  Teat  Results,  The  results  of  the  test 
on  the  model  frame  are  given  in  Table  4*9  and  the  explana¬ 
tion  of  the  method  of  obtaining  these  results  follows. 

Pig.  4* 1.21  shows  the  horizontal  reaction  variation  with 
half  the  total  vertical  load  obtained  experimentally  from 
the  strain  reading  on  the  steel  bar.  The  experimental 
curve  was  used  to  calculate  the  midspan  moment.  Fig.  4*1*22 
is  a  comparison  of  the  moment  curvature  relation  at  midspan 
obtained  from  the  model  test  and  the  prototype  test.  In 
Pig.  4*1*23  the  load  deflection  results  are  compared  from  the 
model  and  the  prototype  results. 

4. 1.2.4  Conclusions  from  the  Portal  Frame  Test, 

The  results  of  this  test  proved  a  very  fundamental  point 
which  was  not  brought  out  convincingly  by  the  beam  tests, 
namely  that  the  full  plastic  moment  could  be  developed  in 
phosphor  bronze  models  as  it  is  in  structural  steel. 

Pig.  4*1*22  shows  that  the  M- ^curves  of  model  and  prototype 
agree  very  well  considering  the  model  was  at  1:15  scale. 

This  close  agreement  as  compared  to  the  discrepancy  in  the 
simply  supported  beams  is  due  to  the  better  instrumentation 
of  the  frame  where  strains  were  based  on  average  values 
rather  than  a  single  reading. 

The  load  deflection  predicted  curve  is 
in  poorer  agreement  in  this  case  but  still  the  error  is 
less  than  about  10%.  Again  it  should  be  emphasized  that 
the  agreement  in  predicting  ultimate  load  is  a  function  of 
the  maximum  deflection  to  be  tolerated  in  the  prototype. 

Once  this  is  known  then  Pult  can  be  read  off  the  predicted 
curve.  Fig.  4.1*23  emphasizes  this  point  since  at  deflections 
of  the  order  of  magnitude  of  4  to  6  inches  the  predicted 
load  is  in  closer  agreement  to  the  prototype  test. 

Prom  the  test  of  the  model  frame  it  was 
found  that  the  horizontal  reaction  measuring  device  was 
too  stiff  and  its  strain  response  was  low.  This  is  believed 
to  be  the  main  reason  for  the  large  deviation  of  the  experi¬ 
mental  curve  from  the  theoretical  at  relatively  low  loads 
as  shown  in  Pig.  4*1*21. 
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VERTICAL  LOAD  P  IN  LBS. 


FIG. 4. 1.21  VERTICAL  LOAD  vs.  HORIZONTAL  REACTION  FOR 
THE  MODEL  FRAME 
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CENTER  LINE  MOMENT  IN  KIP -INCHES 
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FIG. 4. 1.22  MOMENT- CURVATURE  RELATIONSHIP  FOR 
MODEL  FRAME 
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VERTICAL  LOAD  P  IN  KIPS 


MIDSPAN  DEFLECTION  ,8  ,  IN  INCHES 


FIG. 4.1. 23  LOAD  DEFLECTION  RELATIONSHIP  FOR 
MODEL  FRAME 
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TABLE  ij.,9  (CONT. )  RESULTS  OF  PORTAL  FRAME  TESTS 
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Explanation  of  Table  4.3. 

Col.  1  gives  the  load  in  lbs.  at  the  end  of  the  beam 
of  the  loading  frame. 

Col.  2  gives  the  observed  deflection  of  the  model  frame 
at  midspan. 

Col.  3  shows  the  average  strain  from  two  gage  readings 
on  the  top  flange  at  midspan. 

Col.  4  gives  the  strain  in  on  the  bottom  flange  at 
midspan.  One  gage  reading  only. 

Col.  5  is  the  average  of  two  gage  readings  at  the 
upper  quarter  depth  of  the  web  at  midspan. 

Col.  6  gives  the  average  of  two  readings  of  the  lower 
quarter  of  the  web  at  midspan. 

Col.  7  gives  the  curvature  or  rate  of  change  of  the 
slope  at  midspan.  It  is  calculated  from  the  followiog 
expressions.  (See  Fig.  D.l  ) 

1  „  -SCl  *  2*2  - 

R  4d 

as  derived  in  the  Appendix  (D.2) 

Col.  8  is  the  values  oi  Col.  7  multiplied  by  the  scale 
factor  of  1/15. 

Col.  9  gives  the  total  load,  2P,  acting  on  the  model 
in  lbs.  It  is  obtained  by  multiplying  Col.  1  by  5.58  the 
ratio  of  lever  lengths  for  the  particular  test. 

Col.  10  is  the  predicted  total  load  of  a  prototype 
structure  15  times  larger  than  the  model.  It  is  obtained 
from  the  following  relation: 

pp  ■  Pm  *  V2  1  V1  (G-F-) 
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where  f  is  the  length  scale  and  is  the  ratio  of  stresses 
as  previously  defined.  G.F.  is  a  geometric  factor  to  account 
for  the  difference  in  moment  of  inertia  as  explained  in  Ap¬ 
pendix  C. 

fr‘2  x  fff_1  x  G.P.  =  (15) 2  x  2-1.45  =  652 

Col.  11  is  the  value  of  one  half  of  Col.  10. 

Col.  12  is  Col.  2  times  15  or  f  *  the  length  scale. 

Col.  13  is  the  horizontal  reaction  induced  by  the 
vertical  load  on  the  frame.  It  is  obtained  from  measurements 
of  strain  on  a  stiff  bar.  The  plot  of  P  vs  H  is  given  in 
Fig.  4.1.21. 

Col.  14  gives  the  midspan  moment  of  the  frame  in  lb. -in. 
It  is  obtained  from  the  expression: 

M  =  4-2P  -  5.6H 

as  derived  in  Appendix  F. 

Col.  15  gives  the  predicted  midspan  moment  in  Kip 
inches.  It  is  obtained  from  Col.  14  by  multiplying  with 
the  factor  .652  as  found  for  Col.  10  times  15  because  the 
units  of  moment  are  length  x  force.  The  expression  there¬ 
fore  is 

«p  -  Mm  *  fr'3  1  V1  * 

The  factors  are  the  same  as  stated  previously. 
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4,1.3  Connections.  Welded  steel  structures  have  become 
increasingly  more  popular  because  of  the  continuity  of 
framed  structures  which  ensues  when  joints  and  connections 
are  carefully  welded.  In  the  model  studies  in  phosphor 
bronze  it  was  demonstrated  that  the  full  plastic  moment  of 
the  prototype  could  be  predicted  from  model  studies  st  reduced 
scales.  This  can  occur  in  a  framed  structure  only  if  the 
connections  can  withstand  the  relatively  large  rotations 
Imposed  by  plastic  design.  To  verify  the  adequacy  of  various 
joints  model  connections  were  made  and  tested  in  a  manner 
similar  to  the  tests  on  connections  at  Lehigh  University^^) , 
These  were  not  model  reproductions  of  the  prototype  tests 
however  since  an  8B13  section  was  used  in  the  "prototype". 

This  did  not  allow  correlation  with  the  approximate  phosphor 
bronze  section  used  in  the  model  connections. 

4.1.3. 1  Prototype  Tests.  Since  the  sections  of  the 
prototype  tests  could  not  be  duplicated  in  the  model  only 

the  method  of  testing  was  of  real  interest  in  the  Lehigh 
test  series.  Some  of  the  model  joints  were  made  in  a  similar 
fashion  as  the  prototype.  The  loading  arrangement  shown  in 
Fig.  4.1.24  was  modified  and  used  in  the  model  tests. 

4. 1.3. 2  The  Model  Tests.  The  Joints  which  were 
modeled  are  shown  In  Figs.  4.1.25  and  4.1.26.  The  joint 

in  Fig.  4.1*26  is  a  typical  beam  to  column  connection  whereas 
the  Joints  of  Fig.  4.1.25.  show  the  same  section  in  beam 
and  column. 

The  method  of  making  the  joints  is  illus¬ 
trated  in  Fig.  4*1.27  which  shows  the  various  pieces  which 
are  silver  soldered  together  to  make  connection  K.  Figs. 
4.1.28  -  4*1.30  show  the  model  Joints  after  testing. 

Fig.  4.1.31  shows  the  method  of  loading 
and  support  of  the  model  frames.  The  horizontal  bar  gave 
the  necessary  alignment  so  that  the  load  could  be  applied 
concentrically.  In  Fig.  4.1.32  the  whole  supporting  frame  is 
shown  during  a  test.  Only  deflections  were  measured  in 
these  tests.  From  these  the  rotation  of  the  joint  could  be 
approximated.  Fig.  4.1.33  shows  a  closeup  of  a  te3t  specimen 
in  which "the  method  of  support  Is  clearly  shown. 


PIG.  4.1 
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.24  TEST  SETUP  FOR  CONNECTION  P  (REP.  4*>  ) 
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2.67 


FIG.  4.1.25  DIMENSIONS  OF  THE  MODEL  JOINTS  TESTED 
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PIG.  4.1.27  PARTS  WHICH  GO  INTO  THE  MAKING  OP 
A  MODEL  JOINT. 


FIG.  4.1.28  MODEL  JOINTS  AFTER  FAILURE.  NOTE  THE 
MODE  OF  FAILURE  IN  EACH. 
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PIG.  U.1.29  TYPES  A  AND  C  JOINTS  AFTER  FAILURE.  TYPE  C 
JOINTS  WERE  USED  IN  THE  PORTAL  FRAMES. 


FIG.  4.1.30  TYPE  P  JOINT  AFTER  FAILURE. 
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Load  Application 


FIG.  4.1.31  TYPE  V  JOINT  IN  SUPPORTING  FRAME 
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PIG.  U.1.32  MODEL.  JOINT  IN  THE  SUPPORTING  FRAME  DURING 
TEST. 


FIG.  4.1.33  CLOSEUP  OF  MODEL  JOINT  K-l  AFTER  FAILURE. 

NOTE  THE  BUCKLING  OF  THE  LOWER  BEAM  FLANGE. 
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4.1.3. 3  Test  Results,  Tables  4.10  to  4.17  give 
the  test  results  on  the  8  tests  made  on  model  joints. 

Explanations  of  the  values  in  the  various  columns  of  the 
tables  are  also  given.  The  notations  used  for  these  calcu¬ 
lations  is  given  in  Pig.  4.1.34.  The  moment  vs.  joint 
rotation  curves  are  plotted  in  Pig.  4.1.35  for  all  the 
joints  except  P-1.  All  the  joints  depicted  in  this  figure 
had  the  same  dimensions,  the  same  cross  sectional  area  and 
were  tested  in  the  same  manner.  The  only  variable  was  the 
configuration  of  the  corner.  In  Pig.  4*1.36  the  moment 
rotation  curve  for  Joint  P-1  is  given.  On  the  same  plot 
the  moment  curvature  relation  for  the  phosphor  bronze  sec¬ 
tion  is  given.  Joint  P-1  failed  on  the  column  side  which 
was  apparently  much  weaker.  This  can  be  observed  in  Pig.  4.1.30. 

4.1. 3-4  Conclusions  from  the  Connection  Tests. 

Looking  at  Pig.  4.1.35  a  number  of  conclusions  may  be  drawn. 

1.  Except  for  connection  A-l  all  the  joints  can  mobi¬ 
lize  considerable  joint  rotations.  Joint  A-l  was  the 
unstiffened  Joint. 

2.  The  scatter  from  only  two  tests  of  each  Joint  is 
not  bad  in  the  case  of  joint  M  and  C  but  it  is  bigger  in 
the  case  of  joint  K.  This  is  due  mainly  to  variation  in 
workmanship. 

3.  The  strongest  joints  were  type  M  which  had  the 
most  stiffeners. 

4.  Joints  of  type  C  show  a  lower  value  of  moment 
capacity  than  type  M  but  since  they  have  proved  adequate 
in  the  case  of  the  frame  tests  they  should  be  recommended 
since  they  are  considerably  easier  to  manufacture. 

5.  Type  K  joints  are  considerably  stronger  than  C 
but  these  are  as  hard  to  mske  as  type  M  as  shown  in  Pig. 

4.1.25  so  they  have  no  decisive  advantage  over  type  C. 
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N  AXIAL  COMPONET  OF  P 
V  SHEARING  COMPONENT  OF  P 
SP  VERTICAL  DEFLECTION 
8V  DEFLECTION  ALONG  V 
8n  DEFLECTION  ALONG  N 
e  HALF  THE  JOINT  ROTATION 
d  MOMENT  ARM  TO  JOINT  CENTER 
Mp  JOINT  MOMENT  CAUSED  BY  P 

FIG.  4.1.34  TYPICAL  CONNECTION  BEFORE  AND  AFTER 
BENDING  DEFORMATION  CAUSED 
BY  TWO  COLINEAR  FORCES 
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0.02  0.04  0.06  0.08  0.10  0.12 

JOINT  ROTATION  IN  RADIANS 

FIG.  4.1.35  MOMENT  ROTATION  CURVES  OF  MODEL  JOINTS 
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2 


0.12 


0  0.0 


0.04  0.06  0.08  0.10 

JOINT  ROTATION  IN  RADIANS 


FIG. 4.1. 36  MOMENT  ROTATION  RELATIONS  FOR  JOINT  P'l 


TABLE  4.10  TEST  RESULTS  OP  MODEL  JOINT  K-l 


1 

2 

3 

k 

5 

6 

Load  at 

Vertical 

Total 

p 

Moment 

End  of 

Deflection 

Joint 

Load  on 

at  Joint 

Beam  in 

SD 

Rotation 

Joint 

( Inches) 

( Inches) 

( Inches) 

( Radians) 

(Lbs.) 

(Lbs-In) 

1 

.001*5 

.00318 

. 00146 

5.91 

9.1 

2 

.006 

.001*21* 

.00195 

11.82 

18.2 

3 

.008 

.00565 

. 00259 

17.73 

27.3 

4 

.009 

.00636 

. 00292 

23.64 

36.4 

5 

.Oil 

.00777 

.00354 

29.55 

45.5 

6 

.012 

.  0081*6 

.00388 

35.46 

54.6 

7 

.011*2 

.0100 

. 00459 

41.37 

63.7 

8 

.0165 

.01165 

.00535 

47.28 

72.8 

9 

.0188 

.0133 

.0061 

53.19 

82.0 

10 

.021 

.011*82 

.0068 

59.10 

91.0 

11 

.023 

.01625 

.00746 

65.01 

100.0 

12 

.026 

.01835 

.00832 

70.92 

109.0 

13 

.0285 

.02015 

.0092 

76.83 

118.0 

14 

.0311* 

.0222 

.0102 

82.74 

127.2 

15 

.031* 

.024 

.Oil 

30.65 

136.5 

16 

.038 

.02682 

94.56 

145.8 

17 

.01*6  5 

.03282 

.0151 

100.47 

154.8 

18 

.0535 

.0378 

106.38 

164 

19 

.072 

.0508 

.0234 

112.29 

3.73 

20 

.127 

.090 

.0414 

118.20 

182 

21 

.25 

.177 

.0811 

324.11 

191 

193 


Explanation  of  Table  4.10 

Col.  1  Is  the  load  at  the  end  of  the  beam  of  the 
loading  frame.  This  load  was  applied  In  increments  of 
one  pound  and  In  such  a  manner  that  enough  time  was  given 
to  the  specimen  to  allow  for  creep  beyond  the  yield  point. 

Col.  2  is  the  deflection  in  the  direction  of  the 
applied  load  P  on  the  joint.  This  was  measured  with  an 
Ames  dial  gage  which  reads  to  the  nearest  one  thousandth 
of  an  inch.  Since  the  bottom  leg  of  the  joint  was 
prevented  from  moving  vertically  at  Point  B  (Pig.  4.1.34) 
the  upper  leg  at  the  point  of  load  application  A  (Pig.  4*1. 34) 
moved  down  2  x  p. 

Col.  3  is  the  component  of  the  vertical  deflection 
of  Col.  2  in  the  direction  of  the  shearing  component  of 
force  V. 

* 

Col.  4  is  the  total  joint  rotation  or  2  9  as  shown 
on  the  Pig.  3.1.34.  This  was  obtained  approximately  by 
dividing  the  deflection^  by  the  moment  arm  a. 

Col.  5  is  the  load  P  on  the  joint  which  is  obtained 
from  the  load  at  the  end  of  the  beam  by  multiplication  with 
the  factor  of  the  loading  ‘frame  depending  on  the  position 
of  the  specimen  along  the  beam. 

Col.  6  is  the  moment  at  the  knee  or  the  intersection 
of  the  2  center  lines  of  the  members  forming  the  joint.  It 
is  obtained  by  multi plying  the  shearing  component  of  P  by 
the  moment  arm  a.  Note  that  the  axial  force  N  has  no  moment 
component  at  the  knee. 

M  =  (.707P)  (2.18)  =  1.54P  (in-lb.) 


TABLE  4.11  TEST  RESULTS  OP  MODEL  JOINT  K-2 

_ , _  -  - —  -  --- » 


1 

2 

* 

4 

5 

6 

Load  at 

Vertical 

* 

Total 

p 

Moment 

End  of 

Deflection 

0  V 

Joint 

Load  on 

at  Joint 

Beam  In 

op 

Rotation 

Joint 

(Lbs. ) 

( Inches) 

( Inches) : 

(Radians) 

(Lbs.) 

(Lbs-In) 

1 

.0018 

.00127  I 

.000582 

5.91 

9.1 

2 

.004 

.00282  | 

.001295 

11.82 

18.2 

3 

.007 

. 00494  1 

.00227 

17.73 

27.3 

4 

.009 

.00635 

.00292 

23.64 

36.4 

5 

.0112 

.0079 

.00362 

29.55 

45.5 

6 

.0138 

.00975 

.00448 

35.46 

54.6 

7 

.0174 

.0123 

.00565 

41.37 

63.7 

l  8 

.£-94 

.0137 

.00629 

47.28 

72.8 

9 

.0218 

.0154 

.00707 

53.19 

82.0 

10 

.024 

.017 

.0078 

59.10 

91.0 

11 

.028 

.0198 

.0091 

65;  01 

100.0 

12 

.032 

.02262 

.0104 

70.92 

109.0 

13 

.036 

. 02545 

.0117 

76.83 

118.0 

14 

.0406 

.0287 

.0132 

82.74 

127.2 

15 

.050 

.0353 

.01622 

88.65 

136.5  i 

16 

.0595 

.0421 

.01935 

94.56 

145.8  I 

17 

.0735 

.052 

.0239 

100.47 

154.8  ! 

18 

.111 

.0785 

.0361 

106.38 

164  { 

19 

.258 

.1825 

.0836 

- 1 

112.29 

■”  j 

195 


TABLE  4.12  TEST  RESULTS  OP  MODEL  JOINT  M-l 


1 

Load  at 
End  of 
Beam  in 
(Lbs. ) 

2 

Vertical 

Deflection 

op 

( Inches) 

3 

Sv 

( Inches ) 

k 

Total 

Joint 

Rotation 

(Radians) 

5 

p 

Load  on 
Joint 
(Lbs. ) 

6 

Moment 
at  Joint 

(Lbs. -In) 

1 

.001; 

.00283 

.0013 

5.91 

9.1 

2 

.0075 

.0053 

.00243 

11.82 

18.2 

3 

.010 

.00707 

.00324 

17.73 

27.3 

4 

.011 

.00777 

.00356 

23.64 

36.4 

5 

.Oik 

.0099 

.00454 

29.55 

45.5 

6 

.0165 

. 01165 

.00535 

35.46 

54.6 

7 

.019 

.0134 

.00615 

41.37 

63.7 

8 

.0215 

.0152 

.00698 

47.28 

72.8 

9 

.0235 

.0166 

.00762 

53.19 

81.9 

10 

.0265 

.0187 

.00859 

59.10 

91.0 

11 

.029 

.0205 

.0094 

65.01 

100.0 

12 

.0315 

.0222 

.0102 

70.92 

109.0 

13 

.03k5 

.024k 

.0112 

76.83  ■ 

118.0 

14 

.0375 

.0265 

.01215 

87.74 

127.2 

15 

.oki 

.029 

.0133 

88.65 

136.5 

16 

.0k7 

.0332 

.01522 

94.56 

145.8 

17 

.060 

.0424 

.01945 

100.47 

154.8 

18 

.  07k5 

.0526 

.02415 

106.38 

164 

19 

.109 

.077 

.0354 

112.29 

173 

20 

.ik5 

.1025 

.047 

118.20 

182 

21 

.227 

.161 

.074 

124.11 

191 

22 

.390 

.276 

.1265 

130.02 

200.5 

196 


TABLE  4.1 3  TEST  RESULTS  OP  MODEL  JOINT  M-2 


1 

2 

3 

Load  at 

Vertical 

0 

End  of 

dv 

Beam  in 
(Lbs. ) 

Deflection 
( Inches) 

( Inches) 

1 

.0025 

.00177 

2 

.0035 

.0021+8 

3 

.005 

.00354 

4 

.007 

. 00495 

5 

.0085 

.0060 

6 

.0100 

.00707 

7 

.0115 

.00815 

8 

.0128 

.00905 

9 

.0135 

.00955 

10 

.011+5 

.01025 

11 

.0159 

.01125 

12 

.017 

.012 

13 

.018 

.01275 

14 

.019 

.01345 

15 

.0205 

.0145 

16 

.022 

.01555 

17 

.023 

.01625 

18 

.0245 

.01732 

19 

.0258 

.01825 

20 

.027 

.0191 

21 

.0285 

.0202 

22 

.030 

.0212 

23 

.032 

.02263 

2k 

.0336 

.0238 

25 

.0352 

.0249 

26 

.0375 

.0265 

27 

.0395 

.0280 

28 

.0424 

.03 

29 

.045 

.0318 

30 

.050 

.0354 

31 

.054 

.0382 

32 

.060 

.0424 
.  046o 

33 

.066 

3k 

.0775 

.0548 

35 

.087 

.0615 

36 

.100 

:X7 

37 

.1215 

38 

.146 

.1032 

39 

.280 

.198 

TABLE  4.14  VEST  RESULTS  OF  h  ■)LEL  JOTKT  c-1 


— 

1 

’  "1 

2  ; 

-  i 

‘> 

2 

j  5 

6 

Load  at 

Vertical 

Total  1 

li  p 

Moment 

End  of 

Deflection  1 

t  v 

Joint 

%\Load  on 

at  Joint 

Beam  in 

Sp  1 

Rotation 

W.'i  Joint 

(Lbs.) 

(Inches) 

( It  chea) 

( Radians ) 

iv  LbJ’l 

(Lbs-In) 

1 

.001 

.000707 

.000324 

V  3 . 22 

4.96 

2 

.002 

..0,414 

.00065 

%  *44 

9.82 

3 

.0035 

.  0348 

.00114 

,  66 

14.9 

4 

.005 

.0.:53 

. 00162 

11%  88 

19.83 

5 

.0065 

.  OCli  .6 

.00211 

16.\L0 

24.80 

6 

.0082 

.  O'  ‘  8 

.00266 

19. 

29.80 

7 

.01 

.Ou  07 

.00324 

22. 54\ 

1  34.70 

8 

.0115 

. 00.12 

.00372 

25. 76  \ 

39.70 

9 

.0128 

.00905 

.00415 

28.98  1 

54.6o 

10 

.  0148 

.01048 

.00481 

32.6 

X  43.60 

11 

.0162 

.31  5 

. 00528 

35.42 

X*  ,60 

12 

.018 

.01272 

. 00584 

38.64 

59  Xp 

13 

.  020 

.01.04 

.0065 

M.  86 

64.5^  4 

14 

.0228 

.0165 

.00747 

45.08 

69.5  1 

15 

.025 

.01  vi 

.00812 

48.30 

74.5 

16 

.  028 

.  Oj  98 

. 00907 

51.52 

79.5 

17 

.on 

07195 

.0102 

64.74 

84.4 

18 

.O: 

.  0248 

.0114 

57.96 

89.2 

19 

.  on:-' 

. C2725 

.0125 

61.18 

93.4 

20 

.041+5 

.0315 

.01445 

64.40 

99.2 

21 

.0525 

.0372 

.0171 

67.62 

104.2  ; 

22 

.0615 

.  0:4.35 

.  0200 

70.84 

109 

23 

.072 

.on 

.0234 

74.06 

114.2 

24 

j  . 0865 

.0612 

.0281 

77.28 

119 

25 

.102 

.  0’722 

.031 

80.50 

124 

26 

.1365 

.0965 

.0433 

83.72 

129 

27 

250 

.163 

.075 

86.94 

133.8 

TABLE  4.1 5  TEST  RESULTS  OP  MODEL  JOINT  C-2 


l 

2 

1 

k 

1 

6 

•mm 

Load  at 

Vertical 

Total 

P 

Moment 

End  of 

Deflection 

dv 

Joint 

Load  on 

at  Joint 

Bean  in 

op 

Rotation 

Joint 

(Lbs. ) 

(Inches) 

( Inches) 

(Radians) 

(Lbs.) 

(Lbs-In) 

3 

.001 

.000707 

.000325 

3.22 

4.96 

C. 

.0015 

. 00106 

. 000486 

6.44 

9.82 

“5 

.0022 

.00155 

.000712 

9.66 

14.9 

4 

.003 

.00212 

.00974 

12.88 

19.83 

5 

.oolj. 

.00282 

.00129 

16.10 

24.80 

6 

.0052 

.00367 

.001685 

19.32 

29.80 

I 

.0065 

.00459 

.00211 

22.54 

34.70 

8 

.008 

. 00565 

. 00259 

25.76 

39.70 

9 

.01 

.00707 

.0032 

28.98 

44.60 

10 

.0115 

.OO815 

.00374 

32.6 

49.60 

11 

.0135 

. 00955 

. 00438 

35.42 

54.60 

12 

.0155 

.011 

.00505 

38.64 

59.60 

13 

.0178 

.0126 

.00579 

41.86 

64.5 

14 

.020 

.01414 

.0065 

45.08 

69.5 

15 

.0242 

.0171 

.00785 

48.30 

74.5 

16 

.025 

.0177 

.00813 

51.52 

79.5 

1? 

.0266 

.0188 

.00862 

54.74 

84.4 

18 

.0292 

.02062 

.00946 

57.96 

89.2 

19 

.032 

.0226 

.0104 

61.18 

93.4 

20 

.0356 

.0252 

.0115 

64*4° 

99.2  I 

21 

.042 

.0297 

.0136 

67.62 

104.2 

22 

.0502 

.0355 

.0165 

70.84 

109 

H 

.0618 

.075 

.0438 

.053 

.201 

.02k3 

74.06 

77.28 

114.2 

119 

25 

.0872 

.0616 

.0283 

80.50 

124 

26 

.106 

.075 

.0344 

83/72 

129 

.127 

.0898 

.r4l2 

86.94 

133.8 

28 

.151 

.1068 

.049 

90.16 

141 

29 

.1635 

.13 

.0596 

93.38 

143.6 

30 

.2585 

.183 

.084 

96.60 

148.8 

31 

.300 

.212 

.0972 

99.82 

153.5 

199 


TABLE  i|..l6  TEST  RESULTS  OP  MODEL  JOINT  A-l 


1 

2 

2 

— 

4 

1 

6 

Load  at 

Vortical 

Total 

P 

Moment 

End  of 

Deflection 

by 

Joint 

Load  on 

at  Joint 

Beam  in 

Sp 

Rotation 

Joint 

(Lbs.) 

(Inches) 

(Inches) 

(Radians) 

(Lbs.) 

(Lbs-In) 

1 

.002 

.001414 

.00065 

3.22 

4.96 

2 

.0055 

.0039 

.00179 

6.44 

9.82 

3 

.0076 

.00536 

.00246 

9.66 

1.4.9 

4 

.0102 

.0072 

.0033 

12.88 

19.83 

5 

.0128 

.00905 

.00416 

16.10 

24.30 

6 

.'15 

.0106 

.0046 

19.32 

29.80 

7 

.017 

.012 

.0055 

22.54 

34.70 

8 

.19 

.0134 

.00615 

25.76 

39.70 

9 

.021 

.GI482 

.0068 

28.98 

44.70 

.0232 

.0164 

.00753 

32.6 

49.60 

.026 

.0184 

.00845 

35.42 

54.60 

.0285 

.0201 

.00923 

38.o4 

59.60 

.0316 

.0224 

.0103 

41.86 

64.5 

.035 

.0247 

.0113 

45.08 

69.5 

.046 

.0325 

.0149 

48.30 

74.5 

.054 

.0382 

.0175 

51.52 

79.5 

.066 

.0466 

.0214 

54.74 

84.4 

,082 

.058 

.0266 

57.96 

89.2 

19 

.110 

.0777 

.0357 

61.18 

93.7 
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TABLE  4.17  TEST  RESULTS  OP  MODEL  JOINT  P-1 


Load  at 
End  of 
Beam  In 
(Lbs. ) 


2 

Vertical 

Deflection 

op 

{ Inches) 


.0065 

.009 

.0105 

.013 

.0155 

.018 

.0202 

.0232 

.027 

.0315 

.0365 

.0425 

.050 

.0715 

.136 

.239 

.430 


1 

( Inches) 


.0046 
. 00635 
. 00742 
.0092 
.01095 
.01270 
. 01425 
.0164 
.0191 
.0222 
.0258 
.03 
.0353 
.0505 
.096 
.169 
.304 


k 

1 

6 

Total 

P 

Moment 

Joint 

Load  on 

at  Joint 

Rotation 

Joint 

(Radians) 

(Lbs.) 

(Lbs-In) 

.00211 

3.22 

11.7 

.00291 
. OO341 

23.4 

35.2 

.00422 

12.88 

46.8 

.00502 

16.10 

58.5 

. 00582 

19.32 

70.2 

.00655 

22.54 

02 

.00755 

25.76 

93.6 

.00875 

28.98 

105.2 

.0102 

32.6 

118.5 

.01182 

35.42 

129.0 

.01375 

.0162 

38.6k 

41.86 

140.8 

152.2 

.02318 

45.08 

164 

.044 

48.30 

175.8 

.0775 

51.52 

187.2 

.1395  | 

54.74 

199 
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14..  2  DYNAMIC  TESTS  ON  MODEL  PORTAL  FRAMES  MADE  OF  PHOSPHOR 
BRONZE 

Tests  on  fixed  ended  portal  frames  were  performed 
to  check  the  dynamic  response  of  simple  structures  and  com¬ 
pare  the. Increase  in  strength  to  that  of  a  static  test. 

In  these  dynamic  tests  the  impulsive  loading  machine  used 
is  described  in  Section*  3. 4* 2.  The  static  test  was  performed 
in  the  beam  loading  apparatus. 

4*2.1  Type  of  Pulse  Possible.  The  impulsive  loading 
machine  has  the  capability  of  giving  rise  times  down  to 
3-4  milliseconds.  However,  these  are  more  impactive  rather 
than  impulsive  type  loads  since  to  achieve  these  small  rise 
times  the  ram  must  not  be  in  contact  with  the  structure  but 
must  have  an  initial  travel.  A  typical  load  pulse  experienced 
by  a  frame  is  shown  in  Fig.  4*2.1  where  a  rise  time  of  about 
3  milliseconds  was  achieved. 

4*2.2  The  Model.  The  portal  frame  was  made  of  the 
same  phosphor  bronze  section  which  approximates  the  8WF40 
rolled  section  at  1:15  scale.  The  dimensions  and  the  details 
of  obtaining  the  necessary  fixity  are  shown  in  Fig.  4*2.2. 

The  model  was  supported  in  a  special  frame  shown  in  Fig.  3.4*11 
which  was  placed  in  the  machine  in  such  a  way  so  that  the 
vertical  travel  of  the  machine  ram  produced  the  lateral  load 
on  the  model.  Load  was  measured  by  the  load  cell  described 
earlier  and  the  lateral  displacement  by  an  L.V.D. T.  Both 
of  these  outputs  were  displayed  on  a  cithode  ray  oscilloscope 
and  a  permanent  record  was  taken  photographically.  Fig.  4*2*3 
shows  the  model  frame  before  and  after  the  test. 

4*2.3  The  Test  Results.  The  load  time  variation  as 
experienced  by  the  model  frame  is  plotted  from  the  photo¬ 
graphic  record  in  Fig.  4*2.4*  Fig.  4*2.5  gives  the  displace  - 
meat  time  variation  for  the  same  frame.  From  these  two  plots 
we  can  find  the  load  displacement  curve  of  the  frame  as 
shown  in  Fig.  4*2.6.  This  is  the  dynamic  response  of  the 
model  frame  covering  the  early  part  of  the  response. 
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1 


LOAD  IN  LBS 


FIG. 4.2. 1  TYPICAL  TIME  LOAD  PULSE  EXPERIENCED 
BY  FRAME 
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a.  DIMENSIONS  OF  MODEL  PHOSPHOR  BRONZE  FRAME 
AT  CIS  SCALE 


d.  DETAIL  AT  BASE  OF  COLUMN  TO  INSURE  FIXITY 


FIG.4.2.2DRAWING  OF  1:15  SCALE  MODEL  SHOWING  CONSTRUCTION  DETAILS 
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a. 


BEAM  AND  COLUMN  PARTS  BEFORE  ASSEMBLY 


b.  WELDED  FRAME  BEFORE  TESTING 


c.  FRAME  #1,  AFTER  TESTING 

PIC.  4.2.3  PICTURES  SHOWING  THE  MODEL  BEFORE  AND  AFTER 
THE  DYNAMIC  TEST 
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DYNAMIC  HORIZONTAL  LOAD  P(t)  IN  LBS. 


FIG.  4.2.4  LOAD  TIME  RELATION  FOR  MODEL  FRAME 
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DYNAMIC  DISPLACEMENT  (x)  IN  INCHES 


FIG.  4.2.5  DISPLACEMENT  TIME  RELATION  FOR  MODEL 
FRAME 
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HORIZONTAL  LOAD  (LBS) 


FIG.  4.2.6  LOAD  DEFLECTION  CURVE  FOR  MODEL 
FRAME,  DYNAMIC  TEST 
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7  A  static  tast  was  also  made  on  a  similar  frame 

and  the  results  are  given  in  Table  l(..l8.  The  load  deflection 
curve  for  this  model  is  given  in  Fig.  it-.  2.7. 

k»2*k  Conclusions  from  the  Dynamic  Tests.  Ho  correla¬ 
tion  between  model  and  prototype  could  be  obtained  in  this 
case  sinoe  there  were  no  known  full  scale  tests  to  model. 
However,  some  conclusions  could  be  drawn  from  this  limited 
study  of  dynamic  behavior. 

1.  It  can  be  seen  that  the  dynamic  load  taken  by  the 
model  frame  is  considerably  larger  than  the  static  test 
value.  The  dynamic  Pult.  is  about  220  lbs.  and  the  static 
is  about  i30  lbs.  where  equal  deflections  were  used  as  a 
criterion  of  finding  the  static  Pult. 

2.  The  mode  of  failure  in  both  static  and  dynamic  tests 
was  essentially  the  same. 
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HORIZONTA 


HORIZONTAL  DEFLECTION  IN  INCHES 

FIG.  4.2.7  LOAD  DEFLECTION  CURVE  FOR  MODEL 
FRAME,  STATIC  TEST 
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APPENDIX  A 


Consider  a  linearly  elastic  beam  of  length  L, 
stiffness  El,  coefficient  of  damping  4,  and  mass  3  per 
unit  length.  The  beam  is  subjected  to  a  solid  body 
acceleration  i  by  the  application  of  constant  forces 

p  =  1  nil 

t  *  to  its  two  ends.  It  is  desired  to  deter¬ 

mine  the  deviation  at  time  t  of  the  bending  mopent  of  points 
along  the  beam  from  their  steady  state  value 


t  F*^iL 

MAM  {/UNIT  LKNSTM 

ctiwnbm  ex 

•  ■  u  ™ 

tL. 

Pig.  A, 1 

Let  the  x  -  and  y-  axes  of  Figure  A,1  move  with  the  beam. 
Then  the  governing  equation  for  dynamic  equilibrium  of  the 
beam  is 


b  y 

FT  " 

=  r* 

b  y 

A  . 

b  y 

(Al) 

5  b  t* 

□  ■ 

bt 

Initial  conditions: 

y  * 

^  y 
dt 

• 

H 

O 

for 

t  = 

0 

(A  2) 

Bound,  conditions: 

y  * 

V  y 
fcx* 

-  =  0 

for 

and 

X  = 

X  = 

0 

L 

(A3) 
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Since  the  governing  equation  la  not  homogeneous * 


substitute  y(x,t)  = 

Y(x,t)  + 

f  (x) 

and  Eqn. 

(Al)  becomes 

El 

fVY  + 

1 

•• 

-  ti  > 

* 

L  Bx4 

**4  J 

Equation 

(A5)  is  homogeneous  if 

8*  ♦(«) 

=  $ 

dx4 

EI 

El  =  . 

B*Y 

_  -  d 

By 

l.e. 

2x4 

8  fct 

e 

Bt 

-  cl 


(A4) 

W  U5) 


Bt 

2  ( A6) 

(A7) 


The  boundary  conditions  of  the  new  governing  eqn.  (A7) 
are  homogeneous  if,  according  to  (A3)  and  (AI4.) 


Bx« 


=  0  for  x=0  and  x=I, 


( A8) 


Equation  (A6)  with  boundary  conditions  (A8)  represents  a 
uniformly  loaded  simply  supported  beam  so  that 


±L 

2kEI 


[ 


__  3  a 

x4  ~  2Lx  +  L  x 


] 


(A9) 


It  remains  to  solve  (A7)  subject  to  the  initial  conditions 

Y  <x,o>  •-♦(*)  (Ai0) 

«ru'0)  =0 

and  boundary  conditions 

Y  (o,t)  =  Y(L,t)  *  0 

£*-< o,t).sq_  (l>t)  =  0 

A  trial  solution  Y  =  f(x)  g(t)  leads  to 

f""_  e?-  «"  +  -ir  e' 


(All) 


•■•Y- 


8 


^  c^ '  cos  wt  +  c'  sin  wtj  [cx 

+  c-  cosh  kx  +  c  ainK  kxl 
a  4  J 


-  *  k4  where  k  is  real,  an 

f  is  differentiated  with  resp 
to  x,  g  with  respect  to  t. 


cos  kx  +  c  sin  kx 
9 
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According  to  the  boundary  conditions  (All) 


c  =  c  =  c  =  0 

18  4 


and  sin  kL  =  0 

.  • .  k  =  J23L. 

L 


also  w  = 


Generally  thus 
Y 


VwVe 

Y  UlUB  j 

-  [n 


El  d« 

4$* 


n  =  0,  ly  2j . .. 

(A12) 


coswt  +  c  ainwt 
2 


J  sin 


(A13) 

L 


ntx 


Prom  the  second  of  initial  conditions  (A10): 

c  =  — c 
•  2$w  1 

so  that  (A13)  becomes 

Y  =  cn  •  ^  Jjjoswt  sin  wtj  sin 


The  first  initial  condition  yields 

M  • 

YL  cn  *ln  ~  =  -  9{x)  =  24ST  [  2L*S‘  L*x] 

n=l 

The  Peurler  coefficients,  determined  in  the  normal 


fashion,  are 


Cn  = 
n  El 


The  general  solution  for  the  deflection  of  the  beam  is 
thus 


is  [i>  i**  ai*y  [(-if-  i]  "!•* 

—  •  T -  *  2T+P  —  * 


(cos  wt  +  sin  wt)  sin  n  — 

ZV  l 
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The  first  three  terms  Inside  the  bracket  represent  the 
static  deflection  of  the  beam,  and  the  last  term  is  an 
oscillation  the  amplitude  of  which  is  decreasing  with 
time. 

The  moment  in  the  beam  at  time  t  is 


U 


It  is  desired  to  give  the  variation  of  M  with  time  for 
various  percentages  p  of  critical  damping.  Prom  (A12) 


dcr  ~ 

2  n*v* 

L8 

V 

SEI  ' 

U 

w  = 

n«  ir* 

■  V 

-s  (1  -  p) 1 

u 

L8 

V 

% 

d 

-  M 

n8  T8 

/IT 

(A 

L8 

Substitution  of  (A16)  and  (A17)  into  (Allj.)  yields  all  the 
desired  results. 
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Table  B.l 


Appendix  B:  Computations  for  Section  2,7*1 


X 

l 

b  =  0.1 fc 

y  * 

b=0. 25  ft 

D 

b=0.50  ft 

b=0.1Cft 

y  =  t 

b=0. 25  ft 

b=0.50ft 

0 

-0.0041 

-0.0093 

-0.0128 

-1.128 

-0.753 

-0.559 

0.125 

-0.0038 

-0.0085 

-0.0118 

-0.447 

-0.704 

-0.642 

0.250 

-0.0028 

-0.0064 

-0.0090 

0.119 

-0.506 

-0.721 

0.375 

-0.0014 

-0.0034 

-0.0049 

0.228 

-0.106 

-0.516 

0.500 

0.0001 

0.0002 

0 

0.050 

0.361 

0 

0.625 

0.0016 

0.0036 

0.0049 

0.070 

0.616 

0.516 

o.75o 

0.0028 

0.0064 

0. 0090 

0.201 

0.509 

0.721 

0.875 

0.0036 

0.0083 

0.0118 

0.149 

0.149 

0.642 

1.000 

0.0039 

0.0089 

0.0128 

O.O63 

0.031 

0.559 

Values  of  the  ratio 


f  x 

P 


for  t  =  l. 
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Table  B.2 


X 

l 

*x_ 

P 

a. 

p 

p 

>.8  t, 

* 

b  = 
tmax 

P 

0.1  t 

£1 

P 

ff, 

P 

0 

-0.1851 

0.2490 

0 

90 

0.317 

0.249 

-0.185 

0.125 

-0.1342 

0.1990 

0.0968 

105 

0.193 

0.225 

-0.160 

0.250 

-0.0388 

0.1113 

0.1171 

119 

0.139 

0.180 

-0.098 

0.375 

0.0187 

0.0624 

0.0839 

12? 

0.087 

0.127 

-0.046 

0.500 

0.0362 

0.051U 

0.0656 

132 

0.066 

0.110 

-0.022 

0.625 

0. 0531 

0.0164 

0.0535 

145 

0.057 

0.091 

-0.021 

0.750 

0.0682 

0.0189 

0.0317 

154 

0.040 

0.084 

0.003 

0.875 

0.0642 

0.0213 

0.0072 

170 

0.022 

0.064 

0.020 

y  =  c 

).8  t. 

b  =  ( 

5.25  * 

X 

*x 

W 

Tmax 

** 

l 

P 

p 

IK. 

P 

P 

P 

0 

-0.291 

0.456 

0 

90 

0.373 

0.456 

Bsa 

0.125 

-0.259 

0.429 

0.198 

105 

0.397 

0.482 

0.250 

-0.169 

0.353 

0.195 

108 

0.325 

0.417 

MU 

0.375 

-0.039 

0.240 

0.234. 

120 

0.273 

0.373 

-0.173 

0.500 

0.087 

0.125 

0.208 

132 

0.209 

0.315 

-0.102 

0.625 

0.162 

0.055 

0.131 

146 

0.142 

0.250 

-0.033 

0.750 

0.169 

0.047 

0.048 

0.067 

0.175 

0.041 

0.875 

0.136 

0.076 

0.005 

0.030 

0.136 

0.076 

1.000 

0.056 

0.094 

0 

1 

C.019 

0.094 

0.056 
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.2896 
.2802 
0.2389 
0.11+31+  0.5316 
0  0.1+000 
0.11+31+  0.2681+ 
0.2389  0.1821 
0.2802  0.11+98 
0.2896  0.11+36 


y  =  0.8  t,  b  =  0.50  l 

1m  «  — 

p  p  p 


0 

90 

iiPiTCn 

0.0786 

95 

0.1711 

101 

0.2577 

109 

0.1+21+ 

0.291+1+ 

118 

0.355 

0.2577 

128 

0.265 

0.17H 

138 

0.191+ 

0.0786 

155 

0.101+ 

0 

. _ 

180 

0.073 

-0.230 

-0.155 

-0.059 

0.017 

0.111 

0.11+1+ 


0 

.125 
.250 
.375 
0.500 
0.625 
0.750 
0.875 
1.000 


0.031+0 

0.0083 

0.0117 

0.021+7 

0.0335 

0.0395 

0.01+18 


90 

0.0221+  98 
0.0361  105 
0.0382  113 
0.0327  120 
0.021+8  126 
0.0171  131+ 
0.0090  11+6 
0 


0.053 


0.038  0.069 

0.026  0.059 
0.017  0.052 
0.010  0.01+6 


0.012  0.01+2 


-0.01+1+ 

-0.025 

-0.007 

0.007 

0.017 

0.026 

0.030 


Table  B.4 


b  = 
*  max 


P 


0.1334 
0.1191 
0.0828 1  0.2116 


0 

.125 

.250 

.375 

0.500 

0.625 

0.750 

0.875 

1.000 


.1593 

.i486 

■0.1176 

■0.0653 

0 

0.0653 

0.1176 

0.1486 

0.1593 


90 

96 

0.0708 

103 

0.0859 

110 

0.0822 

118 

0.0628 

127 

0.0387 

133 

0.0172 

151 

0.249 

0.244 

0.229 

0.203 

0.177 

0.149 

0.125 

0.107 


.4169  0  90 

.4068  0.0403  94 
.3867  0.0778  98 

.3485  0.1062  104 
0.3000  0.1171  109 

0.2515  0.1062  114 
0.2133  0.0778  119 

0.1932  0.0403  121 
0.1832  0  90 


0.264 

0.233 

0.190 

0.140 

0.091 

0.046 

0.012 


.133 

.124 

.099 

.063 

.023 

.016 

.047 


.417 
.409 
.398 
.375 
0.340 
0.298 
0.257 
0.217 
0.183 


-0.129 

-0.091 

-0.040 

0.018 

0.074 

0.125 

0.159 


0.006? 

0.0117 
0.0141  105  0.029 
0.0139  109  0.022 
0.0115  113  0.016 
0.0079  Ilk  0.011 


0.029  0.01*8 

0.022  0.01*3 
0.016  0.038 
0.011  0.033 


0.001*1  111  0.006  0.029 

0  90  0.001*  0.027 


-0.018 

-0.010 

-0.002 

0.005 

0.012 

0.017 

0.019 


.125 
.250 
.375 
0.500 
0.625 
0.750 
0.875 
1.000 


371* 

338  0. 

21*0  0. 

0.0222  0.1106  0. 

0.001*1*  0.0966  0.0315 

0.0222  0.081*6  0.0267 
0.0352  0.0760  0.0188 
0.01*30  0.0710  0.0097 

0.01*56  0.0691*  0 


90 
95 
99 

103  0.073 

107  0.056  0.106  -0.005 

110  0.01*1  0.094  0.012 

112  0.028  0.083  0.028 

108  0.017  0.072  0.040 

90  0.012  0.069  0.046 
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Table  B.6 


X 

l 

ffx 

p 

lx 
*  j 

y  = 

> 

0.4 

oC 

t  b 

r 

max 

P 

=  0.25 

Ti 

«• 

F“ 

-0.0708 

0.2482 

0 

0.1595 

0.248 

-0.071 

-0. 06 52 

0.2442 

0.0170 

93 

0.1550 

0.245 

-0.066 

0.250 

-0.0498 

0.2339 

0.0314 

96 

0.1450 

0.237 

-0.053 

0.375 

-0.0265 

0.2183 

0.0411 

99 

0.1295 

0.225 

-0.034 

0.500 

0 

0,2000 

0. 0466 

102 

0.1100 

0.210 

-0.010 

0.625 

0. 0265 

0.1818 

0.0411 

104 

0.0875 

•0.192 

0.017 

0.750 

0.0498 

0.1661 

0.0314 

0.0659 

0.174 

0.042 

0.375 

0.0652 

0.1557 

0.0170 

0.0481 

0.154 

0.057 

1.000 

0.0708 

0.1518 

0 

0.0405 

0.152 

0.071 

APPENDIX  C 


As  discussed  earlier,  the  cross  section  chosen 
for  the  phosphor  bronze  models  was  only  an  approximation 
of  the  prototype  at  the  length  scale  used.  This  was  done 
in  order  to  faciliatate  the  machining  operations  and  to 
cut  down  on  the  cost.  It  was  also  assumed  that  geometric 
corrections  could  be  easily  applied  to  the  model  results. 

The  most  significant  geometric  property  involved 
in  the  type  of  tests  performed  was  the  moment  of  inertia. 

This  was  the  case  because  bending  moments  were  of  primary 
concern.  The  effects  of  axial  forces  and  shearing  distor¬ 
tions  were  neglected. 

We  define  the  geometric  factor,  6.P.,  as  follows: 

(j  p  -  Geometric  property  of  Exact  Model 

Geometric  property  of  Actual  Model  (C.l) 


For  cases  Involving  the  moment  of  inertia  of  the 
section  the  following  correction  should  be  used  to  the 
values  predicted  from  the  model  test: 


G.F. 


1  x  t« 
P  ** 


m 


(C.2) 


where, 

I  =  moment  of  inertia  of  prototype  about  the 
p  strong  axis. 

I  =  moment  of  Inertia  of  model  section  about 
strong  axis. 

f  =  length  scale. 

In  the  case  of  the  8WFI4.O  section  used  as  a  prototype  at  a 
scale  of  1:15  the  model  section  had  an  I  =  .002  in^-  and 

Dl 

thus  the  G.F.  is 

G.F.  =  • 2-  *■  ( =  I.45 

.002 
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This  was  used  in  all  computations  of  moment  and  load  predic¬ 
tion  from  the  phosphor  bronze  model  tests. 
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APPENDIX  D 


In  a  section  under  constant  nonent  we  can  find 
the  curvature  or  the  change  In  slope  of  the  elastic  curve 
per  unit  length  from  the  measurement  of  strains  across 
the  critical  section. 

Prom  elementary  beam  theory  we  know  that  using 
the  convention  of  Fig.  D.l 

1  =_M_  (D.l) 

R  El 

in  the  elastic  range  and  where  l/R  is  usually  taken  as 

a 

dx  .  If  we  measure  the  strains  across  the  critical  section 
and  assume  that  plane  sections  before  bending  remain  plane 
after  bending  we  get: 

i  .  i  .  ii  (d.2) 

R  tor  dx 

where  ^  is  the  slope  to  the  elastic  curve  or  the  rotation 
of  the  critical  section. 

Assuming  that  strain  measurements  are  taken  at 
the  top  and  bottom  flanges  and  at  the  upper  and  lower 
quarter  points  of  the  web  we  have  looking  at  Fig.  D.l 

*  top  flange  strain  (-ve) 

€  2  s  upper  quarter  point  strain  in  web 

€:  ^  “  lower  quarter  point  strain  in  web 

€  *  bottom  flange  strain  (+ve) 

d  *  depth  of  section 

For  an  infinitesimal  length  of  beam  at  the  critical  section, 
say  dx,  we  have, 


22k 


it 

dx 


.  *3  'S' 

d/2 

.€jl  =h 


(D.3) 


dA 


dA 


The  average  of  these  values  is  taken  as  the  curvature: 

1  =  .  £h  ■  £1  +  2(63  '  €2}  +  ^(62  ~  gl* 

R  d*  4d  i|.d  l|d 

-  63)  (DA) 


1 

R 


Ud 

-  2€2  1  2g3  -  *\ 
4d 
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APPENDIX  E 


If  deflections  ara  measured  along  a  section  of 
a  beam  in  which  the  moment  is  constant  or  approximately  so, 
then  the  curvature  could  be  determined.  Using  the  notation 
of  Fig.  E. 1  ,  we  can  derive  an  expression  for  the 
second  derivative  of  the  deflection  at  point  2  as  follows. 


1 

R 


(E.l) 


Slopes  -  Slope-,  o 

(Second  Slope)  ? 

Ax 

where 


S1°P®l-2  = 

y2  -  yl 

1 

Slope = 

y3  ‘  y2 

.  9 

Ax 

(Second 

Slope) 2  «  (l/R)2 

y  3  - y  2  •  y2  -  yi 

_  AX  *  AX 

Ax 

y3  “  ^2  +  yl 


(E.2) 


(E.3) 
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I 

k 

ft 


R 3-  P.I  Nofatioh  far  s-fra/yj  reacfi'njs 


Ax 

Ax 

Ax 

r  n 

F'r  £•/  SJo  fa  Hon  for  def  lection  s 


227 


APPENDIX  P 


The  singly  redundant  portal  frame  shown  in 
Pig,  P, 1  was  tested  beyond  the  elastic  limit.  In  order 
to  know  the  bending  moment  at  the  center  line,  the 
horizontal  reaction  must  be  measured  during  the  entire 
loading  cycle. 


The  bending  moment  at  the  centerline  section 
is  given  by: 

M  =  P  |  -  Hh  -  ,1 

*  p  (|£)  -  Hh  (P,1) 

8 

Substituting  the  values  of  £  =  11.2",  h  =  5.6" 

M  *  4.2  p-  5.6H  (P.2) 

The  values  of  H  were  taken  from  Pig.  4.1.21. 
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